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Abstract 

Since its commissioning in 1991, the ep collider HERA has been running suc- 
cessfully for almost a decade without stopping improving its performance. In 
this report, the inclusive cross section and structure function measurements for 
the deep inelastic scattering of neutral and charged current processes in the full 
HERA kinematic domain are reviewed. The results are compared with the Stan- 
dard Model expectations for the deep inelastic scattering processes. The new 
insights into the proton structure and on the underlying strong and electroweak 
interactions are discussed. 
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Chapter 1 
Introduction 



Scattering experiments have proven to be a very effective technique for probing 
the structure of matter. The large angle scattering of a particles (a few MeV) off 
a gold foil led Rutherford to conjecture around 1910 that the bulk of the mass 
of the atom must be contained in a very small, positively charged nucleus of 
the atom. The energy available in such a collision is related to a certain wave- 
length. The higher the energy, the shorter the wavelength and thus the smaller 
the distances which can be resolved. To investigate atomic nuclei with higher 
accuracy it is therefore necessary to bombard them with particles of higher en- 
ergy. In the 1950s and 1960s, starting with electron scattering experiments at 
Stanford |l|, experiments were performed at Darmstadt, Orsay, Yale, CEA and 
DESY. At the energies available, the experiments were however restricted to elas- 
tic scattering or excitation of the low-lying resonances. As higher energy electron 
beams (up to 16 GeV) became available at SLAG in the late 1960s, inelastic scat- 
tering experiments could be performed. These experiments 0, Q together with 
those at DESY |^ showed that at large four-momentum transfers, the inelas- 
tic nucleon structure functions were independent of any dimensional quantity, a 
phenomenon known as scaling. This result was interpreted as evidence for the 
existence of point-like constituents in the nucleons. These constituents are now 
known as quarks. Subsequent fixed-target lepton-nucleon scattering experiments^ 
at still higher energies have established the existence of violations of the scaling 
behavior j^, ^] which has been one of the most dramatic successes of perturbative 
Quantum Chromodynamics (QCD). 

Apart from the important roles deep inelastic scattering (DIS) has played 
in understanding nucleon structure, it also contributed to our understanding in 
other ways (see e.g. [0]). Neutral-current phenomena, which provided the first 

^The use of high energy leptons (electrons, muons and neutrinos, along with their antipar- 
ticles) as probes of nucleon structure is unique because the leptons do not interact strongly, 
so that they are able to penetrate the nuclear surface; their short wavelength implies that the 
leptons collide with individual charged or weakly interacting constituents Q. 
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demonstration of the SU{2) x U{1) unification of weak and electromagnetic inter- 
actions, were discovered and corroborated 0| in neutrino-nucleon experiments; 
interference between neutral weak and electromagnetic propagators was demon- 



strated first in an electron-nucleon inelastic scattering experiment |T^. The most 
precise measurement of the weak mixing angle 6w, parametrizing the weak neu- 
tral current coupling, used to be measured from neutrino-nucleon data. Today, 
the most precise DIS data are being widely used to determine the momentum dis- 
tributions of partons in hadrons allowing to predict cross sections in high energy 
hadron collisions. 

The fixed-target experiments, though precise, are limited in their kinematic 
range. The "Hadron Elekton Ring Anlage" (HERA), the first electron proton 
collider ever built (Sec.|2^), running at a center-of-mass energy of up to y/s ~ 
320 GeV, equivalent to a fixed-target experiment with a lepton beam of 50 TeV, 
has significantly extended the kinematic domain explored so far. The proton 
structure are being probed at 10 times smaller distances (down to ~ 10~^^cm, 
which is one per mill of the proton radius) than previously accessible. Partons 
can be studied down to very small fractional proton momenta, Bjorken x ~ 10~^. 
Two general purpose experiments, HI and ZEUS, are dedicated to the study of 
this physics. 

Already the first 0.02 pb~^ of luminosity delivered by HERA in 1992 has 
brought striking results which have opened new interest in QCD. Since then 
the luminosity has been increased by three orders of magnitude and the latest re- 
sults achieved are comparable in precision with fixed-target experiments but in a 
different kinematic regime. At high momentum transfers, based on an early data 
sample of 1994-1996, an excess of events with respect to the expectation of the 



standard DIS has been reported by both HI |TT| and ZEUS |T^. This has again 
initiated a considerable amount of theoretical interest in possible explanations 
within both the Standard Model(SM) and models beyond the SM [^. Results 
obtained using full data taken from 1994 to 2000 will be reviewed in this report. 
The emphasis will however be put on the experimental understanding of the de- 
tector performance by showing in some detail how some of the measurements 
have been carried out. 

The report is organized as follows. In Chapter ^ the experimental aspects 
concerning the HERA machine, the HI detector and the kinematic reconstruction 
are described. In Chapter ^ the theoretical framework is briefly reviewed. In 
Chapter ^, the structure function results at low momentum transfers and the 
interpretation within QCD are given. In Chapter ^, events at high momentum 
transfers are studied. Also discussed are their impact both within the SM for 
valence quark distributions and electroweak tests and beyond the SM for new 
physics search. A summary and an outlook are presented in Chapter |^. 



2 



Chapter 2 

Experimental setups and 
techniques 

Deep inelastic scattering (DIS) experiments can roughly be regrouped into three 
generations. The first generation of fixed-target experiments between 1968 and 
the mid-1970s has played an essential role in establishing the point-like sub- 
structure of the nucleons and in discovering the underlying weak neutral current 
interactions. The second generation of fixed-target experiments of higher energies 
and precision has helped in establishing and testing the quantum chromodynam- 
ics (QCD) and in providing data for extracting a first set of parton distribution 
functions. A list of selected DIS experiments is briefly presented in Sec. |2.2| after 
having introduced the kinematics in Sec. |2.1| . The third generation experiments 
at HERA ep collider are described in more detail in the later sections. 



2.1 Kinematics of deep inelastic scattering 

The kinematics of the inclusive DIS processes for neutral current (NC) interaction 
Ip IX and charged current (CC) interaction ep — > vX or vp IX (Fig. |2.1|) 
at the given center-of-mass energy squared, s = {k + P)^, is determined by two 
independent Lorentz invariant variables, which could be any of following variables: 

Q^ = -q^ = {k-k')\ g^e[0,s], (2.1) 
x=^, XG[0,1], (2.2) 

2/ = !^, ye [0,1], (2.3) 

W^ = {q + Pf, W^e[Ml,s\. (2.4) 

The square of the four momentum transfer (the invariant mass of the exchanged 
vector boson), (f < 0, is space-like and determines the hardness of the interaction. 
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/ (k) I (k') 




Figure 2.1: Schematic diagrams oflepton (I) and proton (p) scattering via photon 
and Z exchange for NC and W exchange for CC. The four momentum vectors of 
the particles or particle systems are given in parentheses. 



or in other words, the resolving power of the interaction. The Bjorken variable x is 
interpreted as the fraction of the proton momentum carried by the struck parton 
in the quark parton model (QPM). The variable y measures the inelasticity of 
the interaction and its distribution reflects the spin structure of the interaction 
in the rest frame of the target. 

The variable u 



is often used in fixed-target experiments where v = E — E' and M is the target 
mass. 
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2.2 The fixed- target experiments 



Fixed-target DIS experiments may be divided into two classes depending on the 
nature of the probe used, which in turn determines the force involved. In elec- 
troproduction, electrons or muons are scattered off the target nucleon and the 
force involved is electromagnetic. The leading process of the scattering is that of 
single-photon (7) exchange. The second class of processes is called neutrinopro- 
duction and, in this, neutrinos are scattered off the target nucleons by the weak 
force. The leading process is that of single-VT-boson exchange for charged current 
interactions and of single- Z°-boson exchange for neutral current interactions. 

The earlier electroproduction experiments used electrons as probe. However 
due to their large synchrotron radiation, the largest electron energy used for scat- 
tering experiments was limited to about 12GeV (Cornell) and 26 GeV (SLAC). 
The kinematic range was extended by muon-nucleon scattering to higher energies 
(20 < ^ < 500 GeV). 

The neutrinopro duct ion experiments complement the charged lepton scat- 
tering experiments and comparisons between neutrino and charged lepton deep 
inelastic experiments provide tests of universality of the parton density functions. 
In addition, the neutrinoproduction experiments are unique as they are able to 
distinguish quark flavors within target nucleons (Sec. |3.1.^ ). On the other hand. 



these experiments were in practice difficult to realize as the electrically neutral, 
weakly interacting neutrinos cannot be directed by electric and magnetic fields, 
as can electrons; and building a usable neutrino beam is a complicated process. 
First, a primary beam of protons is accelerated to high energy and is made to 
collide with a stationary target such as a piece of iron. From these collisions, a 
host of secondary particles, mainly mesons, will emerge in the general direction of 
the incident proton beam, although with somewhat smaller energies. These sec- 
ondary mesons can then decay into neutrinos or antineutrinos and various other 
particles by decays such as vr^ — > /x^ -|- i^^{T'^). Because the muon-decay mode of 
the mesons is generally the most common, it is mainly muon-type neutrinos which 
make up the beam. Finally, the neutrinos are isolated by guiding the secondary 
beam through a barrier of steel and rock equivalent of, perhaps, 0.5 km of earth. 
Only the weakly interacting neutrinos can pass through this amount of matter 
and so the beam emerging from the far side of the barrier is a pure neutrino 
beam with a typical intensity of about 10^ particles per second. Two types of 
high energy beams are commonly used: narrow band beams with a momentum 
and charge selection of the secondary beam, and wide band beams without such 
selection. The narrow band beams provide a measurement of the incident neu- 
trino energy and fiux. Furthermore, the beams are either almost purely neutrino 
or almost purely antineutrino, with little cross contamination. Wide band beams 
provide an order of magnitude more neutrinos; however, there is no direct check 
on the event neutrino energy and no direct measure of the fiux. 
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Table 2.1: A list of selected Gxed-target DIS experiments. Only the main 
targets used are indicated with p, d, and Fe standing respectively for pro- 
ton, deuterium, and iron targets (the use of heavy nuclear targets have the 
advantage of high rate). 



Experiment 


Year 


Reaction 


Beam Energy (GeV) 


Reference 


SLAC-MIT 


1968 


ep, ed 


4.5-20 


B |, EH 


DESY 


1969 


ep 


6 


Em 


Gargamelle* 


1971-1976 




< 10 


i 


SLAG 


1974-1975 


jj, Fe 


56.3, 150 


i 


BEBC* 


1975-1983 




< 200 




Fermilab 


1977-1978 


fip 


47, 96, 219 


i 


CDHS 


1979-1990 


v^,{V^) Fe 


< 200 


m 


1 


CCFR 


1979-1988 


v^,{V^,) Fe 


< 600 


m 


1 


BCDMS 


1981-1985 


jip, lid 


100,120,200,280 




1 


NMC 


1986-1989 


HP, fid 


90,120,200,280 


m 


1 


E665 


1987-1992 


fip, fid 


470 


m 


1 



* The Gargamelle and BEBC (Big European Bubble Chamber) exper- 
iments used bubble chambers while other cited experiments used elec- 
tronic detectors. The bubble chambers have opted for wide-band beams 
that maximize neutrino flux. 



A few selected fixed-target DIS experiments are listed in Table p.l| . The 
early fixed-target experiments between 1968 and the mid-1970s have discovered 
a wealth of information on the structure of the proton, to name a few: 

• Point-like constituents: The approximate independence of the measured 
structure functions of p|, H, H, [TH] indicated scattering off point-like con- 



stituents analogous to the classic Rutherford experiment on atomic struc- 
ture. 



Quark spin: The near vanishing of Fl in electron-nucleon scattering [0, 
[16| supported the assignment of half-integer spin for the quarks. 

Fractional charge: The comparison of structure functions in electron and 
neutrino scattering reactions supported the assignment of fractional charges 
to the quarks. 

Scaling violations and QCD: The muon scattering experiments at high 
energies provided first evidences for scaling violations 0. The scaling 
violations of the DIS data have helped in establishing the theory of QCD 
(Sec.pD. 
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Gluon: The momentum sum rules in both electron- and neutrino-proton 



scattering |Q suggested that quarks carry only about half of the total 
proton momentum. The other half was thought to be carried by neutral 
gluons, the quanta of the interquark field force. 

The neutrino scattering experiments have played other particular roles: 

• Discovery of neutral currents: The first observations of neutral cur- 
rents were reported by the Gargamelle Collaboration [Q. This was the first 
experimental success of the electroweak theory of the Standard Model. 

• Precise measurement of sin Ow- Before the advent of LEP and SLC, the 
electroweak mixing parameter sm6w was best measured by the neutrino 
scattering experiments with a first result (sin6'vy = 0.28 ±0.05) p5[ already 
in 1976 leading to a W mass value of 70 GeV well before the direct 
observation of the W boson in early 1980s. 

• First experimental hint of charm: Neutrino productions of opposite- 
sign dimuon in subprocesses + (ci, s) ^ fi~ + c + X and z7^ + 
((i, s) fi^ + c + X (the second muon arises from the semileptonic decay 
of the charmed hadrons emerged from the charged quark c) provided a first 



experimental hint of charm [27 



• Strange quark sea and evidence for an SU{3) asymmetric flavor 
sea: The subprocess ly^ + s /i"*" + c+X offered a unique probe to measure 
the strange component of the nuclear sea. The subprocess -|- {d, s) — >■ 
H~ + c + X provided a measure of the momentum fraction carried by the 
strange sea quarks relative to that carried by the non-strange sea s/(u + d) 
revealing an asymmetric strange sea. 

• Unique direct measurement of \Vcd\'- The Cabibbo suppressed subpro- 
cess + d ^ fi~ + c + X has been the only direct measurement of the 
Kobayashi-Maskawa matrix element \Vcd\ [ P^ . 

During the last decade, several precise data have been published and a few old 
data have been reanalyzed. A number of previously observed discrepancies be- 
tween different experiments have been resolved. For a recent review, see Ref. ||2l 
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2.3 The HERA accelerator 



The HERA project was authorized in April 1984 and the construction was com- 
pleted late 1990 in accordance with the original time schedule. The layout of 
the HERA accelerator and its preaccelerator facilities is shown schematically in 



Fig. p.2| . The electrons and protons are injected into HERA from PETRA (a pre- 
vious e^e~ collider) with 14 GeV and 40 GeV respectively. HERA is a double 
ring collider. The electron ring is made of supercoducting cavities and normal 




Figure 2.2: The layout of the HERA accelerator and its preaccelerator facilities. 



conducting magnets for accelerating the electrons of 14 GeV up to the nominal 
electron energy of 30 GeV. The electron ring provides either electron or positron 
beam|^. The year-dependent lepton-beam charges and energies are shown in Ta- 
ble p.2| . The actual maximum electron beam energy of 27.5 GeV is limited by 
the maximum available radio frequency (rf) voltage. The proton ring consists 
of conventional RF cavities and superconducting magnets of 4.68Tesla magnetic 
field for accelerating the protons from 40 GeV to 820 GeV, maintaining them at 
this energy and keeping them on orbit. The quality of the magnets is such that it 
is possible to rise the strength of the magnetic field up to 5.8Tesla, which would 
correspond to 1 TeV proton beam. 

HERA has four straight sections (interaction points) spaced evenly around 
its 6.3 km circumference. The electron and proton beams collide head-on in two 
of these points occupied by two general purpose detectors HI p9[ (Sec. |2.4|) and 
ZEUS pg. Two fixed-target experiments HERMES and HERA-B [0 are 



located at the other two interaction regions. At HERMES spin rotator provide 
the longitudinally polarized electron beam in collision with a polarized gas target 



^In the following, the generic name electron is used for electrons as well as for positrons 
unless stated otherwise. 
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Table 2.2: The year-dependent lepton-beam charges and 
energies and the corresponding center- of- mass energy y/s 



of the HERA collider. 



Year 


beam 


Ee (GeV) 


Ep (GeV) 


(GeV) 


1992-93 




26.7 


820 


296 


1994 


e+ 


27.5 


820 


300 


1995-97 


e+ 


27.5 


820 


300 


1998 




27.5 


920 


320 


1999 


e+ 


27.5 


920 


320 


2000 


e+ 


27.5 


920 


320 



Table 2.3: The main HERA parameters 
achieved so far compared with the design 
values. 



Parameter 


Design 


Achieved 


le (mA) 


58 


51.8 


Ip (mA) 


160 


109 


jj bunches 


210 


189 


cr^ at IP (/im) 


280 


179 


cr^ at IP (/im) 


50 


48 


aP at IP (cm) 


11 


11 


Ci^st (10^1 cm-2s-i) 


1.5 


1.8 



for studying spin dependent structure functions, while HERA-B uses the halo of 
the proton beam with wire targets in an attempt to detecte CP violation in the 
B system. 

HERA has made steady progress since 1992. This can be seen from Fig. |2.3| 
in which the time evolution of the following information is shown: (1) the peak 
luminosity, (2) the integrated luminosity collected by HI, (3) the mean as well 
as maximum lepton and proton beam currents, and (4) the number of collid- 
ing bunches. A comparison of the best achieved and design values of the main 



HERA parameters is shown in Table p.3| . The improvement of the luminosity as 



a function of time can be understood from the formula: 



Coc^. (2.6) 



The improvement in 1993 with respect to 1992 was mainly due to the increases 
in the beam currents; from 0.94 to 10.8 mA for protons and from 1.33 to 7.7mA 
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7992 799J 7994 7995 7996 7997 799S 7999 2000 

Year 



Figure 2.3: The time evolution of the following information: (1) the peak lu- 
minosity, (2) the integrated luminosity collected by HI, (3) the mean as well 
as maximum lepton and proton beam currents, and (4 ) the number of colliding 
bunches. 

for electrons. This was achieved mainly by increasing the number of colliding 
bunches from 9 in 1992 up to 84 in 1993. During the summer of 1994, the elec- 
tron beam was replaced by a positron beam which has considerably improved the 
beam lifetime^ at high currents. The number of colliding bunches have increased 
subsequently to 153 in 1994 and then to about 175 since 1995 with a maximum 
number of 189. In addition to these colliding bunches, some bunches are left un- 
paired (so-called pilot bunches, i.e. the corresponding bunch in the other beam 

^The electron lifetime was severely limited due to presumably its interaction with positively 
ionized impurities in the beam pipe. To cure the problem, the original ion getter pumps of the 
electron ring have been replaced in the 1997/1998 shutdown by passive non-evaporating getter 
pumps (adsorption pumps without high voltage that do not accelerate dust particles into the 
beam vacuum). 
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is empty). These pilot bunches are used for evaluating the beam related back- 
grounds both for the luminosity determination and for physics analyses. Two 
successive bunches are separated in time by 96 ns. This is to be compared with 
the collision frequency of 22 fis and 25 ns respectively at LEPQ and LHC 

Other improvements include 

• the full exploitation in the aperture margin of the machine in order to 
squeeze the beam cross sections at the collision point down to 179 /xm x 
48 /im, which are two times smaller than the design value, 

• significant progress in the proton beam intensity by improving controls and 
beam handling in the injector chain. The maximum proton beam current 
achieved in ep collisions is now 109mA (or ~ 8 x 10^° protons per bunch). 

The maximum electron current obtained is 51.8 mA. The beam current is limited 
by the available RF power. Therefore, although the beam currents still fall short of 
expectations, the maximum peak luminosity achieved 1.8 x 10^^ cm~^s~^ exceeds 
the design goal of 1.5 x 10'^^ cm"^s~^ by focusing very tightly the beams at the 
interaction points. 

An important ungrading program has been prepared and has started since 
September 2000. After the upgrade, an increase in the luminosity of more than 
a factor of 5 is expected P3[ . 



■^The number corresponds to four bunches per beam during the early (late) run of 1989- 
1992 (1996-2000), the collision frequency was reduced by a factor of two or more by increasing 
the number of bunches for the years between 1992 and 1995. 
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2.4 The HI detector 



The HI detector ||29| (Fig. p.4|) is nearly hermetic multi-purpose apparatus built 
to investigate ep interactions at HERA. The main components are the tracking 



THE H1 DETECTOR 

Forward Muon System Central Muon Sj'stem Forward Tracker 




Figure 2.4: An overview of the HI detector. 



detectors (See p. 4.1)) , the calorimetry (Sec. |2.4.^ ), the superconducting magnet 



the muon detectors (Sec. p. 4. 3D , the very forward detectors (Sec. p. 4. 4] ) and the 



luminosity detectors and various electron taggers (Sec. |2.4.5| ). In many aspects 



the HI detector^ does not differ strongly from the detectors at e~^e~ or pp colliders. 
Specific at HERA is however the imbalance in the energy of the two colliding 
beams, which requires an asymmetric detector. 

*The ZEUS detector pOl is rather similar to the HI detector in many aspects. The main 



difference hes in the main calorimeter; as described in Sec. 2.4.2, the liquid-argon calorimeter 
of HI has a fine granularity and a good resolution for electromagnetic objects whereas the 
uranium-scintillator calorimeter of ZEUS is compensating, giving equal response to hadrons 
and electrons, and has a good resolution for hadrons {a{E)/E ~ 35%/ Ve 2%). The ZEUS 
calorimeter has in addition a good time resolution, which is better than 1 ns for energy deposits 
greater than 4.5 GeV thus providing a fast time information for triggers and background rejec- 
tions. The location of the superconducting coil providing magnetic field for trackers is different 
as well; the coil of HI is placed outside of the calorimeter thus minimizing the amount of dead 
material in front of the calorimeter, while ZEUS had placed it between the calorimeter and the 
tracker. 
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The coordinate system convention for the experiment defines x pointing to 
the center of the HERA ring, y the upward direction, and the forward, positive z 
direction as being that of the proton beam. The polar angle 9 is defined relative 
to this axis such that pseudo-rapidity, 77 = — In tan 6*72, is positive in the forward 
region. 

2.4.1 The tracking detectors 

The tracking detectors consist of central jet chambers (CJCl, CJC2), central 
trackers for measuring the z coordinate (CIZ, COZ), central multiwire propor- 
tional chambers for fast triggering (CIP, COP), forward tracking detector (FTD), 
backward tracking detector (BPC, BDC), and central and backward silicon mi- 
crovertex detectors (CST, BST). 

The CJCl and CJC2, covering a polar angular range from 15° to 165°, are 
two large, concentric drift chambers. The inner chamber, CJCl, has 24 layers of 
sense wires arranged in 30 phi cells, while CJC2 has 32 layers of sense wires in 60 
phi cells. The cells are at a 30° angle to the radial direction. The point resolution 
is 170 /xm in the r — direction. The z coordinate is measured by charge division 
and has an accuracy of 22 mm. A superconducting solenoid, which surrounds 
both the tracking system and the liquid argon (LAr) calorimeter, provides a 
uniform magnetic field of 1.15Tesla. The momentum of charged particles may 
be determined from their track curvature in the magnetic field with a transverse 
momentum resolution of (Tpj./pT < O.OlpyCeV. The dE/dx resolution for a well 
measured track is better than 7%. 

Two thin cylindrical drift chambers CIZ and COZ have sense wires perpendic- 
ular to the beam axis, and therefore complement the accurate r — measurement 
provided by the CJC with precise z coordinates. The CIZ is located at 18 cm 
inside the CJCl, while COZ is located at 47cm between CJCl and CJC2. These 
two chambers deliver track elements with typically 300 yum resolution in z. To 
each of the z chambers a proportional chambers (CIP/COP) is attached for trig- 
gering. 

The FTD, covering an angular range from 7° to 25°, are integrated assemblies 
of three supermodules, each including, in order of increasing z, three different 
orientations of planar wire drift chambers (each rotated by 60° to each other in 
azimuth), a multiwire proportional chamber (FPC), a transition radiation detec- 
tor and a radial wire drift chamber. 

The backward proportional chamber BPC, located just in front on the back- 
ward calorimeter, is made of four planes of wires with vertical, horizontal and 
±45° orientations. The wires are strung every 2.5 mm, and signals from two 
wires are fed to one preamplifier. Three out of four planes are required in coin- 
cidence in order to reconstruct a space point with a spatial resolution of about 
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1.5 mm in the transverse plane. The BPC provided an angular measurement of 
the electron in the range from 155° to 174°, together with the vertex given by 
the main tracking detectors with a precision better than 1 mrad. This detector 
has been replaced in the 1994/95 shutdown by an eight layer drift chamber BDC 
with an extended acceptance between 155.1° and 177.5°. 

2.4.2 Calorimetry 

The calorimetry system consists of the LAr calorimeter, the backward electro- 
magnetic calorimeter (BEMC/SPACAL), the forward calorimeter (PLUG) and 
the outer calorimeter, the so-called tail-catcher. 

The emphasis is put on the electron recognition and energy measurement. 
This led to placing the LAr inside the magnetic coil in order to minimize the 
passive material. LAr was chosen because of its good stability, ease of calibration, 
possibility of fine granularity and homogeneity of response. The LAr calorimeter 
covers the polar angle range between 3° and 155°. The calorimeter along the 
beam axis is segmented into 8 "wheels" , with each wheel being further segmented 
into octants in (Fig. |2.5| ). The structure of the electromagnetic section (EMC) 
consists of a pile of GlO-Pb-GlO sandwiches separated by spacers defining the 
LAr gaps. The hadronic section (HAC) is made of stainless steel absorber plates 
with independent readout cells inserted between the plates. The orientation of 
the plates varies with z such that particles always impact with angles greater 
than 45°. The granularity ranges from 10 — > 100 cm^ in the EMC section and 
to 50 — > 2000 cm^ in the HAC section. Longitudinal segmentation is 3-4 layers 
in the EMC over 20-30 radiation lengths (^o) and 4-6 layers in the HAC. The 
total depth of both sections varies between 4.5 and 8 interaction lengths (A). 
The most backward part of the LAr calorimeter is a smaller electromagnetic 
calorimeter (BBE) which covers the polar angle ranging from 146° to 155°. The 
LAr calorimeter has a total of 45 000 readout cells. The noise per cell ranges from 
10 to 30MeV. The resolution measured in the test beam is 0.12/^E(GeV) ©0.01 
for an electromagnetic shower and 0.5/ \/ -E'(GeV)©0.02 for a hadronic shower (the 
symbol © denotes addition in quadrature). The electromagnetic and hadronic 
energy scale uncertainty is respectively 1-3% and 2-4%. The calorimeter is non- 
compensating, with the response to hadrons about 30% lower than the response 
to electron of the same energy. An offline weighting technique is used to equalize 
the response and provide the optimal energy resolution P^ , . 

In the backward region, the LAr was complemented in the angular range from 
151° to 176° by the BEMC, a conventional lead-scintillator sandwich calorimeter. 
The BEMC was replaced in the 1994/95 shutdown by a lead/scintillating fiber 
calorimeter (SPACAL). The calorimeter BEMC is made of 88 stacks with a size 
of 16 X 16 cm^ and has a depth of 21.7Xo, or approximately 1 A, which on average 
contains 45% of the energy of a hadronic shower. The four inner stacks around the 
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beam pipe are of triangle shape. A 1.5 cm spatial resolution of the lateral shower 
position is achieved using four photodiodes which detect the wavelength shifted 
light from each of the scintillator stacks. The electromagnetic energy resolution is 
a/E = 10%/^E(GeV) ©1.7%. A scintillator hodoscope (TOF) situated behind 
the BEMC is used to veto proton induced background events based on their early 
time of arrival compared with nominal ep collision. The new calorimeter SPACAL 
has both electromagnetic and hadronic sections. The electromagnetic energy 
resolution is 7 .5% / ^^/^{GeV) (B 2.5%, and the hadronic section has 2 A and an 
integrated timing function to veto proton beam induced background interactions 
(the old TOF system could thus be removed). The angular region covered is 
extended (from 153° to 177.5°) compared to the BEMC, and the calorimeter has 
a very high granularity (1192 cells) yielding a spatial resolution of about 2 mm. 

The LAr and BEMC/SPACAL calorimeters are surrounded by the iron re- 
turn yoke, which is instrumented with 16 layers of limited streamer tubes (LST). 
Eleven of the 16 layers are equipped with readout electrodes (pads). From ion- 
ization energies of particles passing through the chambers the energy of tails of 
hadronic showers leaking out of the LAr are measured in the analog readout sys- 
tem. Therefore the system is also called tail-catcher (TC). The energy resolution 
is a/E ^ 100%/v/^(GeV). 

In the forward direction around the beam pipe, the angular range from 0.3° 
to 3.3° is covered by the PLUG, a sampling calorimeter consisting of nine cop- 
per absorber plates interleaved with eight sensitive layers of large area silicon 
detectors. 

2.4.3 Muon detectors 

Recognition of muons is very important in the study of heavy quarks, heavy 
vector mesons, W^-production and in the search for exotic physics. 

Muons in the central region are identified by looking for particles penetrating 
the calorimeter and coil and leaving signals in the TC. Three of the 16 instru- 
mented LST layers are located before the first iron plate, and three after the last 
iron plate. There is a double layer after four iron plates, and eight single layers in 
the remaining gaps between the iron sheets. A minimum muon energy of 1.2 GeV 
is needed to reach the first LST, while 2 GeV muons just penetrate the iron. 

In the very forward direction, a spectrometer composed of drift chambers sur- 
rounding a toroidal magnet with a field of 1.6Tesla is used to measure muons. 
This spectrometer measures muons in the momentum range between 5 to 200 GeV. 
The lower limit is determined by the amount of material traversed, while beyond 
the upper limit the muon charge can no longer be measured unambiguously. 
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2.4.4 Very forward detectors 

HI has spectrometers downstream of the main detectors in the proton beam 
direction to measure high energy protons, as well as calorimeters at zero degrees 
to measure high energy neutrons. These are used in the study of diffractive 
scattering as well as in the study of leading-particle production. 

A forward proton spectrometer (FPS) has been installed since 1995 at 81 and 
90m away from the interaction point, which detects leading protons in the mo- 
mentum range from 580 to 740 GeV and scattering angles below 1 mrad. The FPS 
has been extended with stations at 80 and 63 m since 1997. In all stations the pro- 
tons are detected with scintillating fiber hodoscopes. The detector elements are 
mounted inside plunger vessels, so called Roman Pots, which are retracted during 
injection and are brought close to the beam after stable luminosity conditions are 
reached. 

The forward neutron calorimeter (FNC) is located at 107 m downstream of the 
interaction point. The calorimeter consists of interleaved layers of lead and scin- 
tillating fibers. The calorimeter has a total depth of 9.5 A and has an acceptance 
> 90% for neutrons with a production angle below 1 mrad. 



2.4.5 Luminosity detector and electron taggers 

At HERA, the luminosity is determined from the rate of the bremsstrahlung 
(Bethe-Heitler) process ep ep'j (Sec. p.4.2|) using a small luminosity system 



(the example from HI is indicated schematically in Fig. 2^ and shown in a larger 



view in Fig. |2.6|) . The process has a large and precisely known cross section. In 




p 
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Figure 2.6: The layout of the HI luminosity system consisting of an electron 
tagger (ET=etag) and a photon detector (PD). Note that the z axis coincides 
with the proton beam direction. 

the HI area, the final state electron, deflected by a set of low-beta quadruples and 
a bending magnet located in the region —23.8 m < z < —5.8 m, passes an exit 
window at z = —27.3 m and hits the electron tagger (etag) a.t z = —33.4 m. The 
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etag is located beside the electron beam pipe and is made of 49 crystals covering 
a total area of 154 x 154 mm^. The photon leaves the proton beam pipe through 
an exit window at z = —92.3 m where the proton beam pipe bends upwards. The 
photon detector (PD), situated at z = —102.9m on the z axis, is built out of 
25 crystals with a total surface of 100 x 100 mm^ and a depth of 22 radiation 
lengths (^o)- A 2Xo Pb filter followed by a IXo water Cerenkov veto counter 
(VC) in front of the PD protects the detector from the high synchrotron radiation 
flux. The VC eliminates events with photons interacting in the filter. The whole 
system provides a fast relative luminosity measurement with a statistical precision 
of about 2% at nominal beam conditions. 

The luminosity system serves in addition several other purposes. It provides 

• electron beam monitoring for the HERA machine, 

• absolute luminosity measurement in the interaction region with an accuracy 
of better than 5%, 

• tagging of photoproduction events, 

• energy measurement for electrons scattered under small angles and 

• energy measurement for photons from initial state radiation (Sec. O). 



The etag at 33 m measures electrons in the energy range 0.25i?e < E'^ < 
0.75Ee with an average acceptance of 48%. The acceptance to higher energy 
(0.76£'e < -Eg < 0.96-E'e) is covered by a new tagger at 44 m installed since 1995. 
The lower energy is covered by yet another tagger at 4 m. 



2.4.6 Trigger and event reconstruction 

The purpose of the trigger system is a fast separation of the interesting physics 
events from background events. The main background sources, common as those 
presented to other accelerator experiments, are synchrotron radiation from the 
electron beam, proton gas interaction in the beam pipe vacuum of about 10~^ hPa 
and stray protons, which produce particle showers by hitting the beam tube 
and other materials around the accelerator (beam-gas and beam- wall). Beam 
halo muons and muons from cosmic radiation also contribute. The rate of the 
background is up to 10^ times higher than the rate of a variety of physics processes 
under study in ep collisions (Table |2.4|) . The rate of physics processes extends 
from photoproduction, where the visible ep cross section of several fih implies an 



event rate of 20-30 Hz at design luminosity (see Table p.3|) , towards W production 
expected to occur a few times per week. 

The trigger system is based on 4 levels in order to filter the intersting physics 
events, followed by the offline reconstruction of the kept events. 
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Table 2.4: A comparison of cross sections and/ or rates between the 
main background sources and various physics processes at design 
luminosity. 



Source/process 


Cross section 


Rate 


Beam- gas (wall) interaction 




50 kHz 


Cosmic fi in barrel 




700 Hz 


Tagged 


1.6 fih 


25 Hz 


cc total 


1 fih 


15 Hz 


DIS low g2 


150 nb 


2.2 Hz 


DIS high Q2 (e in LAr) 


1.5 nb 


1.4min~^ 


Charged current DIS {pt > 25 GeV) 


50 pb 


3h-i 


W production 


0.4 pb 


0.5 d^^ 



The trigger level 1 (LI) system makes a decision within 2 yus on whether to 
accept or to reject an event using information provided by different subdetec- 
tors ("trigger elements (TE)"). The central trigger logic (CTL) combines these 
trigger elements into 128 "subtrigger (ST)". Not all subdetectors can provide 
this information fast enough to make a decision after each bunch crossing (BC) 
immediately. Therefore the information is sent into pipelines where it is kept 
until all relevant subdetectors have provided their trigger elements. The delay 
of 2 fis (24 BC) is necessary as some subdetectors are relatively slow: the CJC 
takes 11 BC due to the longest drift time of 1 /is and the LAr takes 13 BC due to 
long integration time of the preamplifier^. If any of the ST conditions is fulfilled 
by the event, the pipeline is stopped immediately and the signal is passed to the 
next trigger level^ 

The level 2 (L2) and level 3 (L3)[] triggers operate during the primary dead 
time of the readout of about 1.5ms. They work on the same data as the LI, 
and reach a decision within 20 /is and 800 /is, respectively. While at the LI, only 
minimum correlation between different subdetectors (e.g. between the MWPC 
and LAr) is used, the L2 and L3 make full use of the detailed, high granularity 
trigger data of most subsystems. The L2 system includes a complex topological 
correlation and a neural network approach. 

The level 4 (L4) trigger is based on full event reconstruction in MIPS R3000 
based processor boards. Algorithms similar to the ones used for the offline analysis 
are used to select valid events. The events accepted by L4 are written to tapes 

^The typical drift time of an ionized electron in the LAr gap is about 200ns/mm, a gap of 
2.35mm between two absorber plates results already in about 500ns, i.e. about SBC. 

^Some subtriggers are prescaled such that not every event fulfilling the subtrigger conditions 
is kept. 

^This trigger has not yet been used in HI. 
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with a rate of about 15 Hz. 

The accepted raw data or those simulated Monte Carlo data are then fully 
reconstructed and assigned into different physics event classes. The data infor- 
mation is written in a compressed format to Data Summary Tapes (DST) which 
are the basis for physics analyses. 



2.5 Monte Carlo technique and detector simu- 
lation 

Monte Carlo (MC) technique has proven to be indispensable for the extraction of 
physical quantities from the measurements. For the cross section measurement, 
MC programs can be used to determine corrections for acceptance, efficiencies, 
background contamination, and resolution effects of the detector system. Some 
of these corrections can be obtained directly from the data though often with 
a limited statistical precision. The MC, which provides in principle unlimited 
event samples, allows therefore to model the data with a better precision. The 
simulated MC events are also very useful in defining variables and tuning cuts 
for selecting signal events from various background contributions. 

More information concerning Monte Carlo generators used in different analy- 
ses will be given later together with the analyses. Once an event (either a signal 
or a background event) is generated, the HI detector response to the particles 
generated in the event is simulated in detail using the HI simulation package. 



HlsiM, which makes use of the geant program The parameters used by 

this program were determined in test beam measurements and optimized during 
ep data taking. For the simulation of the energy response of the calorimeter a 
fast parametrization is used for the development of electromagnetic and hadronic 
showers to save computing time. These simulated events are then subject to the 



same reconstruction program (Sec. |2.4.^ ) as the data and the same analysis chain. 
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2.6 Kinematics reconstruction and coverage 



2.6.1 Reconstruction of kinematic variables 

For NC events, the kinematics is over-constrained as the HERA experiments 
measure both the scattered electron and the hadronic final state. Here are a few 
commonly used methods: 

Electron method uses the energy, E'^, and the polar angle. Of,, of the scattered 
electron measured relative to the proton beam direction: 

= 4E,El cos' ( I") (2.7) 



2 

Ee'"" \2 

Xe = ^. (2.9) 

sye 

Since the incident beam energy appears both in Ql and in y^, this method 
is thus sensitive to the initial state radiation (the effective beam energy after 
radiation is smaller than Ee, see Sec. p.4| , and Eg. ( [4. 14] ) in Sec. [4. 4]) . 



Double angle (DA) method [37] uses the electron polar angle, 6e, and the 



inclusive hadronic polar angle, 9h, which is the polar angle of the scattered 
quark in the QPM (Sec. |3.1.^ ) with massless quarks 

QIa = ^ (2.10) 

EpOeiae + ah) 

yDA = ^^, (2.11) 

+ Oih 



Q 



XDA = ^ , (2.12) 
syuA 



with 



tte = tan 




V 2 ; pt, 

= K- Vzfi , 



where the summations are over all particles of the hadronic final state. The 
DA method is also sensitive to photon emission of the primary electron. 
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On the other hand, it is, to a good approximation, insensitive to the energy 
scale uncertainties. It has thus been widely used in various analyses to 
check and improve the energy scales of the scattered electron and of the 
hadronic system. 



Sigma (E) method [38| is constructed such that no electron beam energy Eg 

Ql = ' , (2.16) 



is used directly in and y: 

E? sin^ 



1 - 2/s 

y^ = ^^^. (2.17) 

XT = ^. (2.18) 
This method is thus less sensitive to the radiative effects. 

The fact that different methods can be used to reconstruct the event kinematics 
has two consequences: 

• The kinematic domain is fully explored; for example the electron method 
is good at high y while the S method can be used at sufficiently low y. 

• Important cross-checks of systematic effects are possible between different 
methods. 

For a CC event, however, the kinematics can only be reconstructed from the 
hadronic final state: 

Ql = P^, (2.19) 
1 - l/h 

= ^ , (2.20) 

XH = ^. (2.21) 

sVh 

The understanding of the hadronic energy scale and its scale uncertainty is thus 
crucial for all CC measurements. 



2.6.2 Coverage of the kinematic phase space 

The HERA measurements in DIS cover a vast kinematic range as shown in 
Fig. p.T) . The HERA measurements extend those of fixed-target experiments by 
more than 2 orders of magnitudes in both x and Q^. At the highest values of Q^, 
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the measurements are still statistically limited (Chapter The measurements 
for ^ 4 GeV^ have been made possible by using different techniques and by 
upgrading experimental apparatus: 

• Using data in which the interaction point is shifted from its nominal position 
to the forward direction by ~ 70cm (so-called shifted vertex events). For 
the same backward calorimeter, the angular acceptance is increased when 
the interaction vertex is moved away from the nominal interaction point. 

• Using data with hard initial state radiation. A large fraction (~ 30%) of the 
photons from events with hard initial state radiation are measured in the 
photon calorimeters of the luminosity system (Sec.^.4.5D. The resultant ep 



collisions therefore occur at lower center-of-mass energies, and lower and 
higher x values can be attained. As an example, the HI analysis is shown 



in some detail in Sec. 4.4 to illustrate how the measurement is realized in an 



extended kinematic region. 

Upgraded apparatus. The angular acceptance of the scattered electron has 
been significantly extended in the shutdown period 1994-95 by replacing 
the old backward calorimeter and proportional chamber, BEMC and BPC, 
with the new calorimeter and drift chamber, SPACAL and BDC (Sec. p.4| ). 
The ZEUS experiment has added a small angle electron calorimeter, the 
beam pipe calorimeter, in the same shutdown period. 
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Figure 2.7: Coverage in the kinematic plane x — Q^ of the HERA experiments, HI 
and ZEUS, compared with that of the fixed-target experiments from BCDMS |2l|/ , 
NMC J2|/, E665 and SLAC /||/. The corresponding scale (A) probed by a 
virtual boson is also indicated. The points indicate the location where cross 
sections or structure functions have been measured. 
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Chapter 3 



Theoretical framework 



3.1 Deep inelastic scattering and the quark par- 
ton model 

3.1.1 Cross sections and structure functions 

The neutral current cross section for charged leptons scattering off nucleons 
l^N l^X studied in all fixed-target experiments and in the low part of 
the HERA kinematic domain arises to a good approximation from the photon 
exchange only. Taking the proton as an example, the cross section is directly 
proportional to the amplitude squared jiOl 



with 

A=\j,J'' (3.2) 

where the term 1/g^ corresponds to the propagator of the exchanged photon, 
and are respectively the lepton and proton currents. In Eg. ( |3.1|) , L^^, and 
W^" are respectively the leptonic and hadronic tensors. The leptonic part is 
completely determined by QED, whereas the hadronic one contains all the infor- 
mation about the interaction of the electromagnetic current with the proton 
target. By writing down the most-general possible combinations of all the mo- 
menta appearing in the reaction in L^j, and W^'^ and then simplifying the result 
using general theoretical assumptions such as parity and time-reversal invariance, 
the cross section can finally be expressed as 



dxdy Q'^ 



M'^xy \ y 



2 



(3.3) 
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where a is the fine structure constant, M is the proton target mass, and -^2(2;, Q^) 
and 2xFi{x,Q'^) are so-called structure functions, to be determined by the deep 
inelastic experiments (Chapter There are two structure functions because 
there are two essential degrees of freedom. The exchanged virtual photon can 
be transversely or longitudinally polarized. If one defines a^^l' as the total 7*p 
cross section, then it can be decomposed into the longitudinal part ctl and the 
transverse one ar as 

(TL + ar, (3.4) 



with 



a. 



IP 

tot 



xMK 
xMK^ 



2xFi 



(3.5) 
(3.6) 



For virtual photon, the fiux factor K cannot be unambiguously defined. The 
longitudinal structure function Fl is a combination of F2 and 2xFi: 



1 + 



2xFi . 



(3.7) 



The measurement is carried out either directly on F^ or on R, the cross-section 
ratio of the longitudinally polarized over transversely polarized virtual photon: 



R 



1 



AM^x"' 



2xFi 



- 1 



(3.8) 



so that the arbitrary fiux factor K does not enter. 

At HERA, high regime (Q^ > M|) is reached, where the contribution 
from pure Z° exchange and its interference with 7 exchange becomes increasely 
important. As the weak neutral current does not conserve parity, an additional 
term {xF^lXjQ"^)) has to be taken into account in Eq.( 



dxdQ"^ 



2 / 2 

[l - y)F2+^-2xF,T [y - \ ] xF^ 



Y+F2 - v^Fl T Y.xFs 



(3.9) 



where Y± = 1 ± {1 — y)"^ are helicity functions and the term proportional to 
has been neglected. Note that the structure functions F2{x,Q^), 2xFi{x,Q'^), 
and Fl[x,Q'^) have been replaced respectively by F2{x,Q'^), 2xFi(x, Q^), and 
Fl{x,Q'^) to account for the contributions from the interference and Z° 
exchange terms (Eqs.(|3A|),(g)). 
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The cross section for charged current processes, 



(Fixed-target experiments) (3.10) 
e p VjX /jjERA experiments) (3.11) 



can be derived in a similar way and reads 

rfV^c _ GlM^{s-M^) ( 1 ^' 



dxctQ"^ 2TIXS \Q'^ + Ml, 



X (3.12) 



where is the mass of the exchanged W boson. Again since the weak force 
does not respect parity invariance, a similar but different weak structure function 
xF^^i^x.Q'^) has to be introduced (compare Eqs.(CT) and (|3:27|),(^)). 



3.1.2 Structure functions in the quark parton model 

Two ideas, both put forward in 1969, have played an important role in the de- 
velopment of the experiments and in our understanding of them. The two ideas 
are those of the parton model by Feynman and of scaling by Bjorken ^ . 

The parton model is simply a formal statement of the notion that the nucleon 
is made up of smaller constituents, the partons. No initial assumptions about the 
partons are necessary as it is the purpose of the experiments to determine their 
nature. 

The scaling prediction states that when the momentum carried by the probe 
becomes very large, the dependence of the cross section on parameters such as 
the energy u = E — E' and momentum squared g^, transferred by the photon, 
becomes simple. In the parton model, the onset of this simple scattering behavior 
has a straightforward interpretation. The complicated scattering of the probe off 
a nucleon of finite spatial extent has been replaced by the scattering of the probe 
off a point-like parton. The photon ceases to scatter off the nucleon as a coherent 
object and, instead, scatter off the individual point-like partons incoherently. 

When it became clear that the hypothesized point-like constituents, partons 
had properties characteristic of quarks, the parton model came to be called the 
Quark Parton Model (QPM). 

While the cross sections shown in Eqs.( p.3| ),( |3.9| ), (|3.12| ) make no assumptions 
about the underlying structure of the hadron involved in the interaction, in the 
QPM, the structure functions F2 ^, xFg ^ can be expressed in terms of parton 
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distribution functions: 

Ff ^(x, Q') = M>^, Q')^ hi^) + U^)] , (3.13) 

i 

xFf^(x,g2) = J2B^i\.Q')x[q^{x)-q,{^)] , (3.14) 

i 

Fl = F2-2xFi = 0, (3.15) 

where the sum runs over all quark flavors with xq{x) specifying the probability 
of flnding a parton q carrying a momentum fraction x of the proton's momentum 
in a frame where the proton's momentum is large. The last relation is so called 



Callan-Gross relation and is a direct consequence of quark spin 1/2 due to 
the fact that a spin- 1/2 quark cannot absorb a longitudinally polarized vector 
boson. The couplings of the fermions to the currents depend on the lepton-beam 
polarization A: 

A{X, Q') = ^-^Afm + ^^AfiQ') , (3.16) 
5.(A, Q') = ^-^B^iQ') + (Q^) , (3.17) 

with 

Af'^(Q2) = - 2e,ivi ± ai)v,Pz + {m ± aifiv^ + a2)P| , (3.18) 
5f '^(g^) = ^2e,ivi ± ai)a,Pz ± 2{vi ± ai)\a,Pl , (3.19) 

for l~p scattering, the corresponding structure functions for l~^p scattering are 
obtained by swapping L i? L in Eqs.( p.l8|) , (|3.19|) . The vector and 

axial-vector couplings of the fermions are given by 

Vf = Tsf -2efsm'^9w , af = T^f, (3.20) 

where the deflnition holds for both changed leptons and quarks, T-^f is the third 
component of weak isospin, and 9\y is the Weinberg angle. The Z° propagator 
appears in the quantity Pz as: 

1 

^ Q'^ + Ml sin(2^vi/) ■ ^ ' ' 



For unpolarized lepton beams, Eqs.( |3.13| ) and ( |3.14| ) can be expressed^ as 



F2 = F2 - viPzFf + {v^ + )P|Ff ^ F2 + a^iP^ (3.22) 
xFg = -aiPzxFf + 2aiViPlxFi ^ -aiPzxPf (3.23) 



^Here we have taken the convention that Vp- = — 1/2 + 2siT? 9w and ae+ = ag- 

-1/2. 
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with F2, F^^ and 

F2, Ff, Fij = ^ [el 2e,v„ v] + a^] x (g, + (3.24) 

i 

standing for contributions respectively from pure photon exchange, interfer- 
ence and pure Z° exchange. The terms xF^^ and xF.^ 

xFf, xFi] = J2 [2e^a^, 2iHai] x (g, - (3.25) 



have similar meanings. The dominant correction to F2 arises from the ex- 
change term (the F2 term in Eq.( |3.22|) ) since vi is small with respect to ai. For 
the same reason, the dominant contribution to xF^ is from the ^Z^ term (the 
xF^^ term in Eq.( p.23| )). As shown in Fig. |3.1| , the generalized structure function 
F2 is always larger than the electromagnetic structure function F2 and does not 
depend on the charge of the incident lepton beam, whereas the structure function 
xFs has a negative (positive) interference for a positron (electron) beam and its 
contribution is large than that from the F2 term. As a result, the full cross 
section is expected to be smaller (larger) than that from pure photon exchange 
in the HERA kinematic range. 

The structure functions for charged current processes V^p — s> fi~^X, e~p — > u^X 
with left handed polarization are defined as 

F^^{x) = 2xF^^{x) = 2x [u{x) + c{x) + d{x) + s{x)] , (3.26) 
xF^^{x) = 2x [u{x) + c{x) - d{x) - s{x)] , (3.27) 

and F^'^ = 2xF^'^ = Fg^*^ = for e~p u^X with right handed polarization. 
The corresponding structure functions for processes u^p fJ'~X, e^p V(.X 
with right handed polarization are given by 

F^^{x) = 2xF^^{x) = 2x [d{x) + s{x) + u(x) + c{x)] , (3.28) 
xF^^{x) = 2x [d{x) + s{x) - u{x) - c{x)] , (3.29) 

and F2^ = 2xF^'^ = F^^ = for e^p V^X with left handed polarization. 
Comparing Eqs.( |3.2(^ ),( ^.28D with Eq. ( |3.13| ), one notices an important feature 



of the charged current interactions: they are able to distinguish between the 
partons and the antipartons of the target nucleon. This is because the space- 
time structure of the weak interaction ensures that target partons of differing 
helicities are affected differently. In the relativistic limit, in which the rest mass 
of a particle is regarded as being negligible, the parton and antiparton helicities 
are opposite, so they will interact with the W^-boson probe differently. Also, 
because the H^-boson probe is electrically charged, the target parton must be able 
to absorb the charge. This rules out the participation of some types of parton, 
making the weak interaction a more selective probe of the nucleon's structure 
than the photon probe. 
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Figure 3.1: Comparisons of the generalized structure function F2 with the elec- 
tromagnetic structure function F2 and of the structure function xF^ with F2 for 
e^p and e~p collisions, and of the full cross sections with that from pure photon 
exchange. 
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3.2 Quantum Chromodynamics and parton den- 
sity evolution equations 



The QPM has been successful in describing the early DIS data, a few serious 
problems need, however, to be solved. The first is the absence of free quark, 
which may be ascribed to some confinement mechanism which requires very strong 
binding forces between the quarks, to prevent them getting out. The second has to 
do with the seemingly violated Pauli principle for having two identical u valence 
quarks existing in a same proton. Finally, the successful description of lepton 
scattering by nucleons in terms of elastic scattering by quasi-free quarks equally 
requires very weak binding forces between them, apparently inconsistent with 
the first problem. These problems are naturally understood in the framework of 
the non-Abelian gauge theory of the strong interactions of colored quarks and 
gluons, the Quantum Chromodynamics (QCD), one of the components of the 
SU{3) X SU{2) X f/(l) Standard Model, with the SU{2) x U{1) for the electroweak 
sector of the gauge theory. 

The gauge bosons, gluons, are massless like the photons in the electromagnetic 
interaction but in contrast to the photons, they interact among themselves. This 
results in a strong scale dependence of the coupling strength as, one fundamental 
constant of QCD that must be determined from experiment (Sec.O). At small 



distance (large energy scale) the coupling between the quarks and gluons is small, 
this is known as "asymptotic freedom" . At large distance the coupling increases 
which leads to confinement. Taking as the energy scale, the coupling strength 
in leading order (LO) and next-to-leading order (NLO) is given respectively by 

[Q 



a 



NLO 



1 



1 - 



pi lnln(gVA 

Wo in(gvA2 



2M 



(3.30) 

(3.31) 



/3oln(gVA^) 

where Uf is the number of quarks with mass less than the energy scale Q^, A 
represents the scale at which the coupling would diverge [0, and f3i are functions 
which control the renormalization scale dependence of the coupling: 
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(3.32) 

(3.33) 
(3.34) 
(3.35) 



where contrary to Po and Pi, P2 is scheme dependent (see below). 
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3.2.1 Structure functions in QCD 



In QCD, the structure function definitions of the QPM are modified to accommo- 
date strong interactions between the partons and to include mass effects. Using 
the factorization theorem, the structure function can be expressed as 

fx ij? \ 

where, as illustrated in Fig.p.2|(a), i is the parton label and f ® g denotes a 
convolution integral 

f®9= f-f{y)9{-) , (3.37) 

Jx y \yj 

of the hard vector-boson-parton cross section a and the parton distribution func- 
tion xq. The hard cross section a can be calculated in perturbative QCD: 

a f-, ^, 4' o^s{^^')] = cJo6 (l - -) +a.(/i^)ai f-, ^, ^) +0{a',) , (3.38) 

where (Tq is the contribution in the QPM (Fig.^]2|(b)). In performing calculations 




Figure 3.2: (a) Schematic diagram showing QCD corrections and factorization 
theorem in DIS, (b) diagram in the QPM. 

beyond leading order, various divergences arise and the renormalization scale /i 
is introduced to regulate these divergences. The factorization scale fip serves to 
define the separation of short-distance from long-distance scale. Roughly speak- 
ing, any propagator that is off-shell h j fip or more will contribute to a while it is 
grouped into xq below this scale. Since the physical structure function F2{x, Q"^) 
is independent^ of the scheme and scale dependence, these dependences of a must 

^This is true only when the result is summed to all orders. Truncating the perturbative 
series at a given order spoils the perfect compensation and introduces an artificial dependence 
on the choice of the scale. 
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be compensated by corresponding dependences of the parton distribution func- 
tions xq. In the inclusive DIS analysis, these scales are usually chosen to be Q^, 
the characteristic large momentum scale of the process, to avoid logarithms of 
large ratios ln(Q^//i^) and ln(/i|,//i^). 

The two most commonly used schemes are the modified minimal subtraction 



MS and DIS. The MS scheme is appealing for its theoretical elegance and 
calculational simplicity. The DIS scheme |45|, on the other hand, is appealing 
for its close correspondence to experiment. In this scheme, one demands that, 
order-by-order in perturbation theory, all corrections to the structure functions 
F2 be absorbed into the distributions of the quarks and antiquarks, such that 
the higher order formula for F2 is the same as the leading-order formula. On 
the other hand, both Fl and XF3 acquire nontrivial order Us corrections. For 
example the structure function Fl{x,Q^) reads 



Fl{x,Q'' 



TT 



^ dz /a;\2 
z \z) 



(3.39) 



where c = Nf for neutrino scattering and c = ^ • for charged lepton scat- 
tering. Higher order corrections and possibly non-perturbative contributions are 



sizable [H6 



At small X (x < 10 ) the second terms dominates. In fact for 
z ~ 2.5x, a measurement of Fl is almost a direct measure of the gluon distribu- 
tion E7 . 



3.2.2 Evolution of structure functions 

The most powerful quantitative prediction of perturbative QCD is the breaking of 
Bjorken scaling in DIS. Although the parton distribution functions in the hadron 
cannot be calculated from first principle, their Q"^ dependence can be calculated 
within perturbative QCD. The scale dependence of the parton distribution func- 
tions in QCD has its origin in the interactions of the quarks and gluons via such 
elementary processes (Fig. ^.3| ) as gluon emission from quarks, q ^ qg, the cre- 
ation of quark-antiquark pairs by gluons, g — > qq, and gluon emission by gluons, 
g — *■ gg. In describing the way in which scaling is broken in QCD, it is convenient 
to define nonsinglet and singlet quark distributions: 

q""' = q^- qj , (3.40) 
q' = J2^q^ + q^)■ (3.4i) 

i 

It is understood that the parton distribution functions q^^ and q^ are functions 
of X and Q^. The nonsinglet structure functions have nonzero values of flavor 
quantum numbers such as isospin or baryon number. The variation with of 
these and the gluon distribution function g{x,Q'^) is described by the so-called 
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DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equations valid to all 
orders in a,: 



dq 
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(3.42) 
(3.43) 



The splitting functions Pij{x), representing the probability of a parton j emit- 
ting a parton i with a fraction x of the parent's momentum, is calculable in 
perturbative QCD: 



27r 



2tt 



(3.44) 



The splitting functions in leading order |4^ correspond to contributions shown 
in Fig. |3.3| . The truncation after the first two terms in Eq.( |3.44| ) defines the NLO 





g 



'88 




Figure 3.3: Leading-order diagrams contributing to the splitting function P/-. 



evolution. The splitting functions at NLO have been calculated since long, 
splitting functions at next-to-NLO are however only known partially |^ . 



The 



As we will see in later chapters, the conventional DGLAP evolution equa- 
tions have been working very successful in describing the HERA structure func- 
tion data. However, due to the approximation in which only terms involving 
(a^lnQ^)" are summed to all orders in n in LO and (a^lnQ^)" in NLO, it 
is believed that it may not be applicable at very low x when asln{l/x) ~ 1. 



Indeed, in BFKL (Balitsky-Fadin-Kuraev-Lipatov) equation |51] terms involving 
(«<, ln(l/ x))" are summed instead. The solution of the equation gives a functional 
form for the gluon distribution as 



xg X 



(3.45) 
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with uj = 121n2as/7r ~ 0.5 for ag = 0.19. Recently, the next-to-leading loga- 
rithmic (NLL) corrections have been obtained |5^ giving uj = 2.Q5as — 16. ~ 
—0.12 for the same value. Unfortunately, neither LL nor NLL value is com- 



patible with the data. For the latest development, see [^3], ^ and references 
therein. 

When is large and x is small, namely In ln(l/a;) ~ 1 but In and 
aslia{l/x) are both small, terms involving (a^ InQ^ ln(l/a;))"' have to be summed 
to all orders in n to give the so-called double leading logarithmic (DLLA) approx- 
imation. Theoretical research has been trying to find a more general approxima- 
tion scheme resulting in an equation that is applicable in these different regions. 



For example, the CCFM equation ||55[ gives BFKL at small-x and DGLAP at 
large x. 

Understanding the dynamics of the small x region is one of the fundamental 
problems of QCD. For a relatively low transverse scale of px in a process at a 
future high energy hadron machine with a center-of-mass energy ^/s, the value of 
^ = VtI \fs can be low. When x is small enough, the density of partons becomes 
very large so that partons start to interact and overlap, the perturbative QCD 
will eventually fail not because the strong coupling as is large but because the 
parton density is high. 



3.2.3 Higher twist 

The predictions of QCD discussed so far is at leading twist (twist two). In QCD, 
the structure functions have higher-twist power corrections: 

F2(x, Q2) = Ff\x, Q^){\ + /i(x)/Q2 + ■ ■ ■) (3.46) 

where h[x) has a form (1 — a;)~^ according to a phenomenological analysis of the 
structure function data of BCDMS and SLAC |^ . The higher twist contribution 



is therefore expected to be important at high x and low . 

The higher twist contribution involves presumably reinteraction of the struck 
quark with the proton remnant and thus a full calculation may have to await a 
solution to the problem of confinement. 
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3.3 Parameterizations of parton distribution func- 
tions 



One of the main strands of interest in experiments on DIS is to determine the 
parton distribution functions. This is related to the fact that the parton distri- 
bution functions have a key feature, universahty, i.e. they are independent of the 
structure functions in which they appear, and of the physical processes to which 
they are applied. Therefore, the parton distribution functions can be extracted 
from a quantitative comparison of experimental data from a wide range of phys- 
ical processes with QCD master equations, Eqs. (|3.36|) , (|3.42| ) , ( p.43|) . These can 
then be used in other applications to make predictions as well as to provide strin- 
gent tests of the self-consistency of the perturbative QCD framework itself or the 
Standard Model in general. Since any compelling indications of inconsistency of 
the Standard Model are signs of new physics, and since even direct search for new 
physics must rely heavily on understanding of the background from conventional 
physics, the systematic analysis of parton distributions is intimately tied to all 
these ventures. 



3.3.1 Global analysis of parton distribution functions 

The analysis of structure function data in extracting parton density functions has 
a long history, the earlier parameterizations dated around the late 1970's. These 
together with many new generation parameterizations are available in the pack- 
age PDFLIB |5^. There are many subtle differences among analyses by different 
groups or even a same group at different times, but the technique used can be 
broadly summarized as follows. 

• Initial scale Q^: The initial scale is arbitrary, but should be large enough 
to ensure that asiQl) small enough for perturbative calculations to be 
applicable. 

• Functional forms: Functional forms for parton distributions (valence, sea 
and gluon, or non-singlet, singlet and gluon) are assumed to be valid at Q^: 

xu„{x) = AuX^^{l - xf^Puix) , (3.47) 

xd^{x) = AdX^^il - xf^Pd{x) , (3.48) 

xS{x) = Asx^^il - xf^Psix) , (3.49) 

xg{x) = AgX^'il - xf'Pgix) , (3.50) 

where Pi{x) takes the function form 1 + ei^/x + 7jX for the MRS (Martin, 
Roberts, Stirling) group [^] and 1 + 7,jX^' for CTEQ (the Coordinated 



Theoretical- Experimental project on QCD) group In some (earlier) 
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analyses, the gluon distribution takes a simplified form 1 + ^gX. While the 
difference in a given parton density using these two functional forms does 
usually not exceed (9(1)% level over a large x range when same precise data 
are used, large differences are possible at the kinematic boundaries [30]. 



• Evolution: With a particular value of A,| the DGLAP equations ( p. 42 ) and 
( p.43| ) are then used to evolve the parton distributions up to a different 
value, where they are convoluted with coefficient functions, appropriate to 
the chosen renormalization scheme, in order to make predictions for various 
processes corresponding to the chosen measurements. 

• Data sets and comparison between data and parameterizations 
through minimization: Data sets from various measurements are selected 
and compared with the corresponding predictions. The free parameters are 
obtained by a minimizing procedure (e.g. with the minuit package |6^) and 
by taking into account the experimental errors. Some of the parameters are 
constrained by the flavor counting rules and the momentum sum rule. 

A typical set of parton distributions at = 20 GeV^ obtained by the MRS 
group is shown in Fig.|3.4]. The core constraint for these parton density 



functions comes from the DIS structure function data at low x by the HERA 
experiments, as will see in Chapter ^, and at high x from fixed-target data. As 
an example, four combinations u + u, d + d, u + d, and s can be derived from the 
following structure functions^ in the leading-order form: 

F^P - F^"" = ^x{u + u-d-d) , (3.51) 
^(^pip^pin^ = + ^ + + + , (3.52) 

F^^ = F^^ = x{u + u + d + d + 2s) , (3.53) 
^x [Ff + F^^) = x{u-u + d-d} , (3.54) 

where p, n, N stands respectively for proton, neutron, and isoscalar targets. 

The fact that u valence shape is different from that of d valence has been 
known since the earliest days of neutrino scattering when neutrino and antineu- 



trino scattering data on protons and deuterium were compared [^, The 
recent data which fix these valence shapes have come from taking the difference 
and the ratios of and F^"" from NMC At large x, when only valence 



■^Depending on the analyses, the variable A was either treated as a free parameter or fixed 
using the precise value from other independent measurements [pl[ 



*When deriving Eqs.(3.51)-(3.54). we have neglected the charm contribution and assumed 
s = s. 
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MRST partons = 20 GeV^ 
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distributions are significant, one has in leading order 



Fr ~ 4 + • ^ ^ ^ 

The charge asymmetry at the Tevatron pp coUider provides additional in- 
formation in the region of a; ~ 0.1 and ~ M^. Because the u quarks in the 
proton carry more momentum on average than the d quarks, the bosons tend 
to follow the direction of the incoming proton and the W~ bosons that of the 
antiproton. When x — ^ 1, the behavior of the ratio d^/u^ is, however, largely 
unsettled with model predictions varying between and 0.2 |]66| (Sec. p.3| ). 



In the global analysis of the parton distribution functions, the sea quark is 
assumed to be the same as the anti-quark sea, i.e. qi = q^. In the earlier analyses, 
the flavor symmetry u = d is also assumed. In 1992 the NMC data [F^ together 



with the Gottfried sum rule ^8 



/ — (F| - F^) = - dx{u, -d,) + - / dx{u - d) (3.56) 
Jo ^ 3 Jo 3 Jo 



gave the first evidence that d > u. The recent global analyses |^8|, ^ also use the 



asymmetry of Drell-Yan production in pp and pn collisions first from NA51 
for X = 0.18 and then from E886 for an extended x range (0.04 < x < 0.3) 
to determine directly u — d. The semi-inclusive DIS data from HERMES are 
expected to give confirmation or independent information on the u and d flavor 
asymmetry. 

In Fig.^.4|, the strange quark sea s is shown to be different from the quark 



seas u and d. In fact in the recent global analysis of e.g. MRST, s is assumed to 
have the same x dependence asu + d but suppressed by 50%. This is supported 



by the CCFR dimuon data [^. The suppression is presumably due to the mass 
difference between s and u and d though they all have been treated as massless 
quarks in the global analysis. 

The charm sea c is further suppressed with respect to other light quark seas. 
Until recently different groups have used different procedures. One such proce- 
dure is to treat the charm quark (similarly for the bottom quark) as infinitely 
massive below a threshold = m^, and as massless above the threshold thus 
evolving according to the normal massless evolution equations. Up to NLO in as 
this prescription guarantees that the correct results are obtained asymptotically, 
but it is unsatisfactory near the threshold. An alternative procedure considers 
the charm as being produced from the hard scattering between the electroweak 
boson and a gluon, i.e. the boson gluon fusion process (Fig.p.5|). The latter treat- 
ment incorporates the correct threshold behavior automatically but is unsuitable 
for ^ due to the unsummed potentially large logarithm in Q^/m"^. The 



recent measurements of charm production at HERA |74| have emphasized the 
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Figure 3.5: Diagram of the boson-gluon fusion in deep-inelastic lepton-proton 
scattering. 



importance of having a consistent theoretical framework for heavy flavor produc- 
tion in DIS. This has been achieved and applied in the recent global analyses of 



parton distribution functions |58, 3£]. As a consequence, contrary to light quarks, 
the charm quark density is determined by the other parton distributions and no 
extra parameters are introduced apart from the charm quark mass. 

At values of far above that shown in Fig. 3^, the bottom quark sea be- 



comes increasingly important and eventually all sea quark distributions evolve 
to a common form, concentrated at small values of x, since they are driven by 
g ^ qq transitions. 



3.3.2 Dynamical parton distributions 

Whereas the parameterizations of the global analyses depend crucially on the non- 
perturbative input parameterizations at Ql, the ones of the GRV group [^, |76|, ^ 
are constructed to be less dependent on their inputs. 

The original idea behind these parameterizations is that at some very low 
scale = /i^ with fi ~ 3A, the nucleon is assumed to consist only of constituent 
valence quarks. As increases, one generates the gluon and sea quarks in the 
nucleon dynamically from the valence quarks, through the conventional DGLAP 
equations for the processes q ^ qg, g ^ qq. The resulting predictions turned 
out, however, to be too steep in the small x region and subsequently the parame- 
terizations for the sea quarks and gluons were modified to be valence like instead 
of being null. This modification was made based on the argument that partonic 
quark distributions should rather be identified with the current quark content of 
hadrons instead of the constituent quarks. 

In an earlier version in 1992 [^], only fixed-target DIS data at high x {x > 
0.01) were used in fixing the valence-like input parameters at /i^. The resulting 
prediction at low x was in qualitative agreement with first structure function mea- 
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surements at HERA (see Sec. [4.2| ). Quantitatively, however, it could still deviate 
systematically from the later more precise data from HERA. The parameteriza- 



tions were thus updated first in 1994 [^] and then in 1998 IT^] by including the 
HERA structure function data together with other DIS and non-DIS data. There- 
fore, as far as the used data sets are concerned, the GRV parameterizations do not 
differ from the other global analyses. The only important remaining difference 
lies in the chosen initial scale which is much smaller than in the other analyses. 
It is from this low scale that the gluon and sea quarks acquire the sufficient large 
radiative evolution length In {as{^i'^) / oisiQ'^)) to change dynamically its behavior 
from valence-like (5, > 0,i = g, S, see Eqs.( |3.41^ ) and ( ^.5U| )) to sea-like (_Bj < 0) 
found in the other global analyses choosing Ql ~ 0{4:) GeV^. 



3.3.3 Other model parameterizations 

In addition to the parameterizations described in the previous subsections, there 
are a number of other models or parameterizations . One example is the model 
based on the Regge theory ||79[. The Regge theory has been very successful in 
describing the energy dependence of the soft hadron interactions: 

a ~ s"^-^ (3.57) 

where a-p = 1.08 is the intercept of the soft pomeron (the so-called Regge trajec- 
tory corresponding to the exchange of families of particles with different spin). 
The connection between the small x and the energy is made by W^, the invariant 
mass of hadronic system defined in Eq. ( p. 41) : 

= Q\l/x-l) -Q^x. (3.58) 



In this approach, the gluon distribution is expected to behave therefore as [|80[ 

xg ~ x^-"^ . (3.59) 

When this description is applied to the HERA structure function data, it is 
found unsuccessful as soon as reaches 0(1) GeV^ (see Sec. ^I^) . The model 
has been extended recently to include a hard pomeron in addition to the soft 



one to describe the HERA data ^T[. One interesting finding of the new model 
is that while the hard pomeron is needed to describe the strong rise of the F2 
as X decreases (Sec. [4.2|) , the contribution of the soft pomeron at = 5 GeV^ 



dominates at high x and can be still important for x down to ~ 0.0002, a behavior 
of the higher- twist contribution (Sec. |3.2.;^ ). 
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3.4 Radiative corrections and hard radiative pro- 
cesses 



We have given explicitly in Sec. |3.1.1| the Born cross sections corresponding to the 
contributions from the lowest order lepton-nucleon scattering processes shown in 
Fig. p.2| (b). Experimentally, the inclusive cross sections that are measured include 
higher order electroweak corrections. These corrections, in particular the con- 
ventional electromagnetic (QED) bremsstrahlung corrections, can be very large 
(> 100%) depending on the kinematic phase space and the kinematic reconstruc- 



tion method used ||8^, These QED contributions do not contain additional 

information about the strong and weak interaction part of the complete theory, 
and are generally removed from the measured cross sections. 



3.4.1 Radiative corrections for neutral and charged cur- 
rent DIS processes at HERA 

For the neutral current (NC) process, higher order electroweak contributions can 
be separated into the QED and weak corrections as shown respectively in Figs.pTB 
and |3 



The QED contributions can be further subdivided into 



• the leptonic corrections described by diagrams containing an additional 
photon attached to the electron line, i.e. the photon emission from the elec- 
tron line, Figs. |3.6| (a) and |3.6| (b), and the photonic lepton vertex correction 
combined with the self energies of the external fermion lines, Fig. |3.6| (e); 

• the quarkonic corrections represented by diagrams with an additional pho- 
ton at the quark line, Figs.|3^(c), p^(d) and |3^(f); 

• the interference of bremsstrahlung from the electron and the quark line, 
Figs.p^(a)-(d), and the box diagrams, Figs.|3.6|(g) and p^(h), which connect 
the electron and the quark line by an extra virtual photon. 



The leptonic corrections constitute the bulk of all radiative corrections and 
represent the practically most important contribution. There are two reasons 
why these corrections can be large. First, large logarithmic term of the form 
a/7rln(Q^/mg) are present due to the radiation of photon collinear with the 
emitted lepton. Secondly, the emission of a very energetic photon shifts the mo- 
mentum in the propagator of the exchanged photon to a value which is essentially 
smaller than determined from the energy and momentum of the final electron. 
The second effect is primarily of kinematical nature; it can be reduced in mag- 
nitude by applying suitable cuts (Sec. |4.3.lD and depends on the method used to 
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(e) (f) (g) (h) 

Figure 3.6: Diagrams of QED corrections to the NC DIS process. 




(d) (e) 
Figure 3.7: Diagrams of weak corrections to tlie NC DIS process. 
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reconstruct the event kinematics (Sec. p.6.I| ). The process in which the emission 
of an energetic photon occurs in the incident electron hne (Fig. p.6|) can be fur- 
ther used to measure the structure function F2 in an extended kinematic region 
(Sec.O) and the longitudinal structure function (Sec. 4.7.^) . 

The quarkonic corrections exhibit a different behavior: they are essentially 
smaller in magnitude, typically a few percent, and are flat functions of Q^. The 
main reason for this difference is the absence of the kinematical effect in the pho- 
ton propagator, which dominates the radiation from the lepton. The occurrence 
of quark mass singularities of the type a / T[e'^^hi{Q'^ / m?^ can be absorbed into the 
quark distribution functions. 

The interference corrections are not affected by either lepton or quark masses 
and therefore free of the mass singularities. 

The weak corrections (Fig. |3.7| ) are infrared finite and numerically small (of 
the order of 1%). 

Diagrams of electroweak corrections for the charged current (CC) process 
is shown in Fig.p.8|. In contrast to the NC process the subsets of the real 
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W(Z) 



(h) 



Figure 3.8: Diagrams of electroweak corrections to the CC DIS process. In addi- 
tion to (h), the corresponding crossed diagram also contributes. 



bremsstrahlung diagrams for the CC process are not gauge invariant. More- 
over the appearance of the non-abelian •yWW vertex in the real (Fig.^]^(d)) and 
virtual (Figs. [3^ (f) and p.8| (g)) corrections, not present within the conventional 
QED but typical for a non-abelian gauge theory, indicates that the classification 
convenient for NC processes is less sensible in the CC case. The radiative correc- 
tions for CC are flat functions of Q^, in contrast to the QED corrections in the 
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NC case. This is a consequence of the absence of the photon exchange diagram; 
the kinematical effect of the lowering the exchanged is less important in the 
W propagator due to the presence of the large W mass. 



3.4.2 Hard radiative processes 

In the process, ep — > e7X, where a real photon is emitted (Fig.|3l6|(a) and |3^(b)), 
the virtuality of the intermediate electron Q'^ can be defined 

_ -{k - Kf (Fig.3.6(a)) 

^ " ^ " I -{k' + Kf (Fig.|3(b)), ^"^-^"^ 

in analogy to that of the exchanged virtual boson, Q^, defined in Eq. (|2.1| ), where 
/c, k' and K are respectively the four momentum of the incident electron, the 
scattered electron, and the radiative photon. 

The presence of the real photon introduces in the cross section formula an 
additional propagator due to the intermediate virtual electron. One has therefore 

' ' (3.61) 



Depending on the relative values of and Q'^, one refers to the following pro- 
cesses: 

• Bethe-Heitler and radiative photoproduction processes correspond 
respectively to the elastic and inelastic channel when — ^ and Q'^ — 0. 
The Bethe-Heitler process has the largest cross section and has been used by 
both HI and ZEUS collaborations to provide the luminosity measurement 
(Sec.m)- 



QED Compton process corresponds to the case when — > but 
Q'^ > 0. Since the radiative photon and the scattered electron can both be 
measured in the main detector, these events have been used to check and 



calibrate the energy scale of the electromagnetic calorimeters (Secs. |4A4 
and pXBI) . 



Radiative DIS process corresponds to the case when > and Q'^ 0. 
These events are of particular interest, see Secs. |4.4| and [4.7.2 . 
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Chapter 4 



Measurement of structure 
functions and their interpretation 

4.1 Pre-HERA results and expectations 

Before the advent of the HERA experiments, various structure functions have 
been precisely measured by several fixed-target experiments for x < 0.01 and 
values of Q ranging from C(10) GeV^ to about 100 GeV^ (FigQ. As an exam- 
ple, the proton structure function F2 measured by NMC and BCDMS 



is presented in Fig.|4.1|. These data have been used by various groups to extract 
parton density functions. A few parameterizations of the proton structure func- 
tions are also shown in Fig.[4.1|. These parameterizations, which all described the 



then existing low energy data, differ at x ^ 10 ^ by more than a factor of four. 
The large uncertainty at small x arises for two reasons 

1. theoretically, there were concerns that within the perturbative QCD frame- 
work, the occurrence of powers of ln(l/x) can spoil the conventional (twist- 
2) formalism, 

2. phenomenologically, even within the standard approach, the initial parton 
distributions, which are needed in solving the evolution equations, were 
largely unknown at x < 10~^, because the existing data did not extend into 
this region, and the only constraint, the momentum sum rule, does not fix 
the shape of the distributions. 

All phenomenological analyses of parton distributions based on the usual QCD 
formalism used certain assumed parameterizations of the initial distribution func- 
tions that implicitly determined the extrapolated small-x behavior. For the 
MRS D parameterizations the small x evolution of the gluon density (at 



Ql = 4GeV^) was singular {Bg = -0.5, see Eq.(|33q)) for MRSD-' and con- 
stant (Bg = 0) for MRS DO'. Similarly, for the CTEQ IMS |86| parameterization 
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the gluon density was singular, but the sea quark density was not strongly cou- 
pled to the gluon density, leading to a slower rise of F2 with decreasing x. For 
the GRV ||75| parameterization small x partons were dynamically generated ac- 
cording to the DGLAP equations, starting from "valence like" quark and gluon 
density functions at Ql = 0.3 GeV^. 



2.5 



2.25 



Q^=15 GeV^ 



NMC 



1.75 



1.5 - 



1.25 - 



0.75 - 



0.5 



0.25 - 



BCDMS 



MRSD- 



MRS DO' 



CTEQ IMS 




Figure 4.1: The proton structure function F2 measured by the fixed-target experi- 
ments NMC and BCDMS / |2l|/ and compared with various parameterizations: 
MRS D-', MRS DO , CTEQ IMS |S|/, GRV 
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4.2 First measurements at low x from HERA 



^ of data collected in 1992, the first year of the HERA running, 

P^ , ^ at small x down to 



With about 25 nb 

both HI and ZEUS have made a first measurement 
0.5 X 10~^ at values of comparable with the low energy fixed-target data [^, 



21]. The X dependence of the measured F2 for three selected Q values is shown 



in Fig.| 



Q^=60 GeV^ 



CTEQ IMS 
GRVHO 




Figure 4.2: The first proton structure function F2 measured at HERA by HI / |8^/ 



and ZEUS /||/ compared witli 4 parameterizations: MRSD-', MRS DO j8§/, 
CTEQ IMS GRV 



The data showed a significant rise in F2 towards lower values of x. This is 
a striking feature compared to all previous structure function data obtained at 
fixed target experiments. These latter data can be successfully described 
by a single nonperturbative pomeron with intercept close to 1.08 (Sec. p.3.3|) . 
This is no longer possible for the new data. These new data, though limited in 
precision, have been used by several QCD analysis groups to constrain the parton 
density functions at low x resulting in better parameterizations of MRS H |^0| and 
CTEQ 2 
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The behavior of the strong x dependence has been subsequently 
confirmed with much improved precision already with 1993 data. The same data 
have allowed a first measurement be made at values of beyond those covered 
by fixed-target experiments. This analysis is presented in the next section to 
illustrate as an example how these measurements were performed. 
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4.3 First F2{x,Q'^) measurement at values of 

beyond those covered by fixed-target exper- 
iments 

In 1993, both HI and ZEUS have collected a factor of 10 more data than in 
1992. The integrated luminosity collected by HI was 0.271 pb~^. The kinematic 
region covered by the 1993 data is shown in Fig. |4.3| . The open region on the left 
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Figure 4.3: The kinematic region covered by 1993 data. It is limited on the top 
by the angular acceptance of the hadronic final state (6h > 5°), on the left by 
the angular acceptance of the scattered electron, and from the bottom by the 
cuts on E'^ > S GeV and y < 0.8. The shaded region, analyzed in the following, 
corresponds to the acceptance of the liquid argon (LAr) calorimeter, while the 
open area corresponds to the acceptance of the backward calorimeter (BEMC). 

corresponds to the acceptance of the backward calorimeter (BEMC) in which the 
scattered electron is detected. The first measurement of F2{x, Q"^) based on 1992 
data was made in this region. The increase of the luminosity has significantly 
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extended the kinematic coverage to the higher region (shaded area in Fig. |4.3| ) 
in which the scattered electron is measured in the hquid argon calorimeter (LAr). 
This kinematic region is covered in the following analysis. 



4.3.1 Event selection and background filters 



The selection of DIS events at high was based on an identified scattered 
electron in the LAr calorimeter and additional requirements for background re- 
jection. One of the highest neutral current events measured by the HI detector 
is shown in Fig.H^. 



Run 59384 Event 25101 Class: 8 9 12 M 16 20 22 23 



Date 1/10/1993 



Q = lOBOO GeV x = . 23 y = 0.54 








Figure 4.4: A neutral current event at high measured by the HI detector in 
1993. The event kinematics is also shown. 



The electron identification algorithm developed for this analysis used the 
salient feature of the LAr calorimeter, namely the fine granularity in both lateral 
and longitudinal direction. Requiring the longitudinal profile of an isolated en- 
ergy deposit (cluster) to be consistent with an electromagnetic object, the cluster 
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having the most compact lateral size was assigned as the scattered electron re- 
sulting in a very good efficiency and low misidentification probability . As far 
as the inclusive cross section and structure function measurement is concerned, 
the misidentification is more critical since it results in a wrong kinematic recon- 
struction for a DIS event and brings in additional background contribution from 
other events. The misidentification probability was found to be 1.2 ± 0.5% in 
data which was in good agreement with the simulation 1.5 ± 0.3%. These values 
corresponded to an identification using purely calorimetric information^. A study 
based on the Monte Carlo events showed that if the track-cluster link were used, 
the misidentification probability would have been reduced by a factor of two. 

The main non-ep background in this high analysis was due to muons trav- 
eling off axis parallel to the proton beam. These muons are produced by proton 
beam halo interactions and occasionally generate an electromagnetic shower in 
the LAr calorimeter. Other important background sources included the cosmic 
ray events and beam gas/wall interactions. Most of these backgrounds were re- 
jected by requiring 

• a reconstructed event vertex around the nominal interaction pointQ: 
l^vtx — ^ol < 30 cm with = —5 cm and 

• at least one reconstructed charged track having less than 2 cm of the dis- 
tance of closest approach between the track and the z axis. 



The remaining background events were rejected by a few topological background 
filters developed for this analysis [^. For example, the halo events were elimi- 



nated by searching for energy distributions in the LAr calorimeter which are very 
localized in the r — plane while having large spread in the z direction. These 
background ffiters have been corroborated and complemented later with other 
filters 1^ and are being widely used in the current physics analyses 



The dominant ep background was from photoproduction in which the scat- 
tered electron escaped the detector along the beam pipe and an energy cluster 
from the hadronic final state faked a scattered electron. A sizable fraction of 
the background, concentrated at low energies and high was suppressed by the 
following cuts: 

Ve < 0.8 (4.1) 
E'^>8 GeV (4.2) 
Se + S/, > 30 GeV . (4.3) 



^Thc significant part of the inner eentral tracker was not operating in the whole 1993 data 
taking period due to the presence of broken wires, thus the track-cluster link could not be used 
in the electron identification. 

^ The cut corresponds to about three times of the proton-bunch length (erf — 11cm, see 



Table ^ 
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For below 284 GeV^ the minimum energy requirement is more restrictive than 
the Uf. cut, the latter becomes effective at higher (Fig. ^4.3|) . The quantities Eg 
and Y^h are defined respectively in Eqs. (|2.14|) and (|2.15|) . Note that the sum 



Se + E/i is rather insensitive to energy loss in the forward beam hole as the 
energy and the z component of the momentum essentially cancel, while it is very 
sensitive to energy loss (e.g. the scattered electron of a photoproduction event 
or the initial state radiative photon of a DIS event) in the backward beam hole. 
While the sum for a DIS event should be around 2i?eJ^ a smaller value is expected 
for the photoproduction events as the contribution from the scattered electron 
was not included in the sum. The same cut also rejects DIS events with an 
energetic photon radiated along the electron beam direction, thus reducing the 
radiative corrections to the measurement. 

The final selected sample consisted of 1038 events. The residual non-ep back- 
ground was estimated to be smaller than 1.2% with dominant contribution from 
the cosmic ray event candidates and negligible other non-ep background events 
by analyzing the pilot bunch data and a visual scan. The overall photoproduction 
contribution was less than 1.5% and the largest contribution at high y did not 
exceed 10%. 



4.3.2 Monte Carlo simulation 



For the high analysis, three Monte Carlo samples have been generated for the 



neutral current DIS interaction using the event generator DJANGO and par- 
ton density distribution parameterizations MRSD— MRS DO', and MRSH. The 
DJANGO program is based on HERACLES for the electroweak interaction and 
on LEPTON [0 to simulate the hadronic final state. Heracles includes first 
order radiative corrections, the simulation of real bremsstrahlung photons and 
the longitudinal structure functions. The acceptance corrections were performed 
using the MRS H parameterization, which is constrained to the HERA F2 results 
of 1992. To describe higher order QCD radiative processes lepto uses the color 



dipole model |^ as implemented in ARIADNE which is in good agreement 
with data on the energy flow and other characteristics of the final state as mea- 
sured by HI ||100|| and ZEUS [|101|| . The program jetset ||102|| is then used for the 
fragmentation of the resulting partons into hadrons, and for their decay. Jetset 
is based on the Lund string model of fragmentation | 103 |. 

The "soft" vector meson contribution of the photoproduction interaction was 



simulated using the rayvdm [|104|| program, and the "hard" scattering part using 



the PYTHIA 1 105 1 program. The relative contributions of both were adjusted to 
agree with the total photoproduction cross section analysis | |106 |. 



•^Thc relation Sg 

: + = E(, — Pe.z 



2E, 



E„ - Pr 



is derived from the energy and momentum conservation: 
: 2Ee by neglecting the electron and proton mass. 
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4.3.3 Kinematic reconstruction 



As mentioned in Sec.p.6.1], the kinematics for a neutral current DIS event can 



be redundantly reconstructed. Three methods have been studied and compared 
for this analysis: the electron method, the double angle method and the mixed 
method (Mixed). The first two methods have been defined already in Sec. |2.6.1| . 
In the third method, the is that of the electron method, i.e. QMixed — Ql^ 
while the y is determined from the hadronic system, i.e. yuixcd = Vh- 

For the electron method, the energy and the polar angle of the scattered 
electron were used. For values of > 120 GeV^, the scattered electron is de- 
tected in the LAr calorimeter. The polar angle 6e was defined with the cluster 
center in LAr with the vertex position. With the 1993 data, a 5mrad accuracy 
and a 7mrad resolution were obtained by using the Monte Carlo simulation and 
comparing data and the simulation with the central tracking chambers |^, |107|| . 

Using the redundancy in the reconstruction of the kinematic variables by 
requiring Ql = Qda with Ql and Qda being defined respectively in Eqs. { \i-7\ ) and 
(P.IOI) , the energy of the scattered electron can be predicted 



,DA 



E, 



(4.4) 



where and a/i, defined in Eq.( |2.13|) , are related to the angles of the scattered 
electron and of the hadronic final state. The energy scale determined by test beam 



measurements ||108|| was thus refined using this method. The resulting systematic 
uncertainties in E'^ were smaller than 3% for the barrel part of the calorimeter 



and 5% in the BBE region and the cracks 0, |107 



The mixed method needs the hadronic final state for reconstructing y. The 
hadronic final state is measured in the LAr calorimeter^. The hadronic energy 
scale was known to 6% as determined from studies of the transverse momentum 
balance of DIS events. The test-beam data of pions between 3.7 GeV and 205 GeV 



showed agreement on 3% level with the Monte Carlo description ||109 |. 

The hadronic angle Qh was reconstructed according to Eq.( p.l3D from the 



energy deposits in the calorimeter cells. The Monte Carlo simulation showed that 
the angle was well measured except for the small and large angles; at small angles 



it is sensitive to the calorimeter noise effect (see Sec j5. 1.5] for more explanation), 
while at large angles, hadrons were not well contained in the BEMC. 

The four quantities, E'^,6e,yh and Oh, are compared in Fig.^]5| between data 
and the Monte Carlo simulation. Given the precision, the agreement was reason- 
ably good. 



"^In fact, this is only true for medium and low y. For very high y hadrons are mainly 
reconstructed in the BEMC, but these events were suppressed due to the cuts on E'^ and ?/, see 
Eqs.Q and (O). 
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Since different quantities were used in tlie reconstruction of tfie kinematic vari- 
ables, tlie resulting precision could differ. This is illustrated in Fig .[4 .61. While Ql 
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Figure 4.6: A comparison of the x and resolutions of the electron method (E), 
the mixed method (Mixed) and the double angle method (DA). The dashed lines 
indicate the bin boundary. Three y values are shown with the dotted lines. The 
averaged true kinematics is around the center of each bin. The location of the 
crosses with respect to the corresponding bin center indicates the magnitude of 
the migration of a method and the size of the crosses shows the resolution of the 
reconstructed kinematics. 



resolution is the best over all the kinematic region studied, the Xe resolution de- 
teriorates as y decreases, this can be easily understood by making the derivatives 
of Ql and ye defined respectively in Eqs.( p^.7|) and 



5Ql 



5E' 



Ql 

Sye ^ l-Ve 
Ve Ve 



5K 
EL 



tan 



-1 



(4.5) 
(4.6) 



The Xh resolution is better at lower y than that of Xe but is still worse than that 
of xda- Overall the double angle method alone gives maximum coverage of the 
available kinematic range with reasonable precision. 
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4.3.4 Measurement procedure of the structure function 



In this section, the method used to measure the inclusive cross section and to 
extract the proton structure function F2 is first introduced. Various efficiencies 
and their systematic uncertainties relevant for the measurement are then very 
briefly mentioned only as their precisions have been superseded by high statistics 
data samples to be discussed in Chapter |^. 



Method 



What is measured experimentally is in fact the inclusive cross section. The mea- 
surement of the inclusive cross section and the extraction of the structure function 
^2(0:, Q^) are performed in the following steps: 

• Bin definition: In commensurating with the available statistics and the 
size of migration of the reconstructed kinematic variables due to the finite 
detector resolution, a finite bin size AxAQ^ in (x, Q"^) plane is defined. 

• Averaged cross section in a bin: From the number of observed events in 
data A'data and the number of estimated background events A^bg normalized 
to the integrated luminosity £, an integrated cross section over a defined 
bin AxAQ^ is measured 

Cr|AxAQ2 = (4.7j 

where e is an extra correction factor for those efficiencies which are not 
simulated properly in Monte Carlo (see below). The variable A includes 
the efficiency (acceptance) and unfolding corrections, which in the simplest 
bin-by-bin unfolding method is 

= , (4.8) 

-''^gen •'-'gcn'^gen I AxAQ^ 

where A^mc and A^gcn are respectively the number reconstructed and gen- 
erated events within the bin, and £gcn is the integrated luminosity of the 
generated Monte Carlo (MC). Other sophisticated methods have been pro- 
posed | |110|| to unfold the acceptance correction and the migration effects 



associated with the finite detector resolution, the (strong) variation in the 
DIS cross section, and the QED radiative corrections (Sec. |3.4.l| ). However, 
in practice, since the Monte Carlo from which the unfolding matrix is gen- 
erated describes the data reasonably well (after e.g. the extra correction 
mentioned above and the iteration on the input structure function, see be- 
low) the different methods give similar results for the unfolded data [|111|| . 



56 



To some extent, the migration effects depend on the shape of the structure 
function, it is therefore important that the structure function parameteri- 
zation used in generating the input Monte Carlo cross section a^cn is not 
substantially different from the measured one. Otherwise (this is the case 
for the earlier analyses), an iteration is often necessary by reweightingQ the 
Monte Carlo events by the measured cross section: 



new 



MC" 



C" new I Ax AQ2 
Cgen I Ax AQ2 



(4.9) 



Differential cross section: A differential cross section at a quoted kine- 
matic point (x, Q'^) within a bin is obtained by applying a bin center cor- 
rection 



dxdQ"^ 



cr|AxAQ2 



d 0"(new)gen (-^j Q ) 



C"(ncw)gcn I AxAQ^ 
1 iVdata - A^bg C 



gen 



dxdQ"^ 

d C"(new)gen (-^j Q , 



(new)MC 



dxdQ^ 



(4.10) 
. (4.11) 



Extraction of the structure function F2{x,Q'^): Once the inclusive 
cross section is measured, to extract the structure function F2{x,Q'^) from 
the cross section formula (Eq. (|3^ )): 



dxdQ"^ 



2t\o? 



- v'^Fl + Y_xF^ 



5. 



2(1 -y) + 



R 



(4.12) 
(4.13) 



one needs to know the contribution of the interference and exchange 
to F2 (Eq. (|3.22| )), and the structure functions R{x,Q'^) | and xF3(x, Q^), 
neither of which has been measured so far in the kinematic region covered by 
HERA. However, according to the QCD prediction, R is small. Due to the 
additional suppression factor y^, its contribution was found to be well below 
1% for most of the kinematic region considered and reached up to 4.5% at 



the high y region [^. The related issue will be discussed later in Sec.^]7. 
At values of < 1000 GeV^, the contribution from the exchange and 
'-fZ^ interference is expected to be less than 1% (see Fig. |3.1|) . At higher Q^, 



the structure function xF^ becomes increasingly larger and the structure 
function F2 also deviates from F2 for the one-photon exchange. The largest 



^The reweighting is in practice performed on an event-by-event basis using the true kine- 
matics X and Q^. 



®The structure function R is generalized from the R defined in Eq.(3^) when the 
interference and Z'^ exchange contribution is taken into account. 
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correction was estimated to be 4.3% |91|] for the highest = 1600 GeV^ 
studied in this analysis. The experimental evidence of the xF^ contribution 
at HERA will be discussed in Chapter |^. In addition, a radiative correction 
(5rc should be applied in order to extract the structure function from the mea- 
sured cross section. First-order QED radiative effects were included in the 
Monte Carlo simulation. The dominant contribution, partially suppressed 
by the selection cut (Eq.( [4.3|) ), is originating from the initial state radiation 
(Sec. p. 4.1] ). The corrections are highly correlated with the method used in 
reconstructing the kinematics. Because of small separation angles between 
the final state radiative photon aand the scattered electron, the photon is 
mostly not resolved from the scattered electron by the calorimeter despite 
of its fine granularity. This considerably reduces the kinematical migration 
effects caused by the final state radiation and therefore the radiative effects 
of the final state radiation are less important than those of the initial state 
radiation. 



Efficiencies and uncertainties 



In Eq.( 4.11 ), an efficiency correction e is explicitly shown to be applied. In 
principle, all efficiency corrections could be taken into account in the Monte 
Carlo simulation and if it is properly done, no additional correction is needed. 
In practice, this has to be verified and any systematic uncertainty has to be 
propagated into the measurement of the cross section or the structure function. 
In two cases, an additional correction will be needed: 



1. when there should be a discrepancy between data and the Monte Carlo, e.g. 



the correction of the vertex efficiency in Sec. 4.4 



2. when a correction is determined after the Monte Carlo has been made, e.g. 



the correction of the trigger efficiency in Sees. 4.4 and 5.1.6 



All efficiencies were determined in this analysis from the data and compared 
with the Monte Carlo simulation. Agreement between the experimental and 
the simulated values for the individual efficiencies (trigger, vertex) was found to 
be better than 2%. An overall error of 4% was assigned due to the imperfect 
description of the various efficiencies. 



4.3.5 Results 



The proton structure function F2 measured with the 1993 HI data by this analysis 
is compared in Fig.^]7|with other analyses [112| using the electron method at low x 
and the sigma method at high x. Good agreement is observed with few exceptions. 
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For this comparison, only statistical errors are shown because the measurements 
were highly statistically correlated although there were differences in the electron 
identification and selection cuts, which can result in a systematic difference in 
the F2 measurements. Additional systematic differences can also arise because 
the analyses use different quantities in the kinematic reconstruction and do not 
have the same radiative corrections. Systematic uncertainties have been studied 
107|| . Instead of giving the long list of considered systematic 



in detail 91, 112 



sources here, it is worth simply noting that most of the uncertainties were limited 
either directly or indirectly by the low statistics data sample available then. The 
cross section measurement and its precision have been much improved with the 
increased integrated luminosity, as we will see in Chapter ^ Nevertheless, these 
measurements together with others at low | |107|| , shown in Fig. [4.8| , provided 
the basis for a QCD analysis ||113|| in which a first determination of gluon density 
was performed (Sec. 4. 6.1 ). Also shown in Fig. 4^ are the results from the ZEUS 
experiment ||114||, which are in good agreement with that of HI. 
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Figure 4.7: A comparison of the measured F2{x, Q"^) from this analysis using the 
double angle (DA) method with other analyses ^112^] based on the electron (E) 
and sigma (T, ) methods. The error bars show the statistical errors only. 
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Figure 4.8: Measurement of the proton structure function F2{x, Q'^) based on the 
1993 e~p data from HI. The results of the ZEUS collaboration are also shown 
with open squares (the ZEUS F2 data were shifted to the HI x values by using 
the parameterization in l\114] )- The curves represent a phenomenological fit to 
the HI, NMC and BCDMS data (see Ref^U^). The F2 values are plotted with 
all but normalization errors in a linear scale adding a term c{x) = 0.6{ix — 0.4) 
to F2 where ix is the bin number starting at ix = 1 for x = 0.13. 
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4.4 The structure function F2 at low from 
radiative events 



4.4.1 Extended kinematic domain using radiative events 

When extracting the structure functions from a measured inclusive cross section 
of DIS processes, one of the corrections that one has to make is the radiative 
correction. The effect is dominated for the electron scattering by the emission of 
the real energetic photons from the incident electrons. At HERA, a significant 
fraction (~ 30%) of these initial state radiative (ISR) photons are detected in the 
photon detector of the luminosity system. The photon detector, situated at about 
100 m from the interaction point, has an angular acceptance of about 0.45 mrad. 
These radiative photons can thus be considered as being emitted collinearly with 
respect to the incident electrons. The energy of the electron after the radiation 

Ef = E,-E^ (4.14) 

available for the subsequent deep inelastic interaction has effectively been reduced 
and so does the center-of-mass energy 

T^eff 

s^*^ = 4(E,-E^)Ep = ^s, (4.15) 

where Ep and s are respectively the proton beam energy and the nominal center- 
of-mass energy squared. 

When replacing E^ and s respectively by E*^ and s'^^ in Eqs.( |2.7| )-( |2.9| ), the 
resulting kinematic values can be very different from those of a non-radiative event 
depending on E^. This is illustrated in Fig. [4.9| with a dotted curve originating 



from X = 0.00015 and = 6 GeV and continuing upward to the left as the 
photon energy increases. As a consequence the limited kinematic region for the 
non-radiative events due to mainly the angular acceptance {6e < 174°) of the 
scattered electron in the backward calorimeter BEMC is significantly extended 
to lower and high x. 

A measurement of the structure functions in the extended kinematic region is 
of basic interest: 



• The uncertainty in the size of the radiative correction originates from a pri- 
ori unknown shape of the structure functions in the unmeasured kinematical 
region rather than from the technical aspects of the matrix element inte- 
gration. The measurement of the cross section in the extended kinematical 
region allows thus a direct control of this uncertainty. 

• The extended kinematical region covers an intermediate region between 
the photoproduction processes (at ^ 0) and the DIS regime (Q^ > a few 
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(GeV') 

Figure 4.9: The extended kinematic domain from radiative events (the region 
enclosed by the full lines from the top by the inclusive hadronic angle 9^ > 5° and 
from the bottom by the energy threshold of the scattered electron E'^ > 8 GeV for 
= 4 GeV (Eq. ( ]4.J6| ) ) compared with the corresponding one from non-radiative 
events (the region enclosed by the dashed lines, in addition to the limits on 6h 
and E'^, the main limit is on the angular acceptance of O^. < 174°j. The dotted 
curve and large dots illustrate how the kinematics is modified for different photon 
energies E"^. 
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GeV^). The measurement in this region may shed hght on the underlying 
dynamics of the transition. 



4.4.2 Event selection and background studies 

The analysis shown here is based on the data taken in 1994 by the HI detector. 
The data correspond to an integrated luminosity of 2.7 pb~^, which is about a 
tenfold increase with respect to the earlier data collected in 1992 and 1993, upon 
which first experimental studies of the radiative process were performed [[115| . 



Event selection 



The selection of DIS events with hard photon emission collinear to the incident 
electron is based on an identification of the scattered electron in the backward 
calorimeter BEMC, a measured radiative photon in the photon detector and ad- 
ditional requirements for background rejection. The non-ep background events, 
dominated by interactions of beam protons with residual gas and beam line el- 
ements upstream of the HI detectors, are efficiently rejected at the trigger level 
using a time of flight system consisting of two scintillator planes installed behind 
the BEMC. 

The scattered electron, which is defined to be the most energetic cluster in 
the BEMC, has to satisfy the following criteria: 

• The energy of the cluster measured with the BEMC is larger than 8 GeV. 

• The lateral size of the cluster is required to be smaller than 5 cm, as expected 
for the signature of an electron. 

• The cluster in the BEMC is required to be associated with at least one 
reconstructed space point in the BPC by less than 4 cm. 

An energetic photon detected in the photon detector (PD) is required to have: 

= EpB + Eve > 4 GeV (4.16) 

where EpD and Eye are the energies deposited in the photon detector and in the 
water Cerenkov veto counter (VC). This condition suppresses beam related back- 
ground events and cosmic rays which produce energetic showers in the BEMC. 

The remaining non-ep background events are further rejected by requiring 
an event vertex, reconstructed from tracks in the central and forward tracking 
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chambers, within ±35 cm from the nominal interaction point. The vertex po- 
sition together with the impact point of the scattered electron measured in the 
BPC/BEMCQ also defines the polar angle O^. of the scattered electron. 

The main source of ep background is pile-up events due to accidental coin- 
cidence of DIS and 7p events with a Bethe-Heitler ||116|| (BH) Bremsstrahlung 



event {ep — > e7p, Sec. p. 4. 2]) in a time window of ±5 ns. These pile-up events are 



efficiently rejected with the following cuts: 



Eetag < 2 GeV 



A = f"^'' < 0.5 (4.17) 

where i?etag is the energy deposited in the electron tagger. The quantity A com- 
pares the measured photon energy in the photon detector E^ with the measured 
missing energy -Emiss based on the main detector without including the photon 
detector: 

Emiss = Eeive - Vh) (4.18) 

with Ue and i/h being defined respectively in Eqs.( p.8|) and ( |2.20|) . One expects for 
radiative DIS events A = {E^ = -Emiss) while for pile-up DIS+BH events A = 1 
{Eraiss = 0). More details concerning these cuts and the remaining background 
events are discussed in the following paragraphs. 

To ensure a high trigger efficiency, two fiducial cuts 



^BPC + ^BPC > 15 cm (4.19) 

\xbpc\ + lyBPcl > 18cm (4.20) 

are applied to remove the region around the beam pipe where the efficiency 
degrades (Fig. [4.12| ). 



Background studies 

After having applied the selection cuts described above, a sample of 8229 events 
is selected. It consists mainly of the following types of events^: 

1. Radiative DIS events alone or in random coincidence with BH events (with 
the probability to be discussed below) where the radiated photon is detected 
in the PD, 

^It is the position measured by the BPC that has been used as it has a better spatial 
resolution than the BEMC because of the coarse granularity of the latter. 

^Note that the angular acceptance for photons from radiative DIS events is about 30% while 
for BH events it is about 98% due to the different angular distributions of photon emission. 
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2. Pile-up events due to radiative DIS with BH events where the radiated 
photon from the DIS event is not detected, 

3. Pile-up events due to non-radiative DIS with BH events, 

4. Pile-up events due to 7p with BH events or inelastic BH events alone. 

Events of type 1 constitute the signal since the radiated photons are coUinear 
with respect to the direction of the incident electron. All other events are back- 
ground. The events of type 4 contribute because of the vr" 77 decay in the 
hadronic final state, and of the electromagnetic nature of the BEMC making an 
unambiguous separation between electrons and hadrons at low energies difficult. 

The probability for random coincidence depends on the minimum photon en- 
ergy (-E™") of BH events. According to an analysis of BH events ||117|| , this 
probability is 5.6% for .E™™ = 0.13 GeV, which is consistent with the value de- 
termined using the A distribution described below. 

As mentioned above the variable A can be used to reject pile-up events. The 
distribution of A is shown in Fig. 4.10| for six different event samples. The plots 



on the left show those events which have activity in the electron tagger (-Eetag > 
2 GeV, tagged sample) while the plots on the right correspond to the non-tagged 
sample. The dependence on E'^ is shown in three different energy ranges indicated 
with the axis on the right side. As expected, the peaks around one in the two 
upper plots are due to DIS+BH events. The normalization of this background is 
fixed0 by the upper left plot and confirmed with the upper right one for events at 
high A values (A > 0.8) where the expected signal (peaking around A = 0) does 
not contribute. The difference (4% ) is taken as the uncertainty in the background 
subtraction of pile-up events of types 2 and 3. 

The normalization for the events of type 4 is fixed by the lower left plot after 
the DIS+BH events have been subtracted. In this way, any possible background 
contribution from inelastic BH events (with photons being detected in the PD) is 
also taken into account. In fact, among the estimated background events of type 
4, only about one third is expected from 7p events in random coincidence with 
BH events using the probability given above. A 30% uncertainty is assigned for 
this normalization due to the missing inelastic BH process in the MC and to the 
dependence on the low energy spectrum of the structure function. 

The relative contribution of these background events after all selection cuts is 
smaller than 10% in most of the kinematic region to be measured and the largest 
contribution at high y does not exceed 30%. 



^With this method we found a probability for random coincidence of 6.4% for E™™ = 
0.1 GeV, which corresponds to 6.0% for E'™™ = 0.13 GeV, and is consistent with the value 
mentioned above. The large value of this method may be understood since second order overlaps 
are included. 



66 



S 600 



400 



200 



Eetag>2 GeV 



HI Data 



s 







100 




> 



50 



^ 
I 60 



40 



20 








^ 600 



400 



200 



^ 
+^ 

s 

I 800 



600 



400 



200 



n RadMC 




s 







600 




400 



200 



^ 1 








30 



20 



> 



12 



8 







2 

A 







2 

A 



Figure 4.10: Distribution of A (see Eq.( \i.l7\ ) for definition) for tagged events (left 
plots) and for non-tagged events (right ones) in three different electron energy 
ranges as shown with the axis on the right side. The normalization for the 
simulated DIS+BH events is fixed by the upper left plot and the normalization 
for the '-fp+BH events is fixed by the lower left plot. 
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4.4.3 Monte Carlo simulation 



For this analysis, neutral current DIS events (both radiative events and non- 
radiative events) were generated using django |9! 



The GRV parton density 
parameterization was used for the generation of events because this parameter- 
ization provided a reasonable description of previous HI and ZEUS measure- 
ments |P7| , ^ and was one of few parameterizations that provided structure 
function parameterization at both low and high needed for this analysis. 

The 7p events were simulated using the phojet generator [|118|| , which gener- 
ates the total 7p cross section by taking into account both soft and hard processes. 

The BH events were generated according to the Bethe-Heitler approxima- 
tion ||116|| in which the proton recoil energy is neglected, i.e. -Eg + -E^ ~ E^.. The 
acceptance of the electron tagger was determined directly from the HI data. 



4.4.4 Kinematic reconstruction 

A precise reconstruction of the kinematic variables relies crucially on the mea- 
surement of the angle Qe and the energy E'^ of the scattered electron, of the 
energy i?^ of the radiative photon, and of the hadronic final state. Most of the 
measurement is better studied with the high statistics sample of non-radiative 



events ||119|| . This is the case for 6*6, -Eg, and the measurement of the hadronic 
system. The key points and the methods used for the measurement of these 
quantities are briefly discussed here. The calibration of E^^ a quantity which is 
most relevant for this analysis, is then described. 

The angle of the scattered electron is defined as the straight line between the 
reconstructed vertex and the impact point in the BPC Several effects contribute 
to the precision of the 6'g measurement : 

1. the relative alignment of the BPC with respect to the central trackers, 

2. the precision of the vertex reconstruction, 

3. the amount of dead material the scattered electron encounters on its way to 
the BPC affecting the multiplicity of the BPC hits and thus the resolution 
of the impact point. 

With the 1994 data, a systematic precision of up to Imrad is achieved. 

The energy calibration and the resolution of the BEMC is originally deter- 



mined with test beams [P-2CI|| . Both the energy scale and the resolution are checked 



and improved based on the ep data using the following three methods. 

• The kinematic peak method uses events in the pronounced peak of 
quasi-elastically scattered electrons in the region close to the beam energy. 
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The absolute energy scale of the individual stacks is calibrated. For the 
energy scale of the electron cluster, the effects of the presence of dead ma- 
terial in front of the BEMC (up to 2 Xq) and of cracks between the BEMC 
stacks are taken into account by a Monte Carlo simulation. 

• The double angle method uses the predicted electron energy E^^-da as 
defined in Eq. (|4.4|) to an independent cross check of the energy scale of the 
cluster and of its uniformity over the whole BEMC and allows to perform 
additional corrections in the crack region, which is difficult to calibrate with 
sufficient precision. 

• The QED Compton method uses events from the physical process cor- 
responding to the Compton scattering of a quasi-real photon on an incident 
electron (Sec. p.4.2|) , with the dominant contribution due to the elastic chan- 
nel {ep ep~f). The energy and the angle of the scattered electron and of 
the photon are constrained by the QED theory, making this process well 
suited for calibration of the BEMC pT| . 

The three methods are complementary and allow the energy linearity be studied. 
After the recalibration, the systematic uncertainty on the energy scale is reduced 
from 1.7% obtained from 1993 data to 1%. The dead material description is 
also improved resulting in a comparable energy resolution between data and the 
Monte Carlo. 

For a precise measurement of the energy of the radiated photon it is important 
to know the energy calibration ||117|| of the photon detector with high precision. 
The off-line calibration of the photon detector (and also of the electron tagger) 
has been determined with a sample of BH events. These events were selected 
using the following requirements: 

£^etag>4GeV 

E^>4: GeV 

24 GeV < Ectag + < 31 GeV (4.21) 
l^^ctagl < 65mm and l^ctagl < 65mm 

The last condition, where x^tag and ?/ctag are the coordinates of the impact point 
with respect to the center of the electron tagger, was used to reject events in which 
a large amount of energy leaks over the transverse detector boundaries. The 
crystal calorimeters were calibrated with events having Eye ^ 0.2 GeV making 
use of the kinematic constraint -Ectag + Ej = E^, while the water Cerenkov veto 
counter was calibrated with those events having Eye ^ 0.2 GeV. In Fig. |4.11| a) 



the correlation of E^ and i?etag is shown for those Bremsstrahlung events with 
Eye ^ 0.2 GeV. The stability of the calibration was checked by dividing the 
selected event sample into 7 subsamples. For each subsample a gaussian function 
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Figure 4.11: Correlation of the photon and electron energies for Bremsstrahlung 
events in a). The energy sum is expected to be 27.5 GeV. In b) the mean values 
of the energy sum is plotted for 7 different run intervals. The calibration is stable 
within ±0.5%. 



was iteratively fitted within ±2cr to the energy sum Ectsig+E^- The resulting mean 
values of the fits are given in Fig. [4.11[ b), separately for the two event samples with 
and without energy deposit in the water Cerenkov counter. For both samples the 
calibration could be verified within 0.5%. The relative calibration of and -Eetag 
is determined with a precision of 1.3% from studies of the E^ and Eye dependence 
of the mean + -Eetag- Taking into account a maximum nonlinearity of 1.3% in 
the response of the photon arm the global energy scale is known with a precision 
of 1.5% for E^ > 4GeV. 

The precision of the kinematic variables reconstructed with the electron method 
is similar to that shown in Eqs.( [4.5|) and ( [4.6|) except that, here, one has to take 
into account of an additional term 5E^/ E^ and (1 — ye)/yeSE^/{Ee — E^) respec- 
tively for SQl/Ql and Sye/Ue arising from the energy resolution of the radiative 
photon. Therefore, as for non-radiative events, ye degrades as ye decreases. For 
this reason, the electron method is used only for ?/ > 0.15 and for lower y, the S 
method (Eqs.( p?T6D -( p.l8| )) is used instead, which needs E/i in addition to Eg and 

At lower y, the hadronic final state is measured by the LAr calorimeter. Its 
energy scale is controlled by checking the transverse momentum balance between 
PT,h and Pr.e- A systematic uncertainty of 4% is obtained. 
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4.4.5 Measurement of the structure function ^2(2:, Q^) 



The analysis follows essentially the same method as used in the previous section. 
The only difference is that here, we measure the inclusive radiative cross section 
which is part of the QED radiative corrections to the usual DIS cross section. 
In the range of the present analysis, the effect of exchange is negligible 
and the double differential cross section for single virtual photon exchange in 
DIS, integrated over the solid angle within the acceptance of the photon detector 
9^ < 9a — 0.45 mrad and from to 2^1 = 0.85 corresponding to -E"^ > 4GeV, is 
given by [|122| 



(Pa 
dxdQ"^ 



dz 



a 



dxdQ'^dz 
2(1-2/) + 



(4.22) 



l + R 



F2{x,Q')il + 6Ho) / dzP{z) (4.23) 
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with 



E' 



eff 



Ep — E^ 



E. 



P{z) 



E. 



In 



(4.24) 



ml 



2z 
l-z 



where R, defined in Eg. (|3.8| ), is related to the longitudinal structure function, 
Fl. The higher order (ho) correction to the first order radiative contribution 
is given by 5ho(see Sec. |4.4.6]) . In comparison with the usual DIS cross section, 
the radiative cross section to be measured here is thus suppressed by a factor of 
fJo'dzP{z). 



Efficiencies and their uncertainties 



The electron identification and most of the selection cuts described in Sec .14. 4. 21 are 



similar to the analysis on the sample dominated by non-radiative events. Their 
efficiencies have been studied with this high statistics sample and found to be 
well described by the Monte Carlo simulation ||123| , |1 19|| . What are particular for 
this analysis are the trigger and vertex efficiencies. The trigger efficiency differs 
from the other analysis since a different cut on the minimum electron energy is 
applied (8GeV versus llGeV). The hard radiation studied in this analysis also 
corresponds to a different event topology, on which the vertex efficiency strongly 
depends: a scattered electron at large polar angle beyond the acceptance of the 
central tracker in the backward direction is unlikely to contribute to the vertex 
reconstruction and similarly the efficiency also drops when the hadronic system 
goes very forward beyond the acceptance in the other direction]^. 

^"Thc forward tracker can be and has been used for the vertex reconstruction, its efficiency is 
however worse than that of the central tracker due to the fact that the passive material in front 
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The events studied in this analysis are triggered by requiring a local energy 
deposit, or cluster of more than 4 GeV in the BEMC. The efficiency of this trigger 
was determined as function of the electron energy and the impact position with 
an independent data sample (dominated by non-radiative events) which were 
triggered with an independent tracker-related trigger. 

During the 1994 data taking period the four inner triangle stacks in the BEMC 
did not contribute to the trigger for most of the time ("closed triangle" period). 
During the last weeks of data taking these modules were included in the trigger 
("open triangle" period). It turned out that for the open triangle period the 
trigger efficiency did not depend on the impact position but on the electron 
energy only. For electrons with an energy of 8 GeV an efficiency of ~ 80% was 
found and from 11 GeV on it was found to be ~ 100%. 

For the closed triangle period the trigger efficiency strongly depends on the 
impact position, as can be seen in Fig. [4.12| . To ensure the trigger efficiency to be 
acceptably high an additional cut, max{|a:BPc| , |i/bpc|} > 13.5cm, was applied 
for this period, as it is indicated in Fig. [4.12[ In the region max{|xBPc| , |2/bpc|} > 
17 cm the trigger efficiency depended on the electron energy only, while in the 
region 13.5cm < max{|a;BPc| ? |2/bpc|} < 17cm it was determined as function 
of the electron energy E'^ and max{|xBPc| , ||/bpc|}- More details on the trigger 
efficiencies are given in Refs. ||124| , |125|| . 



The vertex inefficiency has two contributions: either there is no reconstructed 
vertex or there is a vertex but outside the vertex cut (±35 cm). These contribu- 
tions may be determined for radiative MC events in the following way: 

A^(all cuts except \z^t^ — Zo\ < 35 cm) ^ 

^ ^ ATl rill /^n-^-o r-\-\.rr-tr~it~\-i- ^ ^ o I ^ 'v I /-twi 1 

and 



A^(all cuts except 


^vtx 


- Zo 


< 35 cm) 


A^(all cuts except 3 ^^tx and 


^vtx 




< 35 cm) 



(4.26) 



_ iV(all cuts) 

kvtx zo\<35cm ]\f^Q\\ cuts except |zvtx — Zo\ < 35 cm) 

giving an overall value of respectively 82.6 ±0.3% and 94.4 ±0.2% for the selected 
radiative sample. Unfortunately this simple method when applied to the real data 
may bias the efficiency due to the residual background events contained in the 
denominators. 



Alternative method exists ||126|j, which is schematically shown in Fig.[4.13 



The vertex zqip is defined in this method as the crossing point between the z 
axis and a straight line formed with the impact point in the BPC and two CIP 
hits in coincidence in the r — z plane. This method deffers from the standard 
method in that Zvtx is still defined by charged particles from the hadronic final 
state when the electron track is inefficient in a DIS event or when there is no 



of it is relatively more important and that the track density is higher because of the Lorentz 
boost in the proton beam direction. 
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Figure 4.12: Trigger efEciency in the BEMC inner region as functions of the 
impact position (xbpc, Ubpc) for the closed triangle period. The efficiency is pro- 
portional to the box size with the full box size corresponding to 100%. The 
electron energy was required to be far above the trigger threshold to sepa- 
rate the spatial dependence of the trigger efficiency from the energy depen- 
dence. The cuts a) and b) are always applied, see Eqs. (14. J^ ),( |4.20| ). The cut 
c) max{|a;BPc| , |?/bpc|} > 13.5 cm was applied for the closed triangle period only 
to ensure a high trigger efficiency. 
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Figure 4.13: Schematical presentation of the CIP vertex definition, using tlie 
impact point in the BPC and two CIP hits in coincidence (figure taken from 
Ref./ini/;. 



such a track in a background event, while ^cip can not be defined if there is no 
CIP hits associated the electron track. Therefore the CIP vertex is less sensitive 
to the photoproduction background. For DIS events it has been checked that a 



good correlation was obtained when both vertices were defined ||124 . 

By applying the vertex cut using zcip, namely l^cip — Zo\ < 35 cm to both the 
numerators and the denominators in Eqs.( [4.25D and ( [4.26| ), the vertex efficiencies 



are found to be consistent between the data and the MC. Globally, one has 
^3 zvtx = 82.8 ± 0.5% and e|^rtx-^o|<35cm = 94.2 ± 0.4% in data to be compared 
with €3 = 83.8±0.4% and e|^,,^_^„|<35cm = 95.2±0.2% in the MC. In principle, 
it is straightforward to compare the data and the MC more differentially, e.g. in 
terms of (x, Q^) bins. However, due to the limited CIP geometrical acceptance, 
part of the sample in which the electrons scattered at large polar angles cannot 
be checked by this method. Here an improvement in the experimental apparatus 
is highly desirable. The installation of a backward silicon tracker certainly helps 
in this direction. 

Selecting events only at the considered {x, Q^) region, a difference shows up 
between the data and the MC when the efficiencies are plotted as a function of 
the hadronic invariant mass squared in Fig. |4.14| . The difference is more 



^^We have chosen on which the vertex efficiency depends in a similar way as it depends on 
the inclusive hadronic angle Oh- In fact, the smaller the is, the more forward the hadronic 
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Figure 4.14: Vertex efEciency as a function ofW"^ for radiative DIS events within 
the considered (x, Q^) bins in data(points) and in MC(open squares) for the E(a) 
and T,(b) methods. 



pronounced using the E method than using the S method. One explanation is 
that the former is more sensitive to higher order QED radiative contributions 
which are present in the data but are missing in the MC. If we assume that 
the difference between the data and the MC with the S method is entirely due 
to the vertex efficiency, then a correction of ~ 6% is needed at low values 
corresponding to low y values (< 0.14)[5 This correction has been applied for 
both E and S method in the low y region with a systematic error of 10% for the 
whole region using the E method, and of 10% and 4% respectively for the 
low and high values using the S method. 

After having corrected for the trigger and vertex efficiency, data are compared 
in Fig. [4.15| with the MC for various relevant quantities used in the reconstruction 



of the kinematics. The MC events generated with the GRV parameterization have 
been reweighted according to Eq.( [4.9|) by a next-to-leading-order (NLO) fit|^ to 
the measured F2 for > 5 GeV^ based on the HI non-radiative DIS events taken 



also in 1994 (see Sec. |4.6.l| ). The data distributions are found to be fairly well 



described by the MC simulation. 



system goes. The tracker covering the forward region has a smaher efficiency than the central 
tracker, see footnote [lo| . 

^^Using Eq.( 3.5S| ), the inelasticity y is related to as y ~ W^/s'^*^. 

-'^^More precisely, since the fit is done for > 5 GeV^, so for lower values, it is in fact a 
backward extrapolation of the fit. 
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Figure 4.15: A comparison between data and MC for various quantities relevant 
for tlie kinematics reconstruction. 
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4.4.6 Higher order QED correction 



As shown in Eq.( |4.23| ), we measure F2{x,Q^) using events with coUinear hard 
photon radiation in first order. The measured inclusive cross section has thus to 
be corrected for higher order radiative contributions. 

The effect of high order contributions, which is dominated by multiphoton 
emission collinear to the incident electron, is estimated with the event generator 
LESKO ||127|| : 

YFS 

1 + 5,, = (4.27) 



where and cr^^"^ are the two options of the lesko program. The first option 
includes 0{a) QED radiative corrections, and the second one describes multi- 
photon leptonic radiation in a framework of the Yennie-Frautschi-Suura exclusive 
exponentiation procedure ||128|| . The correction (5ho in %) for the E method is 



shown with the upper number in each (x, Q"^) bin (separated with the dashed 



lines) in Fig. |4.16| . The errors are statistical. The correction can be as large as 
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Figure 4.16: Higher order QED correction (5ho (in %) for the E method estimated 
with LESKO l\127\j (upper number, the error is statistical) and with an independent 
method based on hector I\129{I (lower number). The numbers in brackets are 
the assigned systematics errors. 

13% in the lowest x bins considered and becomes small for the high x (or low 
y) region. The correction has been checked with an independent method [|130 
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in which the higher order effect was simulated using the hector paclcage ||129|| 
with the varying incident electron beam energy E'^ (= — E^, E^ being the 
energy spectrum of the radiated photon in first order). The results^ are shown 
with the lower numbers in each bin in Fig. 4.16| . Reasonable agreement is seen 
and occasional large differences are due either to an intrinsic cut on > 1 GeV 
in LESKO, which introduces some edge effect for the first bin at 1.5 GeV^, or 
to a statistical fluctuation in lesko, e.g. as a function of for the x bins in 
the second row from the top. Since both methods give comparable results and 
the latter one has less fluctuation, this latter correction was used with systematic 
uncertainties given with the numbers (in %) in brackets in Fig. |4.16| . The correc- 
tion for the E method is small (< 3%) and is rather uniform in the {x, Q^) bins 
considered. 



4.4.7 F2 results 

The final proton structure function F2 combines the E method at high y (> 0.15) 
with the E method at the lower y region such that the measurable kinematic 
domain is significantly extended. The results of this analysis [|131| , |124|| have 
been cross checked by an independent analysis [ |125|| . The x dependence of F2 
at four lowest values are shown in Fig. 4. 17 and compared with the measured 
F2 from the other analyses using non-radiative events with the interaction point 
shifted in the proton beam direction by ~ 70cm (shifted vertex data), and the 
results of a similar analysis from the ZEUS collaboration ||132|| . The fixed-target 
experiment data from NMC ||133|| and a few parameterizations are also shown. 



Several observations can be made: 



This measurement has allowed to extend the HERA F2 measurement for 
g2 down to 1.5 GeV^ The rise of F2 with decreasing x observed in the 
previous HERA measurements for > 5 GeV^ [0, |107| |114|| persists 
down to = 1.5 GeV^. This extension to lower values is relative and 
remains true even after the angular acceptance of the BEMC is improved 
by the new backward calorimeter SPACAL. Indeed, when the same analysis 
technique was applied on the 1996 data collected with the SPACAL, the 
range has been further extended down to 0.2 GeV^ 



13^ 



This measurement has allowed to fill the gap between the HERA measure- 
ments at low X and that of fixed-target experiments at high x. In the 
common x region, the measurement is in good agreement with other mea- 
surements with a tendency to be slightly lower than e.g. those measured 
with a data sample (58nb~^) in which the interaction point was shifted to 



^^A cut on > 225 GcV^ has been applied as higher order contributions at lower values 
are unlikely to contribute according to Fig. 4.14 when a reconstructed vertex is required. 
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Figure 4.17: A comparison of the measured F2{x,Q'^) from this analysis with 
other measurements from HI (based on non-radiative DIS events with shifted 
vertex) ^131^ , from ZEUS (based on radiative DIS events) l\13^j , from E665 |2^/ , 
and from NMC /|2^/; also shown is an extrapolation to low of the HI QCD fit 



94 l\13]\ l to the data at > 5 GeV , and two model predictions: GRV and 



DOLA / pq/ . The inner error bars show the statistical errors only, while the outer 
ones represent statistical and systematic errors added in quadrature. 
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an averaged value of +67 cm with respect to the nominal interaction point 



in the proton beam direction (labeled "shifted vertex" in Fig. [4.17| ). The 
shifted vertex data sample has a luminosity uncertainty of 3.5% which is 
independent of the uncertainty of 1.5% of this analysis. 



All measurements are compatible with the GRV ||76[ calculation except in 
the very low x region where the measurements seem to be (systematically) 
below the prediction. 



The parameterization DOLA [^|, which is motivated by Regge theory and 
relates the structure function to Reggeon exchange phenomena, are seen to 
be ruled out by the measurement. 

The measurement agrees well with an extrapolation to low of a next- 
to-leading-order QCD fit ||131|| , based on the HI 1994 data and fixed target 
data at Q2 > 5 GeV^. 
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4.5 Precision measurement at HERA 



With an increased data sample taken in 1996-1997, corresponding to an in- 
tegrated luminosity of about 20pb~^, and the new backward apparatus, the 
HI collaboration has recently published a precision measurement on the inclu- 
sive cross section and structure functions covering 1.5 < < 150 GeV^ and 



3 ■ 10~ < X < 0.2 ||135|| . The statistical accuracy of the measured cross section is 
better than 1% for a large part of the data. The systematic precision has reached 
3% apart from the boundary of the covered region. 

Contrary to the previous HERA measurements, the new measurements are 
presented in so-called reduced cross sections, which are related to the double 
differential cross sections defined in Eg. p.9| ) 

with y+ = 1 + {1 — uY being the helicity function. Thus the dominant 
dependence in the double differential cross section due to the propagator of one 
photon exchange is explicitly suppressed in the reduced cross sections so that in 
most of the kinematic region the relation = F2 holds to good approximation. 
The only exception occurs at high y, where can be substantially different from 
F2 due to the sizable contribution from the longitudinal structure function F^. 

The reduced cross sections are shown in Fig. |4.18| together with the fip data 



by the NMC and BCDMS |21| experiments, and a new HI NLO QCD fit to 
the new low data and the published e~^p high data taken in 1994-1997 [Q 
(Sec.|^), for 3.5 < < 3000 GeV^. The cross section rises at low x. This rise 
is observed to be damped at the smallest values of x which is attributed to F^, 
see Sec.[4.7|. 

The data at < 3.5 GeV^ seem to overshoot the backward extrapolation of 
the QCD fit. A fit with Q^^^ = 1.5 GeV^, however, describes the low x data well. 
Theoretically one expects higher order logarithmic and power corrections to be 
larger for — 1 GeV^ such that a NLO DGLAP treatment may be inadequate. 
Further exploration of this interesting effect requires low x, high precision data 
at low Q2 ^ 1 GeVl 

The proton structure function F2{x,Q'^) can be extracted from the reduced 
cross sections by applying a (small) correction for the Fl contribution in Eq. ( [4.28| ). 
The resulting F2 is shown as a function of in Fig. ^.l9| . Also shown are the 



structure functions obtained from HI at high Q^{> 200 GeV^) and from the 
fixed-target data by SLAG H, NMG |3|, and BGDMS |2l| at low and 



high X. To extract the F2{x, Q^) from the inclusive cross sections measured with 

^^In the considered kinematic region, the jZ*^ interference and exchange contributions 
can be safely neglected. 
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Figure 4.18: The precision measurement of the inclusive reduced cross sec- 
tions from HI / |I35| / in comparison with the measurements from NMC l\13!i] and 



BCDMS The full curves represent the corresponding theoretical expectation 
based on a NLO QCD fit to the HI low data and the published high data 
taken in 1994-1997 for 3.5 < < 3000 GeV^ The dashed curves show the 
extrapolation of the fit towards lower Q"^. 
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Figure 4.19: A compilation of the proton structure function F2{x,Q'^) from the 
1996-1997 HI e+p data at low Q^{< 150 GeV^) and low x ^TSBjl , from the 1994- 
1997 HI e+p data at high Q^{> 200 GeV^) (Sec!^, and from the fixed-target 
data by SLAC |3|/, NMC |T3f , and BCDMS The F2 values are plotted in 
linear scale adding a constant c{x) = O.Q{i — 0.4) where i is the x bin number 
starting at i = 1 from x = 0.65. The full lines correspond to the NLO QCD 
fit il35|/ . 
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the 1994-1997 HI e+p data at high H (Sec.|]2|), the increasingly important 
contributions from the '~fZ^ interference and Z° exchange have been subtracted. 

The HI NLO QCD fit, based on the HI data alone, provides a good overall 
description to all data with a few exceptions at kinematic boundaries. The strong 
scaling violation is clearly displayed. It is this data feature and the precision of 
the data which allows the strong coupling constant «s and various parton density 
distributions in particular the gluon distribution to be determined within the 
QCD framework(Secs.^]6|, and [4.8| ). 
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4.6 Current knowledge of the gluon density 



4.6.1 Impact of the HERA F2 data on the determination 
of the gluon density at low x 

About half of the proton's momentum is carried by gluons. Despite this, the 
determination of the density of gluons g{x, Q"^) in the proton has turned out to 
be a difficult task. Before the advent of the HERA experiments, our knowledge 
on the gluon density distribution was very imprecise and limited to the region 
at X > 0.01 ||136| , |137|| . The principle difficulty in measuring the gluon density 
in deep inelastic scattering is that the gluons contribute only in high order pro- 
cesses through gluon bremsstrahlung from quarks and quark pair creation from 
gluons (Fig.^. At small x{< lO'^) the latter process dominates the scaling 
violation The HERA data on the proton structure function F2{x, Q"^) at low 
values of x down to 10~^ can thus be exploited to extract the gluon density in 
the kinematic range. 

Indeed, already with the first F2 measurement in 1992, an extraction has been 
performed by HI ||138|| and ZEUS ||139|| using the following approximation relation 



under the assumption that the quark contribution at low x is negligible 140 



It was derived that the gluon density rises strongly toward low x. 

With the improved F2 measurement in 1993, both HI and ZEUS have per- 
formed for the first time NLO QCD fits to F2 in a similar manner as in the 
global analyses by the MRS and CTEQ groups (Sec.|3.3D. The emphasis is how- 
ever different. While in a global analysis, universal parton density functions are 
extracted, the analyses by the experimental collaborations emphasize the deter- 
mination of the gluon density at low x. The result of the fit to F2 is also used 
to have a better parameterization of the cross section or the structure function 
than the initial one used in the Monte Carlo generation. For this reason, only a 
minimum number of data sets have been used in the analyses; in addition to the 
HERA F2 data at low x, a few data sets at high x from fixed-target experiments 
(e.g. BCDMS and NMC) were used. This is in contrast with the global analyses 
where in addition to the structure function data, other constraints from inclusive 
jet (and/or dijet) cross sections and the prompt photon data have also been used 
as illustrated in Fig. [4.20| . 

The gluon density for = 20 GeV^ obtained by HI [TT^ with the 1993 data 



is shown in Fig. [4.21| (a) in comparison with the result of ZEUS [|141|| and other 
determinations from the global analyses. The early parameterization MRS DO', 
which did not use HERA data, is clearly disfavored by the data as we have seen in 
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^ 22.5 




Figure 4.20: Various experimental measurements used hy global QCD analyses to 
constrain the gluon density xg at different x range. The resulting xg from the MRST 
group j 58] at Q'^ = 20 GeV'^ is shown with the curve. The ranges in x for the various 
measurements are only indicative since they are dependent. For instance, the HERA 
data reach to much large x values than indicated at higher though the statistical 
precision of the data there is still limited. 



Fig. [4.2| . The other par ameterizat ions which used the HERA F2 data are in good 
agreement with the results of HI and ZEUS. The error band represents the ex- 
perimental errors with the statistical and systematic errors added in quadrature. 
The experimental collaborations have the advantage of knowing the correlations 
between their systematic errors so that they could be properly treated [|144 



As the integrated luminosity increases, the accuracy of the gluon density is fur- 
ther reduced by HI (Fig.|;5l|(b)) using their 1994 F2 data [|131|] and by ZEUS 
(Fig. |4.2l| (c)) using both 1994 and 1995 data ||143|| . The most precise determina- 



tion from HI is shown in Fig.^:^(d). It is a result of a NLO DGLAP QCD fit 
to the recent HI ep cross section data (sec. [4.5|) and the BCDMS fip data. An ex- 
perimental precision (the inner error band) of 3% at a; ~ 10~^ and = 20 GeV^ 
has been reached for the first time. The gluon distribution has been obtained 
together with the strong coupling constant (Sec. [4. 81) . The effect of the uncer- 
tainty of the latter is shown with the middle error band. The outer error bands 
represent the uncertainties related to the QCD model and data range. For the 
first time, the gluon distribution has been determined with the HI inclusive data 
only. Both determinations are in good agreement. 
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Figure 4.21: The time evolution of the determination of the gluon density xg{x) at 
HERA (a) the HI result in 1993 ^J^] compared with the ZEUS result |l4i|/ and param- 



eterizations CTEQ3M g4|/, MRS (DC/) and MRS(G), (b) the HI result in 1994 pl|/ 
and (c) the ZEUS result based on their 1994 and 1995 data |J4j/ , and (d) the most 
recent one from HI using data taken in 1996-1997. 
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4.6.2 Uncertainty and future improvements 



For a wide range of theoretical and experimental applications, it is important to 
know the range of uncertainties of the parton density functions. However the 
task of deriving a reliable uncertainty is extremely complex and difficult. The 
difficulty is related to many sources of uncertainty to be considered, to name a 
few: 



The choice of experimental measurements and data sets. As an example, 
the gluon density xg{x) at = 25 GeV^ from the most recent parameter- 
izations MRST |^ and CTEQ5 is compared in Fig. [4.22[ Both groups 
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Figure 4.22: Comparison of the gluon density functions from MRST ^/ (the 
parameterizations MRS98-1,2,3 correspond respectively to the standard, high- 
gluon and low-gluon options ) with those from CTEQ5 J6^/ (the parameterizations 
CTEQ5M and CTEQ5HJ are respectively the standard and large-gluon options, 
the latter one being tailored to accommodate the Tevatron inclusive jet data, 
the high pt tail of which was in excess with respect to the prediction based on 



previous parton density parameterizations). The figure is from Fig. 19 of ^ 6i / 



have used the HERA F2 measurements based on the data recorded in 1995 
and before, which were the only available constraint at low x (see Fig. |4.20|) . 
Consequently, there is little difference at low x. On the other hand, the 
difference is much larger at medium and large xj^ where in addition to 



^^The diff iculty stems also from the fact that the gluon density becomes very small towards 
large x (Fig. [4.2l[) 
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the DIS structure function data from fixed-target experiments, the prompt 



photon data from WA70 ||145|| and E706 ||146|| were used in the analysis of 
MRST, while the inclusive jet data from CDF [|147| and DO |48| were used 
instead in the analysis of CTEQ5. 

The experimental errors. The most precise structure function data are 
dominated by the systematic errors. The non-trivial part of these are the 
correlated systematic errors, which are not always available for the global 
analyses. 

The technical uncertainties due to the freedom in choosing the initial scale 
Qq and the functional form at Qq with varying number of free parameters, 
and the internal correlation among different parton density functions and 
between different x range. 

The theoretical uncertainties due to effects such as the higher-order correc- 
tions, scale- and scheme-dependence, soft-gluon resummation, higher-twist 
effects, and nuclear (deuteron) corrections. 

The evolution uncertainty related to the values of the strong coupling con- 
stant as and the strategies applied. The determination of xg{x,Q'^) is 
strongly coupled with a^. One strategy is to determine as-, gluon and quark 
density functions together, an alternative is to take and its uncertainty 
from other independent measurements. 



These uncertainties are often internally correlated. However one thing is clear, 
namely in order to reduce the uncertainties, both experimental and theoretical 
efforts are needed. On the experimental side, more precise data will directly 
reduce the experimental uncertainty. The precise data in an extended kinemati- 
cal range also verify whether the chosen functional form is adequate within the 
conventional theoretical framework and test the validity of the latter. On the 
theoretical side, one of the most urgent tasks is to have a better understanding 
of the comparison between the measured prompt photon data (E706) with the 
QCD prediction so that the uncertainty of xg{x, Q"^) at large x can be reduced. 

More independent (direct) measurements are highly desirable to either provide 
more constraints or check whether new physics phenomena have been artificially 
absorbed into the chosen functional form. As far as the gluon density is concerned, 
several measurements could be used either as testing ground or as a source of new 
constraints. 

One such measurement is the open charm production at HERA. The pro- 
duction of the charm quark at HERA proceeds in perturbative QCD almost 
exclusively via photon-gluon fusion, where the exchanged photon interacts with 
a gluon in the proton by forming a quark- ant iquark pair (7(7 — > cc, Fig. |3.5| ). This 
holds both for DIS and for photoproduction where the exchanged photon is quasi 
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real. Based on the data collected from 1994 to 1996, the HI measurement [ [149|| is 
shown in Fig . [4 .231 . The determination of the gluon density is in good agreement 
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Figure 4.23: The gluon density at = 25GeV^ determined from the open charm 
data collected by HI in 1994-1996 1149^] compared with the preliminary gluon density 
at = 20GeV^ based on the indirect scaling violation data 1996-1997. 

with the other determinations described above. The errors were however still too 
large to be competitive with the other determinations. Both high luminosity at 
HERA and improved experimental setups (e.g. use of silicon vertex detectors) 
will allow considerable improvements in the future. 

Inelastic J/ ip photoproduction at HERA has been suggested as a measurement 
which could allow one to measure the gluon density. However it appears that 
the perturbative calculation does not behave well in the limit PriJ/'^P) 0, 
and if the small px region is excluded from the analysis the predictions are not 
very sensitive to the small x behavior of the gluon. Elastic (diffractive) J/ip 
production in DIS and in photoproduction is more promising, since the cross 
section depends on {xg{x, Q ))^, where the scale of the process is given hj Q = 
(Q^ + Mjy^)/4 with and Mj/^ being respectively the virtuality of the photon 
and the rest mass of the J/ ||150|| . These data could give information on the gluon 



distribution in the region 10 < x < 10 . At the present time the theoretical 
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framework for extracting the gluon distribution from these measurements is still 
under development and the experimental precision of the data is still fairly low. 

From Eg. ( p. 391 ), one sees that the longitudinal structure function Fl receives 
a direct contribution from the gluon, which dominates at low x. Under certain 
conditions, a measurement of is almost a direct determination of the gluon. 
Therefore it is important to measure the longitudinal structure function at HERA. 
This will be the subject of the next section. 
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4.7 Longitudinal structure function Fi{x,Q'^) 



The longitudinal structure function is an important quantity to measure. The 
earlier results on the smallness of R supported the assignment of half-integer spin 
for the quarks. A precise measurement allows an independent check of the gluon 
density derived from the indirect scaling violation of the structure functions. 
The knowledge of Fl is also needed to extract in a model independent way the 
structure function F2 from the measured cross section. 



4.7.1 Current knowledge of Fl{x,Q 



Using Eqs.( |3.4| )-( |3.8| ) and neglecting the target mass term, the differential cross 
section for one-photon exchange given in Eq. (|3.3| ) can be rewritten in terms of R 
and (Tt as|^ 



a 



dxdQ'^ 2'KxQ'^ 
where y+ = 1 + (1 — y)^, R = crL/aR, and 

2(l-y) 



F+ [1 + eR{x, Q2)] arix, Q^) , 



(4.30) 



(4.31) 



is the polarization of the virtual photon exchanged in the process. Therefore 
for measuring i? at a given (x, Q^) point, it is necessary to vary e, i.e. y or the 
center-of-mass energy squared s because y is related to s as 



xs 



(4.32) 



The difficulty in measuring R or the longitudinal structure function Fl is then 
directly related to a good control of the relative normalization and systematic 
errors at two different energies. This explains why measurements of Fl^x^Q"^) 
are so delicate, and why only a few fixed-target results have been published. 

The most extensive results from early fixed-target experiments were obtained 



by BCDMS and SLAC who has reanalyzed the data ||152|| , covering respec- 
tively the kinematic range 0.07 < a; < 0.65, 15 < < 59 GeV^ by BCDMS 
and 0.01 < x < 0.9, 0.6 < Q"^ < GeV^ by SLAC. The data lie in a region 
where non-perturbative effects are likely to be important. Most recently R has 
been measured by a new SLAC experiment E140X [ |153| |, by CCFR | J154|| and 
NMC [E2i, the latter reached a lowest x value at 0.002, which is still one or two 



orders of magnitude larger than the relevant kinematic region at HERA. 

^^In writing Eq.( 4.30| ), we have used the Hard convention |151| K ^ v = Q'^ / {2xM) for the 
flux of the virtual photon in analogy to the real photon. 



92 



At HERA, as the integrated luminosity increases, the latest cross section data 
become precise enough to see the sensitivity to at the high y region (see 
Fig. [4.18D . Two methods have been used to extract Fl{x,Q'^) by HI. This first 
method [ 155 |, used for > 10 GeV^, assumes that the dominant contribution to 
the cross section from F2 can be described by perturbative QCD by an evolution 
from lower at low y < 0.35, where the contribution of Fl to the cross section 
is small. The second method ||135|| , applied for < 10 GeV^, is derived from the 
cross section derivative 

da, dF2 p.22-y OFl y' , . _ 

FL2y - -wT—^TT (4.33) 



d\\iy dlny Yf dlnyY^ 

taken at fixed for y = Q'^/sx. At high y the second term dominates ||135|| . 



This is in contrast to the F^ influence on the reduced cross sections ar where 
the contribution of F2 dominates for all y. A further important advantage of the 
derivative method is that it can be applied down to very low where the first 
method can no longer be applied as reliable assumptions on the QCD description 
of F2 are prohibited. The most recent extraction is shown in Fig. [4.24| in compari- 
son with the other direct measurements mentioned above. The extracted Fi is in 
good agreement with the expectation from the NLO QCD fit. It should be em- 
phasized however that such an extraction does not represent a real measurement 
but a determination under specific assumptions. 



4.7.2 Feasibility study and a new method for a direct mea- 
surement of Fl at HERA 

A direct measurement can be achieved by running the collider with reduced beam 
energy [FSS], but this procedure has the obvious draw back that a significant 



running time is lost for high energy physics and that the collider is not operated 
in optimal conditions. For the Fl measurement itself, in addition to the relative 
normalization uncertainty mentioned above, a major experimental problem is the 
photoproduction background, when a hadron is wrongly taken as the electron 
candidate. 



Another method proposed by Krasny et al. ||157|| makes use of the radiative 



events that we have discussed in Sec. |4.4| , namely those DIS events in which a 
real photon has been emitted in the direction of the incident electron beam, 
which corresponds to an effective decrease of the beam energy. The spectrum 
of measured photon energies induces, for given x and values, a continuous 
distribution of the y, s and e variables. The relations ( [4. 311 ) and ( [4.32| ) stay the 



same, but the value of s changes from event to event, depending on the photon 
energy. 

The advantages of this method are that it can be used in parallel with normal 
data taking, that it avoids luminosity normalization problems, and that the sta- 
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Figure 4.24: The extracted longitudinal structure function Fl{x,Q ) by HI ll3, 



in comparison with the direct measurements by charged lepton-nucleon hxed-target 



experiments BCDMS gj/, NMC and SLAC p5^ . The inner error bars denote the 
statistical error. The outer error bars represent the systematic and model errors. The 
curves and error bands show the NLO QCD fit to the recent HI low data Ji35| / and 
the high from 1994-1997 ^ for y < 0.35 and > 3.5 GeV^. 
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tistical and systematic precisions increase continuously during data taking. As 
the basic principle of the R{x, Q^) measurement is to perform a linear fit of cross 
section, Eq. (^:30D , function of e. The -R(x, Q'^) is thus determined from the 
slope of the cross section. This procedure is independent of the knowledge of 
F2, however, it requires a very large integrated luminosity (of the order of 200 
pb~^) to provide a significant measurement of R (with a statistical precision of 
20% -50%) p7| |. 



In Ref. ||158|| , a new method was proposed, in which we have considered to 
use the available measurements of F2 and to exploit the dependence on R of the 
shape of the e distribution itself for radiative events. This is possible because 
the structure function F2 has been measured very precisely at HERA, and in 
view of the good understanding of DIS processes with initial state radiation as 



demonstrated by the HI ||131|| (Sec. [4.4|) and ZEUS [ |132|| analyses. 



As an illustration, the e distribution is presented in Fig. [4.25| for i? = (full 
line) and R = 00 (dashed line), for the full kinematic domain covered by the 



five bins in x and Q shown in Fig. 4. 26. Three bins are designed for value of Q 



2 



around 2 GeV^ (with x ranging from 4-10 ^ to 2 ■ 10 ^), and two bins for Q'^ 
around 5 GeV^ (with x ranging from 10~^ to 3 • 10~^). 

In Fig. [4.27| , the e distributions are presented for each of these 5 bins. These 



distributions are obtained from a Monte Carlo simulation using the GRV [76] 
parameterization of the proton structure function F2. 

The simulated integrated luminosity corresponds to 10 pb"^, which is the 
luminosity expected from data at the time of the study. The electron and proton 
beam energy is respectively 27.5 GeV and 820 GeV. The following kinematical 
cuts are applied for the feasibility study: 

E'^>2 GeV , 
Be < 177° , 

E^>4 GeV , (4.34) 

where E[, and 6e are the energy and the polar angle (defined with respect to 
the proton beam direction) of the scattered electron. As will be discussed later 
(see Fig. [4.31j ), the sensitivity of the method depends significantly on the lowest 



energies E'^ which can be accepted. With the new HI backward detectors, these 
are expected to be below 5 GeV and possibly as low as 2 GeV, the value thus 
chosen here. 

The R dependence of the e distribution can be studied using the variable 
p(i?;£o); defined as the ratio of the numbers of events with e smaller or larger 
than a chosen value Eq: 

^^^'^°^ = iV(i?;.>.o) - ^'-''^ 
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Figure 4.25: Distributions of the e parameter of radiative events generated ac- 
cording to Eq. ([4.34|) using the GRV parameterization of the structure function 
F2 for R = (full line) and R = 00 (dashed line), and for the kinematic domain 
covered by the bins shown in Fig. [4.26| . The ratio of the two distributions is shown 
in the inset. The electron and proton beam energy is respectively 27.5 GeV and 
820 GeV. 



Fig. |4.28| shows the p{R; Eq) dependence on R in the selected (x,Q^) bins, for 
the typical input value R = 0.5 as obtained by HI | |155| in a first analysis using one 
of the methods mentioned in Sec. |4.7.1| Since each bin covers a different e range, 
the chosen optimal value So is bin dependent. The dashed curves show the p{R) 
distribution for an input structure function F2 modified by ±10% at x = 10^"^, 
the modification decreasing linearly to ±5% at x = 10~^. This corresponded 
to a conservative estimate of the uncertainty on F2 at the time of study. The 
uncertainty has been improved since. The grey bands correspond to the statistical 
precision of the p{R) measurement for an integrated luminosity of 10pb~^. 

The measurement of R is deduced from the intersection of the grey bands 
with the spread of curves describing the p{R) dependence of the input structure 
function F2. The inner error bars in Fig. |4.29| show the statistical precision of the 
R measurement for the cuts ( |4.34|) , an integrated luminosity of 10 pb~^ and the 
quoted uncertainty on F2. 
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Figure 4.26: Selected bins. 



To estimate the sensitivity of the proposed method to several experimental 
parameters, the variable S is defined as 



S(£o,Ai^2) 



|p(i? = 0;eo) - p{R = oo;eo)\ 



(4.36) 



a 



p,R=0 



+ a 



p,R=oo 



for a given choice of Eq, of the integrated luminosity £ and of the input struc- 
ture function F2, ap^R being the statistical error on p{R), estimated through the 
Monte Carlo simulation. This variable quantifies the possibility of distinguishing 
between the two extreme values of i?: R = and R 



00. 



Fig. |4.3CI| shows that the sensitivity S for each [x, Q"^) bin is only weakly de- 
pendent on the Eo value over a rather large domain in Eq. It is found that it also 
depends little on detector smearing effects. 



On the other hand, as can be seen in Fig. |4.31| , the sensitivity E is strongly 
dependent on the detector acceptance conditions, in particular the electron en- 
ergy threshold which is related to the y and e ranges. For the same luminosity, 
the sensitivity is enhanced by a factor of 2.1 for E'^ decreasing from 6 to 2 GeV. 
A decrease on the photon energy threshold also improves significantly the 
sensitivity. The lowering of the electron energy threshold is a challenge for the 
HERA experiments because of the significant background from photoproduction 
interactions in which low energy hadrons are misidentified as the scattered elec- 
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Figure 4.27: Distributions of the e parameter in the selected bins, for the GRV 
parameterization of the structure function F2 and for R = 0.5. 



tron. 

Studies have been performed of the effects of experimental uncertainties, which 
are in general bin to bin dependent. The detector resolution was simulated using 
realistic smearing functions for the HI experiment; in addition, systematic uncer- 
tainties were taken into account (1% on E'^, Imrad on 6^ and 1.5% on E^). The 
effects of these uncertainties, combined in quadrature with the statistical errors 
and the effects of the uncertainty on the structure function F2, are displayed as 
the outer error bars on the R measurements of Fig. [4.29| . 

The subtraction of the remaining photoproduction background is another im- 
portant source of systematic uncertainty, which affects mostly the lower x bins. 
There, it was found to induce systematic errors of the same order as the errors 
due to detector resolution. One more source of systematic error will be the over- 
lap of non radiative deep inelastic events with bremsstrahlung events for which 
the photon is detected in the photon detector and the scattered electron is not 
detected. An electron tagger with a large energy acceptance is an important tool 
to reduce this background. Finally, as far as the uncertainty on the structure 
function F2 is concerned, it is observed in Fig. 4.28| that it does not imply a large 
systematic uncertainty on R. 

Taking all these effects into account, we concluded that for an integrated lu- 
minosity of 10 pb~^ a statistical precision of > 30% can be achieved under the 
considered experimental conditions and the statistical errors dominate over the 
systematic errors in most of the chosen bins. With increased statistics, a signifi- 
cant improvement of the measurement precision is thus to be expected. Detailed 
optimization studies are also expected to improve the measurement precision. 



98 



0.048 



0.046 



0.044 



0.042 



Bins 




10 10 1 



Figure 4.28: The dependence of p{R) on R for the GRV parameterization (full 
curves) and for the modified parameterizations described in the text (dashed 
curves), in the selected bins. The grey bands correspond to ilex statistical errors 
for an integrated luminosity of 10 pb~^ and for R = 0.5. 
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Figure 4.29: Typical precision of the R measurement in the selected bins, under 
the experimental conditions specified in Eq.( |4.34D . The inner error bars show the 
measurement precision for an integrated luminosity of 10 pb~^, taking also into 
account the uncertainty on the structure function parameterization described in 
the text. The outer error bars include, added in quadrature, the effects of the 
uncertainties on E'^, 6e and E^. 
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Figure 4.30: Dependence of the sensitivity S on the e parameter in the selected 
bins, for the GRV structure function. 
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Figure 4.31: Sensitivity S as a function of a) the integrated luminosity; b) the 
electron energy threshold E'^] c) the electron angular acceptance 6e] d) the pho- 
ton energy threshold E^, for the GRV parameterization under the experimental 
conditions specified in Eg. ( ^.341) and in the total kinematical region covered by 
the selected bins. 
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By now, we have collected by each HERA experiment a factor of ten more 
luminosity than what has been used in the feasibility study, therefore statistically 
there should be more than enough luminosity to provide a first direct measure- 
ment. The challenges are to achieve the best sensitivity with the lowest possible 
electron and photon energy thresholds and to keep the experimental systematic 
uncertainties under control. Currently the method is being used in an attempt 
to realize such a measurement. 
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4.8 The strong coupling constant as 



The strong coupling constant a^, the only fundamental parameter of the strong 
interaction sector of the Standard Model, has been measured in last decades by 



a great variety of processes ||159|, |160|, |161||. Here we briefly mention a few recent 



measurements of as{Mz) based on DIS data, in particular those using the HERA 
data showing the impact of the precision measurement of the structure function 



data [|T3| (Secpj). 

One of the methods used on fixed-target data is the analysis of sum rules. 
One such measurement is derived by using the nonsinglet structure function data 



X-F3 in the Gross-Llewellyn Smith sum rule ||162|| which is known to order a 



[ dx [F3^^(x, Q2) + F^''{x, g2)] = 3 [1 - a, (1 + 3.58a, + 19.0a^ - ^HT)] 
Jo 

(4.37) 

where a, = as/vr, and the higher-twist contribution AHT is estimated to be 
(0.09 ± 0.045)/Q2 in Refs.|T6|, |163 and to be somewhat smaller by Ref. [|l65|] . 



The CCFR experiment [|166|| , combines their data with that from other exper- 
iments [|163 and gives a,(\/3GeV) = 0.28 ± 0.035(expt) ± 0.05(syst)+°;°;^^(th). 



The systematic uncertainty is dominated by the extrapolation of the integral to 
the regions x < 0.01 where no measurements exist, and to x > 0.5 which is sub- 
stituted by F2 data since it is more precise and has little correlation to the poorly 
known gluon density at high x. The error from higher-twist terms dominates the 



theoretical uncertainty. If the higher-twist result of Ref. ||165|] is used, the central 
value increases to 0.31 and corresponds to pl[ 



as{Mz) = 0.118 ± 0.011 . (CCFR, XF3 sum rule) (4.38) 

The original and still one of the most powerful quantitative tests of pertur- 
bative QCD is the breaking of Bjorken scaling in DIS. The earliest and many 
subsequent determinations of a, in DIS were obtained by analyzing the scaling 
violation of the structure function data. As an example, a combined analy- 
sis of SLAC and BCDMS F2 data in a range from 0.5 to 260 GeV^ gives 
as = 0.113 ± 0.005 ||56[. For several years, this result together with the one. 



a,{Mz) = 0.111 ± 0.002(stat) ±0.003(syst) ±0.004(th), obtained by CCFR [|T68 
based on xFaQ (and also F2 at high x for the same reason mentioned above thus 
resulting in a statistically more precise measurement of a,) is the most significant 
from DIS data. However, these measurements are numerically smaller than typ- 
ical values obtained from e+e~ annihilation, as{Mz) ~ 0.120, thus raising spec- 
ulations about possible explanations. These speculations came to a halt when 



^*The advantage of using the less precise structure function xF^ instead of F2 is that it has 
no correlation to the gluon density functions. 
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the CCFR collaboration corrected their previous result, due mainly to effects of 
a new calibration of the detector, to [pO| 



as{Mz) = 0.119 ± 0.002(expt) ± 0.004(th) . (CCFR, xF^) (4.39) 

Recent scaling violation analyses include also HERA structure function data 
extending thus significantly the lever arm in the evolution. In Ref. [[169|| , a 
NLO QCD fit is performed using the early F2 data of HI from 1993 [|10^ (Sec.^) 
giving 

as{Mz) = 0.122 ± 0.004(expt) ± 0.009(th) . (BP, HI F2 93) (4.40) 
The dominant part of the theoretical error is from the scale dependence; errors 



from terms that are suppressed by 1/ ln(l/x) in the quark sector are included ||170 
while those from the gluon sector are not. 

In another analysis ||171|| , where the 1994 F2 data by HI ||131|| (Sec 



and ZEUS ||172|| are used together with the fixed-target data from SLAC |3| 



BCDMS |2T[, and E665 ||173|| in a fit including all known NNLO terms, the re- 



sulting as is 

as{Mz) = 0.1172 ± 0.0017(expt) ± 0.0017(th) , (SY, HERA F2 94) (4.41) 

where the theoretical error includes the uncertainties in the quark masses, higher- 
twist and target-mass corrections, and errors from the gluon distributions. But 
the scale uncertainty is not explicitly taken into account, and therefore following 
Ref.[^ the total error is increased to 0.0045. 



The HI collaboration have used their recent precise F2 data together with the 
BCDMS data |^ to obtain a precise determination of a<j as well as the gluon 
density 1T35[ (Sec ^XTj) : 



as{Mz) = 0.1150 ± 0.0017(expt)^[!;J^|]J^(model) , (HI, F2 96-97) (4.42) 

where the model error includes all uncertainties associated with the construc- 
tion of the QCD model used for the fit. In addition, a rather large theoretical 
uncertainty resulting from the renormalization and factorization scale choices 
should be added. This uncertainty amounts to about ±0.005 as is discussed in 



Refs.|T7|, |T75| 



At HERA, the strong coupling can also be determined in other ways. One 
such measurement is obtained by comparing the rates for (1 + 1) and (2 + 1) jetQ 
processes. A final state of (1 + 1) jet is produced at lowest order in in the ep 
scattering process, while a (2 + 1) jet final state is produced at next order in due 
to photon-gluon fusion (Fig.p3|) and QCD-Compton processes (Fig. |4.32| ). Unlike 

^^The notation refers to the proton remnant jet. 
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Figure 4.32: Diagrams of the QCD-Compton process in deep-inelastic lepton- 
proton scattering. 



the similar measurements from the e^e^ experiments, the ones in ep scattering 
are not at a fixed energy scale thus allowing the running of to be tested in a 
single experiment (see Fig. |4.33| below). The first determination by HI was based 
on the data from 1993 ||1 76|| . A similar determination by ZEUS was obtained 
177|| . In the meantime, the NLO QCD calculation have been 



from the 1994 data 

improved such that the prediction can be used to compare with experimental data 
in an extended kinematic phase space [p.78||. Using the new prediction, the HI 



collaboration has had an update using the 1994-1995 data [|179|] . The latest result 
was determined from the measured inclusive jet cross section dPa^ct/ dETdQ"^ 



m 



the range < < 2500 GeV^ and 150 < < 500O GeV^ using the 1995-1997 



e'^p data | |180|| . Jets were defined by the inclusive kx cluster algorithm | ]181[ | in 
the Breit frame|^ in the pseudorapidity range —1 < i] < 2.5 in the laboratory 
frame. The values determined in different Ex range are shown in Fig. 4.33| and 
are consistent with the scale dependence predicted by the renormalization group 
equation. The combined result choosing Et as the renormalization scale yields 

as{Mz) = 0.1186 ± 0.0030(exp) ± 0.0051(th) . (HI, jet 95-97) (4.43) 

The most recent result from ZEUS was derived from the measured dijet rate in 
the range 470 < < 20 000GeV^ using the 1996-1997 data p2| 



0.1161+°:°°^^(exp)+°:°°^Kth) . (ZEUS, dijet rate 96-97) (4.44) 



The dominant experimental error in both experiments is from the hadronic energy 
scale uncertainty, which is expected to be significantly reduced according to recent 
studies (Sec. ^.l.^) . The theoretical error includes uncertainties arising from scale 
choice, parton density functions, and hadronization correction. 

^"in the Breit frame, the gauge boson exchanged in the t channel is purely space-like with 
four momentum q = {0, 0, 0, —Q} and collides head-on with a parton from the proton. In the 
leading process (1 + 1 jet) the incoming quark is back-scattered and no transverse energy is 
produced. The high Et inclusive jet cross section receives only higher order contribution and 
is therefore directly sensitive to the strong coupling constant as- 



105 



a 0.26 
0.24 
0.22 

0.2 
0.18 
0.16 
0.14 
0.12 

0.1 



H1 



ttg from inclusive jet cross section 

for CTEQ5M1 parton densities 
150<Q^<5000 GeV^ 




10 



inclusive algo. 

• «s(Et) 
as(Mz) 



^^r=ET 



10 



Ej I GeV 



Figure 4.33: Values of as determined from the inclusive jet cross section. The 
results are shown for each Et value (circles). The single values are extrapolated 
to the mass (triangles) . The upper curves show the expected scale dependence 
and the size of the error band indicates the precision of the combined result. 



These results are shown in Fig. 4. 34 together with the new world average value 
as = 0.118 ± 0.002 and one of the most precise measurements from r lepton 



decays by ALEPH ||183| . The comparison shows that the results from the DIS ex- 
periments are consistent among themselves and also with the other independent 
measurement. Furthermore, the DIS measurements are now becoming increas- 
ingly precise, in particular in terms of the experimental uncertainties. The future 
high statistics HERA data are expected to improve further the experimental pre- 
cision. Similar improvement on the theoretical side is also expected when higher 
order QCD calculations for the structure functions and the jet rates will become 
fully available. 
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Figure 4.34: A comparison of a few selected measurements from the DIS data. 
Also shown are one of the most precise measurements from t lepton decays by 
ALEPH and the new world average from (the shaded band). The inner 
and full error bars represent respectively the experimental and total errors. See 
text for the label definitions and the references. 
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Chapter 5 



Neutral and charged current 
interactions at high 



In the previous chapter, the emphasis was put on the measurements of structure 
functions at relatively low and also low x, a new kinematic domain unexplored 
by fixed-target experiments, and their impact on the determination of parton (in 
particular gluon) density functions of the proton. As the integrated luminosity 
increases, the new kinematic domain at high up to a few 10 000 GeV^ has seen 
sufficient data to be studied for the first time. Contrary to low x HERA data, 
these data at high cover a range of x values which have been precisely mea- 
sured by fixed-target experiments. This allows the DGLAP evolution equations 
be tested over four orders of magnitude in . Deviations from the standard 
DIS expectations may be found as the largest value attainable at HERA is 
equivalent to a resolution power of one thousandth of the proton radius. In fact, 
already with the data from years in 1994 to 1996, HI and ZEUS have reported 
a low statistics excess of events at high and large x. In this chapter, the 
analysis on three independent data samples taken since 1994 will be discussed 
in some detail. The first e^j) data sample, corresponding to an integrated lumi- 
nosity of 35.6 pb~^ and a center-of-mass energy of 300 GeV (see Table p.2|) , was 
taken from 1994 to 1997. Both the e"p data of 1998-1999 and the e+p data of 
1999-2000 are taken at a center-of-mass energy of 320 GeV resulting from the 
increased proton energy of 920 GeV. The e~p data have an integrated luminosity 
of 16.4 pb~^. The recent higher energy e+p data of 45.9 pb~^ correspond to the 
data taken until the beginning of June, 2000.Q In Sec.BJ. a few relevant technical 



issues will be discussed first. The cross section results are given in Sec.B^ 



^The data taking continued till the beginning of September 2000 when a long upgrade 
shutdown has taken place. 
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5.1 Technical aspects 



The inclusive neutral (NC) and charged current (CC) cross section analysis on the 
later e~p data of 1998-1999 and e~^p data of 1999-2000 did not differ significantly 
from the one on the earlier e~^p data from 1994 to 1997 apart from a few technical 
improvements. Therefore in the following discussion, no separation is given unless 
stated otherwise. 



5.1.1 Event selections and background studies 

NC event selection 

The selection of NC events is only briefly described here as it did not differ in 



essence from the one on the earlier data of 1993 described in Sec. [4.3| and was based 
on the characteristic feature of the events; an identified scattered electron [ 184 
in the liquid argon (LAr) calorimeter and the hadronic final state which is mainly 
measured in the LAr calorimeter as well. In fact the main selection cuts were 
only slightly modified with respect to the one used in the earlier analysis: 

• A reconstructed event vertex |2;vtx — Zq\ < 35cm. As mentioned already 



in the previous analyses (Secs.OI and |4^), this requirement suppresses 



efficiently the non-ep background contributions. 

• An energy E'^ of the scattered electron measured in the LAr calorimeter 
greater than 11 GeV. Above this threshold the trigger was fully efficient 
and a small fraction of inefficient regions was excluded by applying fiducial 
cuts 1851, M . 

• An inelasticity ye lower than 0.9. This cut becomes more restrictive than 
the cut on E'^ for > 907.5 GeVl 

• A longitudinal momentum conservation verifying Eg + S/, greater than 
35 GeV. For an ideal detector without energy loss, the sum is expected 
to be 2Eo. 

The last three requirements minimize the size of radiative corrections and reduce 
the background from photoproduction. 

The dominant remaining background contribution was from photoproduction 
events where a hadronic final state faked the scattered electron. For the anal- 
ysis on the 1994-1997 e^p data, it was subtracted statistically according to the 
simulated photoproduction Monte Carlo events generated using the pythia ||187 



generator. The size of the subtraction was controlled with a subsample in which 
the genuine scattered electron is tagged in the small angle electron tagger at 
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Figure 5.1: A typical charged current event measured by the HI detector. 



z = —33 m away from the interaction point. This method suffers however from 
a large uncertainty (30%) on the estimated background events to be subtracted. 
For the analysis on the new data recorded since 1998, a new method was applied 
relying on the fact that the track charge associated to the electron candidate in 
the photoproduction event could be either negative or positive instead of being 
negative as expected from the e~p collisions. Therefore, the background can be 
subtracted with two times of positive charged tracks based on the data itself 
thereby reducing considerably the uncertainty in the subtraction. 



CC event selection 



Contrary to the NC selection, the selection of CC events represents a major 
challenge. The CC events are characterized by a missing transverse momentum 



^miss ^ p^^^ ^Yie undetected neutrino in the final state (see Fig. [5.1| for a 

typical CC event). The main selection cut was 

PT,h> 12 GeV, (5.1) 
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which was much lower than the cut pT,h > 25 GeV apphed in the earher stud- 
18^ , |T89| , |190|| based on the 1993-1994 data. While the kinematic coverage 



les 



has been significantly extended to the lower (Fig. |5.2|) , the acceptance to the 
background events were opened up at the same time. In addition to the pT,h cut 
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Figure 5.2: The extended kinematic region when the cut pT,h > 25 GeV is lowered 
to pT.h > 12 GeV. Also indicated are the HERA kinematic limit y = I and the 
actually analyzed y range 0.03 < y < 0.85. 

and the same vertex requirement as for the NC analysis, the CC analysis was 
limited in the inelasticity range 0.03 < y < 0.85 as also indicated in Fig. |5.2| in 
order that the kinematic quantities which can only be reconstructed from the 
hadronic final state are reasonably precise and that the trigger efficiency remains 
reasonably large (Sec. |5.1.^ ). 

Three background sources contribute in the remaining non-ep background 
events. These sources are halo muons, cosmic rays and beam-gas and beam- 
wall interactions. Since the rate of these background events are a few orders of 
magnitude higher than the CC rate (Table |2.4|) , there is a high probability that 
these events pile up among themselves or with other ep events. An example of 



such pile-up events is shown in Fig. 5.3 
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Figure 5.3: A pile-up event from three background sources: halo muon, cosmic 
ray, and beam-gas/wall interaction. 



One efficient rejection was achieved by applying the timing informations pro- 
vided by both the central jet chambers (CJC) and the LAr calorimeter. The CJC 
timing, tcjc? was determined from tracks which cross the sense wire planes, and 
has a resolution around 2 ns. The LAr timing ||191|| , tLAr, was determined from 



energy deposits in the trigger readout, which consists of 512 big towers[|, a coarser 
structure compared to the highly segmented cell granularity. The resolution of 
^LAr is about 10% of a bunch crossing (96 ns). 



The other powerful rejection used a set of topological filters qbgfmar ||93|] , 
which was based on the characteristic signature of cosmic and halo muon events. 
Part of the filters were developed for the earlier analysis described in Sec. [4.3| . The 
efficiency of the background filters was determined to be well above 95% based 
on the pseudo-CC sample (see Sec. ^.l.6[) and was in good agreement with the 
MC simulation. 

The residual background events were further reduced by other track based 
criteria. Two examples are 



small fraction of the big towers are formed 
BEMC/SPACAL and the forward PLUG calorimeter. 



from the backward calorimeter 
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• Track- Cluster link: For the pile-up events, the charged tracks contribut- 
ing to the event vertex have in general no correlation with the energy deposit 
in the LAr calorimeter. These events were efficiently suppressed by requir- 
ing a minimum distance of 0.5 in the <^ and pseudo-rapidity rj plane between 
the tracks and two largest energy deposits in transverse momentum. 

• Track multiplicity: The track multiplicity of a beam-gas or beam-wall 
event can be large. An important fraction of these tracks do not originate 
from a common vertex. These events were rejected by requiring that there 
were 10 or more tracks originating far away from the event vertex by more 
than 35 cm along —z direction. 

The final contamination of non-ep background events in the CC sample was de- 
termined by visual scanning to be 3.7% for the 1994-1997 e^p analysis and was 
improved by a factor of two for the new data and subsequently removed from 
the CC sample. The visual scanning technique as part of the final selection pro- 
cedure was important because it revealed not only often the need for further 
improvements in the background rejection but also occasionally interesting (ex- 
otic) events. In fact, the first events of observed isolated lepton events with large 
missing transverse momentum were found in this way. 

The dominant ep background sources originated from NC events and photo- 
production events due to the finite detector resolution and the limited geometrical 
acceptance (including the crack region in and z) . Nevertheless, in the transverse 
plane the energy flow is expected to be more isotropic in the photoproduction and 
NC events than in CC events. This is quantified by the ratio V of two variables 
Vap and Vp defined ||192|| as: 



Vap = -J2^H^^ for PT,-PT,h<0 (5.2) 



^P = E^^H-^ for PT,-PT,h>0 (5.3) 

standing respectively for the transverse energy flow antiparallel and parallel to the 
direction of the transverse momentum of the event. The sum in ( ^.2|) and (|5.3|) 
extends over all hadronic final state situated in the hemisphere which is opposite 
to or along PT,h- For the analysis on the 1994-1997 e~^p data, the background 
was significantly suppressed by requiring 

^=^<0.15. (5.4) 

Vp 

For the analysis of the new data, the photoproduction rejection was improved 
by introducing a new variable A(f)h,PLVG IP-93|| : 

A(f)h,pL\jG = 10/1 - 0plug| , (5.5) 
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where (j)h and 0plug are the angles of the hadronic final state measured respec- 
tively in the LAr calorimeter and in the forward PLUG calorimeter. While the 
^0h,PLUG distribution is nearly flat for the CC events, the photoproduction events 
are mainly distributed at large angles as illustrated in Fig. p.4| for Px^h < 25 GeV. 
In order to optimize the efficiency and the background contribution, a Pt,/i de- 
pendent cut for Pt,/i < 25 GeV is applied as 

180 - 90 

A0.,PLUG < 90 + y^p^^^_Q 2 ^y - 0-2) , (5.6) 

with 

and A(f)h,PLVG in degrees and Px^h in GeV. At the same time, the V cut is ex- 
tended from 0.15 to 0.2. Consequently the photoproduction contamination has 
been further suppressed by a factor of two for an efficiency which is hardly af- 
fected at Pr,h < 25 GeV and better at higher Pt,/i values (Fig. |5.5|) . The final 
photoproduction contribution, which was statistically subtracted from the data, 
was estimated to be 2.9% and 1.4% respectively for the 1994-1997 e~^p and 1998- 
1999 e~p analysis. The largest contribution is located at low and high y and 
did not exceed 1% in most of the kinematic region studied. 

Part of the NC DIS background contribution was suppressed by the V cut. 
The remaining contribution was further reduced by requiring that if an isolated 
lepton candidate was identified, neither the lepton in the LAr nor the associated 
track should satisfy |0 — > 120°, where and (j)h were respectively the 
angle of the lepton (or of the track) and of the hadronic final state. The final 
NC background was estimated for both e'^p and e~p analyses to be around 0.7% 
of the total selected samples. Again as for the photoproduction background the 
largest contribution was distributed at the low P^^h region and was negligible for 
most of the kinematic region under consideration. 

The final CC samples comprised about 700 events for 1994-1997 e~^p as well 
as 1998-1999 e~p analyses^, and 1000 events for the 1999-2000 e^p analysis. 

5.1.2 Alignment and electron angle measurement 

For the NC analysis, one of the key quantities needed for the kinematic reconstruc- 
tion is the polar angle of the scattered electrons. The polar angle is determined 
both by the trackers and the position measured in the high granularity calorimeter 
together with the event vertex. The two determinations are complementary: 

■^Thc integrated luminosity of the e'^p 1994-1997 data was about a factor of three larger than 
that of the e~p 1998-1999 data. On the other hand, the cross sections are the other way around 
so that the total selected number of events were about the same for the two data samples. 
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Figure 5.4: Distribution of A(f)h,pL\jG versus V (see ( |5.5|) and l \5.4 ) for definitions) 
for CC candidates in tie real data (large dots) and for the simulated photopro- 
duction Monte Carlo events (small dots ) wliich represent a factor of six larger 
integrated luminosity than in data. The full lines indicate the new cut applied 
for two example Pt^h values at 12 GeV and > 25 GeV^ in the e~p analysis while 
the dashed line for the old cut applied in the e~^p analysis. 
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Figure 5.5: Improvement on the CC selection efficiency and the suppression of 
the photoproduction contamination as a function of Px^h for the 1998-1999 e^p 
analysis with the V and A0/i^plug cuts (full lines) compared with the V cut only 
(dashed lines) used in the 1994-1997 e'^p analysis. 
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• When the trackers are fully operational, the former determination gives 
a better precision than the latter one. On the other hand, it suffers from 
occasional inefficiencies of part of the tracking system^ and needs a delicated 
modeling of the fraction of inefficient regions and its time dependence in 
the Monte Carlo simulation. 

• The latter determination though intrinsically less precise is more stable as 
a function of time. However, it is achieved when the calorimeter is well 
aligned to the tracking system. 

The alignment was performed in this analysis based on a subsample where 
the track associated to the scattered electron was well measured by the central 
tracking device (i.e. both the CJC and the inner and outer z chambers were op- 
erational). Under such condition, the polar angle of the track ^trk gave a good 
reference for the angle of the scattered electron Qf.- This is illustrated by Fig.^ 
where the angle 6'trk is found in good agreement with 6'e,gen, the generated angle of 
the scattered electron, independent of the impact position Ztrk of the track on the 
front of the calorimeter. An octagon around the z axis with a minimum radius 
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Figure 5.6: The polar angle difference between the generated scattered electron 
and the associated track as a function of the z impact of the track. 



of 105 cm was taken to approximate the front plane of the calorimeter. The track 
impact point was defined on this octagon after having been extrapolated from 
the event vertex {x^tx, 2/vtx, z^tx) taking into account its curvature measured by 
the tracker within the magnetic field. The impact position of the correspond- 
ing cluster measured in the LAr calorimeter on the same octagon was obtained. 



^In order to maximize the statistical precision of the cross section measurements, which is 
still limited at high , the data sample used includes data in which part of the trackers could 
be temporarily off. 
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Table 5.1: Parameters describing the relative alignment between the LAr 
calorimeter and the central tracking system for the 1997 e~^p, 1998-1999 e~p, 
and 1999-2000 e+p data. 



Year 


XLAr (mm) 


j/LAr (mm) 


ZLAr (mm) 


a (mrad) 


f3 (mrad) 


7 (mrad) 


1997 e+p 

98- 99 e^p 

99- 00 e+p 


-0.5 ±0.2 
0.6 ±0.2 
0.9 ±0.1 


-1.5±0.2 
-0.5 ±0.2 
0.4±0.1 


-4.2 ± .08 
-1.7± .07 
-2.0 ± .03 


0.6 ±0.1 
-0.5 ±0.1 
-0.6 ± .05 


-1.7±0.1 
-0.0 ±0.1 
-0.2 ± .05 


0.7 ± .06 
2.4 ± .06 
0.9 ± .03 



however, with a straight hne between the barycenter of the cluster and the event 
vertex. 

Defining (xlar, Vlat, -^LAr) as the center of the LAr calorimeter in the coordi- 
nator system of CJC with additional parameters a, (3, and 7 being respectively 
for anti-clock-wise rotations around the x, y and z axes, the cluster barycen- 
ter {xe,ye,Ze) of the scattered electron measured in the LAr calorimeter can be 
written in the CJC system as[| 

= Xe- XLAr + - 7?/e , 

y'^ = ye- 2/LAr " O^e + , (5.8) 

Z'^ = Ze- ZLAr + aye - • 



The parameters were determined by minimizing the sum of the difference 
fltrk — ci'e with a = X, y and z. The resulting parameters are shown in Table 



5.1 for the 1997 e+p, 1998-1999 e'p, and 1999-2000 e+p data. As an example. 



Fig. [5.7| shows how the alignment on the 1997 e+p data improves the difference 
Otrk — Before the alignment, there was a significant rotation around the y axis 
which resulted in a tilt in the x — z plane. Around the other axes, there were no 
important rotation, but the shifts in y and in particular in z were sizable. As far 
as the measurement of the scattering angle 9^ is concerned, the quantity Ztvk — z'^ 



e 



is most relevant, which as shown in Fig.5.7 is in better agreement with the Monte 



Carlo simulation after the alignment. The overall shape, which stems partly from 
the pad structure and partly from the energy distribution on the pads from which 
the barycenter can be defined, cannot of course be improved by the alignment. 
But this effect can be and has been corrected for based on the Monte Carlo. 

For the 1994-1997 e+p data analysis, the final angle of the scattered electron 
was determined from the track when it was well measured by the central tracking 
system [Qf, > 35°) and in other cases from cluster measured in the calorimeter. 



^ Since the rotation angles are expected to be small, the following approximations cosx ~ 1 
and sinx ~ x can thus be made. The z component Ze in data should be replaced by z^°^'^ which 
is related to Ze by z'^°^^ = 23.67 + {ze - 23.67)(1.0 - 0.0027). This correction arises because 
during the final step of data processing the LAr geometrical parameters were mistakenly taken 
to correspond to a warm environment instead of the cold LAr in reality. 
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Figure 5.7: The difference of the impact point x, y and z between track and 
cluster as a function of z impact of the track before and after the alignment. 
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The fraction of latter determination was about 40% and 100% respectively for 
9(. > 35° and 6(. < 35°. The corresponding systematic uncertainty from the 
tracker and the calorimeter was respectively 1 mrad and 3 mrad. For the analyses 
of other data samples, the angle was determined from the cluster position in 
the calorimeter together with the event vertex due to a reduced efficiency of the 
central trackers. 



5.1.3 Electron energy measurement 

For the cross section measurements to be shown below, the energy measurement 
of the scattered electron has been improved by an in situ calibration by using the 
increased event sample and the overconstrained kinematic reconstruction. 

Several complementary methods and event samples have been used: 

• The double angle method based on NC DIS events (DIS DA): 

As shown in (|4.4| ), the energy of the scattered electron can be predicted 
from the polar angles of the electron 6e and of the hadronic system 6^. The 
resolution of 9h is sufficiently good when one uses the subsample at y-£ < 0.3 
{ys < 0.5) for 80° < < 153° (40° < < 80°). In this region, the event 
statistics is large enough to improve the energy measurement locally in 
finely segmented z and (p grid defined by the impact position of the electron 
track on the LAr calorimeter. 

• The elastic QED Compton events and exclusive two photon e+e^ 
pair production (QED Compton/e"*"e^): Due to the limited number of 
events from NC DIS process in the forward region {6e < 40°), two additional 
event samples are used to check and obtain the calibration constant together 
with the DIS sample. These samples are the elastic QED Compton events 
and the exclusive two photon e~^e~ pair production. Contrary to the DIS 
events, these samples provide two leptonic final states for the calibration 
when both are measured in the LAr calorimeter. Furthermore since these 
events cover a large energy range, they allow the energy linearity be studied. 

• The uj method based on NC DIS events (u; DIS): This is an approx- 
imate kinematic method 111 9411 which assumes that the relative error in the 



measured T^h and Pt,/i is the same (i.e. (5S/i/E/i = SPx^h/ PT,h)- This method 
was also based on the NC DIS sample and therefore it was not independent 
of the DA method. On the other hand, the u method is less sensitive to 
the effects of initial state QED radiation than the DA method. 

The resulting electromagnetic energy scale as determined for the 1994-1997 



e^p data analysis is shown in Fig. 5.8 together with the quoted systematic uncer- 



tainty as indicated by the error band. 
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Figure 5.8: Comparison of the electromagnetic energy scale as determined by 
different metliods and event samples. Sliown is the mean fractional energy shift 
of the different method from the absolute energy scale, {SE'^/E'^). The shaded 
error band shows the systematic uncertainty on the energy scale quoted for the 
cross section measurement, which varies from 0.7% to 3%, depending on the 
location of the electron in the detector. 
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For the analyses on the new data, the energy scale was checked using the newly 
available event sample. Small variation within the quoted systematic uncertainty 
was found and taken into account in the new analysis. 



5.1.4 Hadronic energy measurement 

A precise hadronic energy measurement is crucial for all measurements relying 
on the hadronic final state for the kinematic reconstruction. In fact many of 
these measurements are now limited by the uncertainty of the hadronic energy 
scale, which was estimated to be about 4%. For instance, such an uncertainty 
dominates the experimental error of the as(M|) measurement of as{M^) based 
on the inclusive jet cross section ||180|| (Sec.| 



For the cross section measurements presented here, a detailed study was car- 



ried out ||195|| in order to improve the uncertainty. While the precision of the 
electron energy calibration was limited in the forward region due to the lim- 
ited statistics, the hadron energy calibration does not have such a limitation as 
the hadronic final state covers rather uniformly all the calorimeter region. The 
calibration method is based on the transverse energy conservation, namely the 
transverse momentum of the hadronic final state Pt^h should be equal to that 
of the scattered electron P^ e, the latter can be again predicted with the double 
angle method as Pt,da- 

A clean and unbiased calibration sample is defined as follows: 

• An identified scattered electron having e > 10 GeV, 

• The predicted Pe,DA is not strongly affected by the initial state QED radia- 
tive events by requiring^ Pr,e/-PT,DA > 0.88, 



• One and only one reconstructed jet based on the cone algorithm ||196|| with 
a cone radius of 1 rad in the and pseudorapidity rj plane and a minimum 
transverse momentum of 4 GeV. In addition, in order that the calibration 
factor to be determined is relevant for the LAr calorimeter, the jet has 
to be well contained in the LAr calorimeter by requiring that there is no 
more than 1% energy or transverse energy deposit in either the backward 
calorimeter or the tail-catcher. The jet is also required to be well away from 
the beam pipe in the forward direction (6'jct > 7°) so that the calibration 
would not correct for the beam energy loss. 

The calibration is then performed in two steps. In the first step, the rela- 
tive difference in PT,h/ Pt,ba between data and the Monte Carlo simulation is 

^As mentioned previously, a cut on Eg + E/i is also very efficient in suppressing the radiative 
effect. On the other hand, since the sum includes a contribution which is directly proportional 
to the energy scale of the hadronic final state, such a cut may bias the calibration. 
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determined in finely segmented and O'-^ctfl S^i*^ account for possible detector 
effects which may be year dependent. In the second step, an absolute hadronic 
energy scale is achieved both in data and simulation by imposing the transverse 
momentum balance between the scattered electron and the hadronic final state. 



Fig. |5.9| shows the performance of the calibration based on a large inclusive 
event sample which was selected with less restrictive cuts as was used for defining 
the calibration sample. Before the calibration, the hadronic scale was found 
to be significantly lower than expected and there was in addition a strong Px^h 
dependence. The difference between data and the simulation was about 4%. After 
the calibration, a precision of 2% was achieved for the first time in HI both on 
the hadronic scale and the systematic difference between data and simulation. 

The performance of the calibration was further checked on a few other aspects: 

• Hadronization model dependence: The check concerns whether there is 
any dependence of the calibration on the usage of two different hadroniza- 
tion models in the Monte Carlo event simulations. One model uses the 
color dipole model in its ARIADNE implementation and the other model 
uses the matrix element plus parton shower model (MEPS) as implemented 
in LEPTO. While some differences were found in each calibration step, the 
overall calibration procedure was found to be independent of the Monte 
Carlo models as expected because of the imposed balance in transverse 
momentum. 

• Independent event sample: While the calibration correction factors were 
obtained using the clean one jet sample, the performance of the calibration 
was checked with an independent sample in which two or more jets are 
reconstructed. 

• Other jet reconstruction algorithms: The calibration factors obtained 
from the cone algorithm were applied to jets which were reconstructed with 
two other jet algorithms to check the dependence on the jet reconstruction. 

In all cases, the results obtained were within the quoted systematic uncertainty. 
In addition, the absolute energy calibration also improved the relative energy 



resolution a{PT,h)/PT,h by 5% at Pt - 35 GeV and 15% at Pt ~ 12 GeV |95 . 

The reduced systematic uncertainty on the hadronic energy scale has an im- 
portant impact on analyses relying on the hadronic final state. One example 
is the measurement of the charged current cross sections. This is illustrated in 
Fig. p.lO| (a). The magnitude and turn-over behavior at around 2000 GeV^ depend 



2 



on the selected event distribution shown in Fig. ^.lO| (b). At the intermediate Q 

''The quantity ^j;,^ differs from 6jct in that the former is calculated always with a fixed vertex 
position at zero instead of with the measured event vertex as in the latter. In this way, the 
angle 6'^^^ is well defined independent of the vertex position. 
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Figure 5.9: The transverse momentum of tie liadronic final state, PT,h, over that 
of the scattered electron, Pt,dA; as predicted with the double angle method, as 
function of Pt.da 3,nd the inclusive angle of the hadronic final state 6^ before and 
after the calibration. 
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Figure 5.10: (a) The resulting uncertainties on the measurement of the e^p CC 
cross section from the improved hadronic energy scale uncertainty of ±2 % (dark 
sharded area enclosed with full points ) in comparison with an uncertainty of ±4 % 
(light shared area enclosed with open points), (b) The selected event distribution 
as a function of upon which the form of the uncertainty depends. 
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the event distribution is rather flat and the energy scale uncertainty on the cross 
section is smallest, as the event migration from one side compensates the mi- 
gration from the other side. At lower the drop in the event distribution is 
completely due to the the selection cuts, while at higher the selection efficiency 
is very high and therefore the shape of the event distribution is intrinsically due 
to the falling cross section as increases. 

5.1.5 Calorimetric noise suppression 

One of the problems in the analyses relying on the liquid argon calorimeter (LAr) 
is the distorted measurement of the kinematic variables at the low y region [|195|| . 
The problem is related to the fact that at low y < 0.05, hadrons are produced 
in the forward direction. The z component Pz,h is comparable with the energy 
Eh thereby resulting in a small value of T^h = Eh — Pz,h{^ 2.76 GeV). Any 
background energy located in the barrel and backward calorimeter if taken as 
part of the hadronic final state will thus distort significantly the total S/j (the 
more backward a given energy Eh is deposited, the larger the resulting Eh —Pz,h)- 
Background energies arise primarily from the following sources: 

• Electronic noise: the noise inherent in the electronics of the calorimeter 
readout, 

• Backscattering: the contribution from the secondary scattering of final 
state particles, 

• Pile-up events: the contribution from the accidental coincidence of a low 
y event with a background process, e.g. a beam halo or a cosmic ray event,^ 

• Hadron/photon separation: the large angle radiative photons measured 
in the LAr are misinterpreted as part of the hadronic final state if they are 
not explicitly identified, 

• Hadron/electron separation: for neutral current events, a wrong sep- 
aration between the scattered electron and the hadronic final state can 
result in a biased measurement for the Eh — Pz,h measurement. This could 
happen when the scattered electron is misidentified. The imperfect clus- 
ter algorithm can give rise to multiple clusters for the scattered electron in 
particular when it hits a (j) crack between two octants, or a z crack between 
two wheels, or the overlapping region between the BBE and the backward 
calorimeter SPACAL. In such a case, part of the electron energy may be 
wrongly attributed to hadrons. Since the electron is located predominantly 

^The LAr time sensitivity of several microseconds allows energy deposits of these background 
events to pile up on real physics events. 
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at large polar angle in the backward region of the calorimeter, a fraction 
of its energy if assigned to hadrons may dominate over the small S/^ in the 
forward region, thus bias the hadronic measurement. 



The last source, which affects only neutral current events, can be eliminated by 
improving either the current electron identification programs or the clustering 
algorithm. The radiative contribution is more difficult to deal with for an in- 
clusive analysis. In order to study the other noise contributions, these last two 
contributions are explicitly removed or suppressed from Monte Carlo samples 
used. 

Part of the electronic noise has been suppressed during the online and the 
offline reconstruction ||197|| . The remaining background contribution was further 



suppressed in the past by various noise suppression algorithms. The basic idea of 
one of these algorithms, to be called "old suppression" | |198|| in the following, is to 



identify and subsequently suppress isolated low energy deposits (clusters) in the 
LAr calorimeter. After the suppression, a significant improvement is achieved in 
particular for y down to ~ 0.05 (Fig. |5.11| (a)). For the region at smaller there 
is still an important bias to the measured yn- 

A closer look at the energy distribution as a function of dh (Fig. |5.1l| (b)) for 
?/ < 0.1 shows that there is an energy surplus towards the large angles. A detailed 
study ||199|| revealed that two distinct types of background energies contribute: 



energetic deposits affect a small fraction of the sample and less energetic ones 
affect the majority of the events. The first type originates from the accidental 
coincidence of a D IS event with a halo or a cosmic muon event. These background 
muons can produce rather energetic showers resulting in very large distortions 
in the measurement of yt (up to an order of magnitude with respect to ytmc)- 
Two specific algorithms were thus developed to identify the background energy 
patterns and subsequently suppress them from the energy measurement. The 
second type is due to other contributions such as the residual electronic noise 
and the backscattering contribution. To suppress these background energies, a 
higher energy threshold was found necessary [ |199|] . After the "new suppression" , 
the improvement is clearly demonstrated not only on the average value of the 
measured yh (Fig. ^.ll] (a)) but also in terms of the energy distribution over the 
detector (Fig. p.ll| (b)). It should be pointed out that the measured and yh 
do not have to coincide with the generated quantities as the absolute calibration 
discussed in the previous section is not applied for this comparison. At very low 
the energy loss in the beam pipe becomes increasingly important (see the energy 
distributions at small angles in Fig. |5.11| (b)), one expects y^ to be smaller than 

I/true • 



It has been checked |p.99|| based on the simulated MC samples that what has 



been suppressed were indeed the various noise contributions discussed above and 
there was no evidence that signal energies were affected. The suppressed noise 
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Figure 5.11: A comparison of the measured yh as a function of y (a), and of the 
measured Eh as a function of 9^ at low ?/ < 0.1 (h) based on the neutral current 
MC files. 



in the data is finally compared with that in the simulation using the NC e~p 
data taken in the years from 1998 to 1999 as an example. First in Fig. ^.l2| (a), 
the event distribution as a function of Snoiso/2-E'e (Snoise = -Enoiso — Pz.noisc) is 
compared. In Fig. |5.12| (b), the relative contribution of the suppressed Snoiso to 
the total measured S including Snoise is shown as a function of yh, which does 
not include the contribution from the suppressed noise. The shaded bands show 
the quoted systematic uncertainty (±25%) on the noise suppression. Within the 
uncertainty quoted, the data are well described by the simulation. 



5.1.6 Trigger efficiency for charged current events 

The trigger efficiency for CC events is one of the most critical parts among the 
various efficiencies needed for the cross section measurements as it was not auto- 
matically simulated in the Monte Carlo production. 
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Figure 5.12: Comparison of (a) the event distribution as a function of Snoise/2-E'e, 
(b) the relative contribution of the suppressed Snoisc to the total measured E(= 
S/j + Snoisc) as a function of yh. The data (points) are the NC events from 
the e~p data taken in 1998-1999, while the histograms show the corresponding 
Monte Carlo simulation. The shaded error bands represent the quoted systematic 
uncertainty ('±25%j on the suppressed noise contribution. 



The logics designed for triggering CC events are based mainly on the missing 
transverse energies provided by the LAr calorimeter and on the timing informa- 
tion from either the proportional chambers or the LAr calorimeter. From the 
detector point of view, the hadronic final state of a CC event is the same as that 
of a NC event. Therefore, the CC trigger efficiency as well as other CC efficien- 
cies can be determined from high statistics NC events from real data when all 
information associated to the scattered electron is removed. This is the so called 
pseudo-CC sample^ mentioned already in Sec. |5.1.1| . 

The CC trigger efficiencies for three different run periods are compared in 
Fig. ^.l3| as functions of the missing transverse momentum P^m. and the kine- 
matic variable yh- The lower cuts on these two quantities in the analysis are 
respectively 12 GeV and 0.03 (Sec. |5.1.T| ). This ensures that the trigger efficiency 
is always above ~ 50%. During the shutdown between 1997 and 1998, part of the 
preamplifiers with large capacity were upgraded by low-noise amplifiers. With 
this upgrade | |200| | and other hardware modifications | |20 1| | , the trigger efficiency 
is expected to be improved as indeed seen at the low y region. Unfortunately, 
during the first part (~ 64%) of the e~p data taking period, the missing trans- 
verse energies provided from the LAr trigger towers were modified to correspond 
to an earlier bunch crossing instead of the nominal one thereby resulting in a 

^These events were further reweighted to the CC cross section such that they behave kine- 
matically just hke real CC events. 
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Figure 5.13: Comparisons of CC trigger efEciencies from tliree different run pe- 
riods as a function of (a) the missing transverse momentum Pf^h and (h) the 
kinematical variable Uh- The selection cuts Pt^h > 12 GeV and 0.03 < Uh < 0.85 
are indicated with the dashed lines. 



significant efficiency loss. The problem was fixed in the second part of the data 
taking. The efficiency shown in Fig. |5.13| for the e~p data is the averaged value 
of the two parts. The inefficiency at the low y region occurs because two big 
towers close to the forward beam pipe are not considered so far for the sum in 
the missing transverse energies. Currently there has been efforts investigating 
the possibility of reopening these trigger towers to improve the situation. 



5.2 Measurement of inclusive cross sections at 
high 

In this section, many cross section results will be given and compared with the 
standard DIS expectations. The cross sections are obtained following the mea- 



surement procedure shown in Sec. [4.3.4 The results are presented in the form of 
reduced cross sections, which are related with the double differential cross sections 
(Eq.(|;9|)) as: 

where o"^ stand for cross sections respectively for e+p and e~p collisions and 
Y± = 1 ± {1 — are the helicity functions. 

Similarly the CC reduced cross sections are related with the double differential 
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cross sections (Eq.( p.l^ ))P'| 

and have the following simple quark flavor decomposition when the LO structure 
function formula in Eqs.( ^.2^ )-( ^.29| ) are used: 

a+c = x[{u + c) + {l- yf {d + s)] (5.11) 
^cc = ^ [(^ + c) + (1 - y)' (rf + s)] . (5.12) 

Before presenting these results, various systematic sources and their effects on 
the cross sections will be discussed. 



5.2.1 Systematic sources and resulting uncertainties on 
the cross section measurements 

The uncertainties in the measurement lead to systematic errors on the cross 
sections which can be point to point correlated or uncorrelated. All the correlated 
systematic errors were checked to be symmetric to a good approximation and 
are assumed so in the following. The correlated systematic errors and main 
uncorrelated systematic errors of the NC and CC cross sections are given below: 

• An uncertainty of the energy of the scattered electron is (see Fig. |5.8| ) 1% 
if the z position]^ of its impact on the calorimeter is in the backward part 
{z < -145cm), 0.7% in the CBl and CB2 wheels (-145 < z < 20cm), 
1.5% for 20 < z < 100cm and 3% in the forward part {z > 100cm). 
These uncertainties are obtained by the quadratic sum of an uncorrelated 
uncertainty and a point to point correlated uncertainty. This correlated 
uncertainty comes mainly from the potential bias of the calibration method 
and is estimated to be 0.5% in the whole LAr calorimeter. The resulting 
correlated (uncorrelated) systematic error on the NC cross sections is < 



3 (5)% except for the measurement at the two highest x values Due 
to the smaller luminosity for the 1998-1999 e~p data, the uncertainty in 
the CBl and CB2 wheels was verified to 1% with the same uncertainty for 



the rest of the calorimeter region |204|. For the high statistics e~^p data 



of 1999-2000, the uncertainty of the electron energy scale is expected to 
be improved in particular in the forward region. For the preliminary cross 
section results which were based on about 2/3 of the full data sample, the 
uncertainty was conservatively quoted to be the same as the e~p data [|202 



^"in the considered kinematic range {Q^ > 300 GeV ), the proton mass term AP can be 
neglected. 

-'^^The variation as a function of z is mainly a reflection of the statistical precision of the data 
samples with which the uncertainty could be checked. 
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A correlated uncertainty of 1 mrad on the electron polar angle, and an 
uncorrelated uncertainty of 2.8 mrad when the angle is determined with the 
position measured in the calorimeter and the event vertex. The resulting 
correlated (uncorrelated) systematic error is small, typically < 1(2)%. 

An uncertainty of 2% on the hadronic energy in the LAr calorimeter which 
is obtained from the quadratic sum of an uncorrelated systematic uncer- 
tainty of 1.7% and a correlated one of 1% originating from the calibration 
method and from the uncertainty of the reference scale of the scattered elec- 
tron. The resulting correlated systematic error increases at low y, and is 
typically < 4% except at high for the CC measurements (see Fig. |5.10| ). 



For the preliminary results of the 1999-2000 e^p data, a conservative un- 
certainty of 3% was quoted. It should be noted that the same uncertainty 
on the hadronic energy measurement can result in different uncertainties 
on the cross section measurements. For CC cross sections at high Q^, it is 
larger for the e^p data at ^/s ~ 300 GeV and smaller for the e~p data at 
~ 320 GeV due to different Q"^ dependences of the cross sections (see 



Sec.5.2. 



An uncertainty of 7(3)% on the energy of the hadronic final state measured 
in the SPACAL (tracking systemQ). Their influence on the cross sections 
is small compared to the uncorrelated uncertainty of the LAr calorimeter 
energy. 

A correlated uncertainty of 25% on the energy identified as noise in the 



LAr calorimeter (see Sec. |5.1.5| ). The resulting systematic error is largest 
at low y, reaching 10 - 15% at x = 0.65 and < 2000 GeV^ in the NC 
measurements, while it remains below 5% for the CC measurements. 

A variation of the anti-photoproduction cuts (see Sec. |5.1.1] ). The resulting 



correlated systematic error reaches a maximum of 12% at low x and in 
the CC analyses. 

An uncertainty of 30% on the subtracted photoproduction background. The 
resulting correlated error is always smaller than 5% in the NC and CC 
analyses and is further reduced in the analyses of the 1998-1999 e~p and 
1999-2000 e^p data due to the improved subtraction method for NC and 
the additional anti-photoproduction cut for CC. 



The other considered uncertainties giving rise to uncorrelated systematic errors 
on the cross sections are: 

^^An improvement in the energy resolution of about 10 — 20%, for events having a Pt,/i 
between 10 to 25 GeV, is obtained by using a combination of the momentum of low transverse 
momentum particles [Pt < 2 GcV) measured in the central tracking detector with the energy 
deposited by other particles of the hadronic final state measured in the calorimeter. 
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An error of 2% (4% aX y > 0.5 and < 599 QeV^) from the electron 
identification efficiency in the NC analyses. 

an error of 1% from the efficiency of the track-cluster link requirement in 
the NC analyses. 

An error of 0.5 (3 — 8)% from the trigger efficiency in the NC (CC) analyses. 

An uncertainty of 1(3)% from the QED radiative corrections in the NC 
(CC) analyses p5| . 



• An error of 3% from the efficiency of the non-ep background filters in the 
CC analyses. 

• An error of 2% (5% for y < 0.1) from the efficiency of the event vertex 
reconstruction in the CC analyses. 

The typical total systematic error for the NC (CC) double differential cross sec- 
tion is about 4(8)%. In addition a normalization error of 1.5%, 1.8%, and 1.7% 
respectively for the 1994-1997 e+p data, the 1998-1999 e'p data and the 1999- 
2000 e^p data has to be taken into account, which is not included in the results 
shown in the following subsections. 



5.2.2 Measurement and comparison of NC and CC re- 
duced cross sections in e^p collisions at two center- 
of-mass energies 

The NC reduced cross sections obtained from 35.6 pb~^ of the e^p data taken 
from 1994 to 1997 at the center-of-mass energy of 300 GeV has recently been 
published by HI |Q. The measurement covers a kinematic region for between 



150 and 30 000GeV2, for x between 0.0032 and 0.65, and for y between 0.007 
and 0.88. This kinematic range has significantly extended the previous HERA 
measurements ||131|| both in (from 5000 to 30 000 GeV^) and towards higher 



X, with measurements at x = 0.65 for between 650 and 20 000GeV^. At 
Q2 < 500 GeV^ the total error is dominated by the systematic uncertainties in 
the energy scale and identification efficiency of the scattered electron and by 
the uncertainty in the energy scale of the hadronic final state. In this region 
the systematic error is typically 4%. At higher the statistical error becomes 
increasingly dominant. 

The new e^p data of 45.9 pb~^ at a higher center-of-mass energy of 320 GeV 
has been analyzed and the preliminary results have been made available for the 
summer conferences [|202|| . The new data with slightly improved statistical preci- 



sion cover the same kinematic range as the 1994-1997 e+p data. The two indepen- 



dent measurements are compared in Fig. 5. 14 and found in good agreement. The 
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Figure 5.14: The e~^p NC reduced cross section ajyc is compared to the NLO 
QCD Fit Shown are the measurements from 1994-1999 data and from 
1999-2000 data M/ and the NLO QCD Fit predictions for different center-of- 
mass energies. The inner error bars represent the statistical error, and the outer 
error bars show the total error. 
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cross sections rise strongly as x decreases as seen in the low data, which can be 
interpreted as the increasing contributions from the gluon and sea quarks at low 
X. The measurements are well described by the standard DIS expectation from a 
NLO QCD Fit |3 to the e+p 1994-1997 dataQ as well as HI structure function 
data at low and fixed-target data by NMC |33l and BCDMS fl^. The 



difference in the center-of-mass energies is also indicated but remains small. 



The similar comparison for the CC reduced cross sections is shown in Fig. |5.15 
The e~^p 1994-1997 double differential CC cross sections cover the kinematic re- 
gion for between 300 and ISOOOGeV^ and for x between 0.013 and 0.4 |94 



The uncertainties on the measurements are dominated by the statistical errors. 
The largest systematic errors come from the uncertainty on the energy scale of 
the hadronic final state at high and from the uncertainty on the trigger effi- 
ciency at low and high x. A similar measurement has been published by the 
ZEUS experiment ||203|1 . The new e'^p data p02| , covering a similar range, are in 



good agreement with the published data and both measurements agree also well 
with the Standard Model (SM) expectations based on the same NLO QCD Fit 
mentioned above. 

Given the good agreement between the two data sets and in order to facilitate 
the comparison of the e'^p data with the e~p date in the later sections, we have 
combined the e~^p data at different center-of-mass energies to the higher energy 
in the following way^ 

(5.13) 



-^300 {(^i^SOo/ ^i,32o) + ^ 



320 



where o"™^^ and c"*^ are respectively the measured and expected cross section. 
The statistical error is determined correspondingly and the systematic error is 
assumed to be fully correlated between the two measurements and taken from 
the new measurement. 



^^The agreement is equally good when the input data used in the fit was limited to < 
150 GeV^ @. The fit is any case independent of the new e~^p and e p data taken since 1998. 

^^The alternative way would be to make the correction to the measured cross section at 
300 GeV from the 1994-1997 e+p data and to combine the corrected cross section with the one 
at 320 GeV from the 1999-2000 e+p data with the standard combination method using error 
weights. Indeed, the combined results obtained from both methods agree in the low region 
where the data are precise. However, the combined results can be substantially different in the 
region of higher and high x where the statistical error of the data is large (> 30%). The 
method using the luminosity weights is believed to be less sensitive to the statistical fluctuation. 
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Figure 5.15: The CC e~^p reduced cross section gqc compared to the NLO QCD 
Fit J9^/ . Shown are the measurements from 1994-1999 data and from 1999- 
2000 data and the NLO QCD Fit predictions respectively for ^ = 320 GeV 
(full curves) and 300 GeV (dashed curves). The inner error bars represent the 
statistical error, and the outer error bars show the total error. 
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5.2.3 Measurement and comparison of NC and CC re- 
duced cross sections in e^p and e j9 collisions 

The combined e^p NC reduced cross sections are compared with the resuhs from 
the e-p 1998-1999 data |0l in Fig.|5l6|. For < 1200 GeV^ where the cross 



section is dominated by the contribution from one-photon exchange, the data 
are essentially independent of the lepton-beam charges. With the increasing 
the e~p cross section becomes larger than the e^p cross section. The data thus 
give the clear evidence for the presence of the structure function xF^ (Eg. ( |5.9| )). 
The standard DIS expectations from the NLO QCD Fit [Q, which was partly 
constrained by the 1994-1997 e^p cross section data, are compared with the 
measured cross sections. The fit is found to give a good description of the x, 
behavior of the data, though in some of the phase space it has a slight tendency 
to be lower than the e~p cross sections. 

The combined e^p CC reduced cross sections are compared with the results of 
the e-p 1998-1999 data in Fig.^^l^. Contrary to the NC data, the CC e-p 



data are everywhere different from the e'^p data. This is the particular feature of 
the CC interaction mentioned already in Sec. |3.1.2] . At low Q^, the difference is less 



pronounced because the sea quark contribution is relatively important, whereas 
at high and high x, the difference can mostly be attributed to the underlying 
quark flavor difference between xu and (1 — y^xd as indicated respectively with 
the dashed and dash-pointed lines in Fig. |5.17| . Therefore the e~p CC cross section 
at high and high x can be used to constrain the u valence quark while the e~^p 
CC cross section provides the constraint for the d valence quark (see Sec . [5. 31 ). The 
expectations from the same NLO QCD Fit are found to give a good description 
of the measured CC cross sections. 



5.2.4 Measurement of the structure function xF^ at high 

Using the measured e~^p and e~p double differential NC cross sections, both 
ZEUS [|206|] and HI ||204|| have measured the structure function^ xF3{x,Q'^) at 
the high region. The measurement was obtained using Eq. ( |5.9|) by subtracting 
the e~^p cross section from the e~p cross section. Since the data are taken at 
slight different center-of-mass energies, a small correction|^ from the longitudinal 

^^It should be pointed out that the structure function xF^{x,Q^) measured at HERA arises 
from the contributions of the jZ'^ interference and Z° exchange with the dominant contribution 



from the former (Eq.3.25), in which the weak neutral current couplings enter. Therefore it is 



different from the corresponding structure function xF^ = x{uy + dy), the sum of the u and d 



valence quarks (Eq.(3.54)), measured by neutrino-fixed-target experiments. 

I'^The correction, estimated with Fl from the NLO QCD Fit Q, is about 10% at the lowest 
X and negligible elsewhere. 
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Figure 5.16: The combined e~^p NC reduced cross section a^c from 1994-2000 is 
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tations from the NLO QCD Fit. The inner error bars represent the statistical 
error, and the outer error bars show the total error. 
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structure function is needed: 
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(5.14) 



where the subscripts 300 and 320 denote the center-of-mass energies. The HI re- 
sults are shown in Fig. 5.18 as a function of x at three values of at 1500 GeV^, 



5000 GeV^, and 12 000GeV^. The data are in good agreement with the expecta- 
tion based on the HI NLO QCD Fit ||9^ with a peak at x ^ 0.1 reflecting the 
valence quark structure of the constituent quarks. The ZEUS results are shown 
in Fig. |5.19| as a function of for two values of x at 0.56 and 0.1. The expecta- 
tion of the SM, evaluated with both the CTEQ4 |0|| and the MRST |5|] parton 
distributions, gives a good description of the data. 



5.2.5 Helicity structure of the CC cross sections 

The measured double differential CC cross sections can be used to test the helicity 
dependence of the electron-quark interaction in the region of approximate Bjorken 
scaling, x ~ 0.1. In Fig. |5.2U| , the reduced cross section^ from the combined e^p 
1994-2000 data and the e~p 1998-1999 data are shown as a function of (1 — 
yY. This kinematic variable is directly related to the scattering angle 6* in the 
electron-quark center-of-mass system through cos^(^*/2) = 1 — y. The e~^p and 
e~p data show a different behavior. This is consistent with the expectations 
from the NLO QCD Fit 0. Indeed, from Eqs.(pll) and (|53^ ), one expects a 
small isotropic distribution from positron-antiquark scattering {u, c) and a strong 
angular dependence on (1 — y)^ from positron-quark scattering {d,s), while for 
the e~p interaction, the isotropic distribution is larger [u, c) and the angular 
dependence is weaker {d,s). In addition, both the isotropic distribution in the 
e~^p cross sections and the angular dependence in the e~p cross sections increase 
as X decreases, this is consistent with an increasing sea quark contribution at 
lower X. 



5.2.6 



Measurement and comparison of the dependence 
of NC and CC cross sections in e+j9 and e^p colUsions 



The single differential cross sections da^c^cQ^ / dQ"^ have also been measured by 
HI 0, |0|, ig and by ZEUS pO|, |o|, |06|, The combined e+p NC and 

CC data are compared with the e~p data respectively in Figs j5. 21] and ^.22|. The 



^''In Refs.[ p4||204| , the so-called structure function term (/)cc was introduced, which is related 
to the reduced CC cross section by a weak correction factor a = <^(1 + (5weak)- However, since 
this correction is small with respect to the uncertainty of the measurements, we choose not to 
distinguish them here. 
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Figure 5.18: A measurement of the structure function xF^ by HI j|204|/ as a 
function of x for three value sof at 1500 GeV^, 5000 GeV^, and UOOOGeV^ 
based on the e^p 1994-1997 and e~p 1998-1999 NC data. The inner error bars 
represent the statistical error, and the outer error bars show the total error. The 
SM expectations for the reduced cross sections and xF^ from the HI NLO QCD 
Fit J9l/ are also shown. 
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Figure 5.19: A measurement of the structure function xF^ by ZEUS j20^ as a 
function of for two values ofx at 0.56 and 0.1 based on the e~^p 1994-1997 and 
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the outer error bars show the total error. The SM expectations from CTEQ4 fO|/ 
and MRST are also shown. 
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Figure 5.20: The CC reduced cross sections dec from the combined e~^p 1994- 
2000 data (left) and from the e~p 1998-1999 data (right) shown as a function 
of (1 — y)^ for three x values 0.08, 0.13 and 0.25. The inner (outer) error bars 
represent the statistical (total) errors. The curves are expectations from the NLO 
QCD Fit 
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Figure 5.21: The combined e~^p NC cross sections da/dQ"^ from 1994-2000 are 
compared with the e~p cross sections from 1998-1999 and the corresponding ex- 
pectations from the NLO QCD Fit. The inner error bars represent the statistical 
error, and the outer error bars show the total error. The shaded error bands 
represent the uncertainties on the SM expectations of the e~p and e^p cross 
sections. 
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Figure 5.22: The combined e~^p CC cross sections da/dQ"^ from 1994-2000 are 
compared with the e~p cross section from 1998-1999 and the corresponding ex- 
pectations from the NLO QCD Fit. The inner error bars represent the statistical 
error, and the outer error bars show the total error. The shaded error bands 
represent the uncertainties on the SM expectations of the e~p and e^p cross 
sections. 
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measurement of the NC cross sections spans more than two orders of magnitude in 
. The cross sections fall with by about six orders of magnitude. The e~p NC 
data have a statistical precision varying between 1.3% and 58% from low to high 
and are dominated by the systematic uncertainty for < 1000 GeV^. The 
combined e+p NC data have a better statistical precision between 0.5% to 38%, 
which dominates over the conservatively estimated systematic uncertainty only 
at > 8000 GeV^. The ratio of the measurements with the SM expectations, 
represented in the middle and lower figures respectively for the e~p and the 
combined e'^p data, show that the data have no significant deviation from the 
SM DIS expectations, albeit may have some structure at the very high region 
in particular for the combined e~^p data. The SM uncertainty, which varies ~ 3% 
to ~ 8% from the low to high values shown, represents the uncertainty of the 
expectation due to the assumptions made in the fit, as well as the uncertainties 
of the experimental data entering the fit |94]. 

The e~p CC data have a statistical precision between 7.5% and 19.2% with the 
best precision at ~ 3000 GeV^. The combined high statistics e~^p CC data have 
about the same precision at high and marginally better precision at other 
values due to the smaller e~^p CC cross sections. The cross section uncertainties 
are mostly dominated by the statistical errors. The CC cross sections covering 
about the similar Q'^ range as the NC data shown in Fig. |5.21| only vary over about 
three orders of magnitude. The difference is due to the dominant contribution 
at low from the photon exchange in the NC process. At high Q^, the NC 



cross sections are comparable with the CC cross sections, as shown in Fig. ^.23 



confirming electroweak unification of bosons exchanged in the t-channel. The SM 
uncertainty is larger for the e~^p data, as at high Q'^ the cross section is dominated 
by the d valence quark which is less well constrained than the u valence quark. 



5.2.7 Measurement and comparison of the x dependence 
of NC and CC cross sections in e'^p and e^p collisions 

The single differential cross sections da^ci^cc) / dx have also been measured by 
HI g, |2g and by ZEUS [gg. In Figs.|;2| the e'p NC cross sections from 



1998-1999 are compared with the e'^p data from 1994-1997. In order to see 
the difference due to the change in the lepton-beam charge, a small correction 
is applied to the e^p data to account for the different center-of-mass energies. 
The HI measurement for > lOOO GeV^ (Q^ > lOOOOGeV^) extends in x 
from 0.02 to 0.65 (from 0.13 to 0.65). The cross sections rise towards low x, a 
behavior seen already from the double differential cross sections (Figs.pTTl and 



5.16| ). The decrease of the cross section at low x edges are due to the kinematic 
cut y < 0.9. The SM expectations for e~p and e~^p collisions are compared with 
the data. While only a small difference between the e~p and e~^p data is seen for 
Q2 > 1000 GeV^ it is much more pronounced for > 10 000 GeV^ in accordance 
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Figure 5.23: The NC and CC cross sections measured by HI and ZEUS for 
e'^p (left) and e~p (right) collisions in comparison with the corresponding SM 
expectations from the NLO QCD Fit The inner error bars represent the 

statistical error, and the outer error bars show the total error. 
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Figure 5.24: The e^p NC cross sections da/dx from 1998-1999 in comparison 
with the e~^p data from 1994-1997 and the corresponding expectations from the 
NLO QCD Fit / p^/ . The dotted hnes represent the expectations from the pure 
photon exchange only. The inner error bars represent the statistical error, and 
the outer error bars show the total error. The shaded error bands represent the 
uncertainties on the SM expectations of the e~p and e^p cross sections. 
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with the constructive and destructive interference from the contribution of xFs 
(see Fig. p.l| and Sec. |5.2.4| ). The measured cross sections clearly disagree with the 
expectation for pure photon exchange only, which is independent of the lepton- 
beam charges, as indicated with the dotted curves in Fig. p.24 . The data are in 
good agreement with the SM expectations with a tendency to be loweif^ than the 
expectation for the e~p data at > lOOOOGeV^, as well as at x = 0.65 for all 



data in the full range. The implication will be discussed in Sec. ^.3| . Also shown 
are the SM uncertainties, which for > 1000 GeV^ varies between ~ 2.5% at 



1000 GeV^ varies between ~ 
X = 0.02 to ~ 7% at X = 0.65 and degrades slightly for > 10 000 GeVl 

The new e~p CC cross sections measured by HI using data taken in 1998- 



1999 pOl is compared in Fig.05| with the e+p HI data from 1994-1997 El. As 
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Figure 5.25: The e~p CC cross sections da/dx from 1998-1999 in comparison with 
the e^p data from 1994-1997 and the corresponding expectations from the NLO 
QCD Fit. The inner error bars represent the statistical error, and the outer error 
bars show the total error. The shaded error bands represent the uncertainties on 
the SM expectations of the e~p and e~^p cross sections. 

for NC, a small correction is applied to account for the different proton beam 
energies. The e~^p cross sections from ZEUS [ ^03| | is also shownQ The larger e~p 
cross sections compared to the e~^p cross sections reflect again the different quark 

-'^^The agreement between the e+p data and the expectation from the NLO QCD Fit is better 
since the data were included in the fit. The new e~p data, which indicate some difference with 
the fit, should bring additional constraints on the various parton distribution functions. 

"'^^Since the original cross sections are measured for > 200 GeV , a (large) correction is 
needed, which amounts to 0.299, 0.594, 0.728, 0.783, 0.812, 0.823, and 0.817 respectively for 
X = 0.015, 0.032, 0.068, 0.13, 0.24, 0.42, and 0.65. 
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flavors probed by the bosons in the e^p scatterings and the helicity structure 
of the CC interaction (see Eqs. (|5.11| ) and ( |5.12|) ). Within the experimental errors, 
the data agree well with the SM expectations represented by the NLO QCD 
Fit p3. The uncertainties of the SM expectations are shown by the shaded 



error bands. As mentioned already in Sec. |5.2.6] , the uncertainty, which becomes 
increasingly large towards high x, is significantly larger for the e^p expectation 
than for e~p. 



5.3 Valence quark distribution at high and 
high X 



As it was demonstrated in Fig. p.lT] (Sec. |5.2.3| ), the CC cross sections at high 
and high x are dominated by the up and down valence quarks xd^) 
contributions. These data can therefore be used to constrain these valence quark 
densities at high and high x. The NC cross sections, though they do not 
possess the same quark flavor discriminant power as the CC cross sections, provide 
an additional constraint on xu^ since its contribution to the NC cross section 
dominates over that of xd^ due to the larger quark charge. 

There has been a renewed interest on the behavior of the ratio d/u aX high x 



in recent years ||209| , piO|| . The ratio d/u is relatively well constrained at x < 0.3 
by DIS structure function data as well as the W asymmetry measurement. At 
higher x (0.3 < x < 0.7), the constraint becomes weaker since it is only available 
from DIS structure function data alone. These data are uncertain because of 



• Experimental systematic uncertainties: The fixed-target data are of- 
ten dominated by various systematic errors. It is not at all trivial how these 
errors are treated properly in a QCD analysis. 

• Higher-twist contribution: The fixed-target data are located at rela- 
tively low and high x. The higher-twist contribution is expected to be 
important as it behaves as 1/[(1 — x)Q'^] with respect to the leading-twist 
contribution. 

• Large nuclear corrections: As far as the x valence quark is concerned, 
it is constrained mainly by the deuteron data, of which the nuclear binding 
correction can be large. This explains why the uncertainty of the SM ex- 
pectation on the e^p CC cross section (dominated by the d valence quark 
at high x) was larger than that on the e~p CC and e^p NC cross sections 
(dominated by the u valence quark). 

At even higher x [x > 0.7) there is no reliable data available and the model 
predictions for u/d vary widely from to about 0.2. In most of the conventional 
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QCD analyses performed up to now, the d/u ratio was assumed to approach zero 
as X goes to one. It has been argued recently that a large value around 0.2 is able 
to give a better description of existing data at high x [p09|| . 

When these parton densities are evolved to high Q^, one has to take into 
account additional uncertainties: 

• The uncertainty on the strong coupling constant ag and higher order per- 
turbative corrections. 



• The so called "feed-down" uncertainty [p 1 1|| . The parton density, in the 
convolution on the right-hand side of the DGLAP equations ( |3.42|) and 
( ^.43| ), is sampled at x' > x (Eq. |3.37D . The evolution is thus susceptible 
to the feed-down of an anomalously large contribution at x ^ 1. Such a 
contribution could escape detection by the fixed-target measurements while 
still influencing the evolution to the high region. 

The HERA data at high as constraints on the valence quarks have thus the 
advantage that most of the these potential problems are absent. In the following 
two independent methods will be discussed to provide a determination of the 
valence quark distributions at high . 



5.3.1 Local extraction of valence quark densities 

The fact that the contribution of the valence quark densities dominates in the 
e^p NC and CC cross sections at high x makes it possible to extract xu^ and xd^ 
locally and directly from the measured cross sections. The first extraction was 
performed by HI for two values of x at 0.25 and 0.4 using the e^p NC and CC 
cross sections from the 1994-1997 data: 

xg„(x, Q"^) = ameas(a;, Q^) | ^^^'(^'^ ) j {5.15) 

where cxmeasl^;, Q^) is the measured NC or CC double differential cross sections, 
and the second factor on the right-hand-side of the equation is the theoretical 
expectation from the previous HI fit Only those points where the xq^ con- 
tribution is greater than 70% of the total cross section are selected. The extracted 
parton densities are thus rather independent of the theoretical input as the uncer- 
tainty on the dominant valence quark contribution and that of the corresponding 
cross section largely cancel in the ratio. 

With the new e~p 1998-1999 and e^p 1999-2000 data, similar extractions can 
be made and are extended to x = 0.65 for xm„. In practice, the d valence quark 
densities are determined from the combined e^p CC cross sections discussed in 
Sec. ^.2.2| (Eq. (|5.13|) ). The u valence quark densities are determined respectively 
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from the combined e~^p NC, e p NC, and e p CC cross sections. The three 
independent determinations of then combined. The resulting xUy and 

xdy represent an improved statistical precision of typically 50% and up to 100% 
at high with respect to the first extraction. The results are shown as the data 
points in Fig. ^.26| . 



5.3.2 Valence quark densities from a new NLO QCD Fit 

The valence quark distributions have also been determined from a new NLO QCD 
Fit0 using all cross sections measured by HI at high Q^. The fit was performed 
with the NLO DGLAP evolution equations using the MS renormalization and 
factorization scheme and treating the heavy flavors as massless quarks. Five quark 
and gluon components are parameterized in an MRS-like form at Qq = 15 GeV^: 

xuy = AuX^^{l - x)'^^{l +auVx) (5.16) 

xdy = AdX^''{l-x)^''{l + adVx) (5.17) 

xg = AgX^' {1 - xf' (5.18) 

x(u + c) =Aix^'{l-xf'{l + aix) (5.19) 

x(d + s + b) = A2X^^{l-xf'{l + a2x) (5.20) 

The number of chosen parameters for each component are dictated by the input 
data sets and their precision. In order to have a reliable constraint on the gluon 
and sea quarks, the new low HI data [|135|| (Sec. |4.5| ) are also included. The 



minimum cut on the data is Qmin — 20 GeV^. The parameters Au and A^ are 
constrained by the valence quark counting rules 

1 .1 

Uydx = 2 and / d^dx = 1 . (5-21) 

'0 Jo 

The momentum sum rule allows the determination of one further normalization 
parameter, taken to be Ag. The rest of the parameters are constrained by the fit 
which uses the minuit program |^ to minimize the defined as ||135|| 



(T^ -I- (T^ 

I. Stat ' i,uncor 



(5.22) 



^"Thc fit differs from tlic one in |Q in that the new NC and CC cross sections measured 
from the 1998-1999 e^p and 1999-2000 e+p data at high shown in the previous sections are 
included as well as the new low precision data from 1996-1997 mentioned in Se c .[4.51 , and at 
the same time no fixed targed data are used. The fit also differs from the one in |135f| in that 
the latters uses only the e+p data at \/s = 300 GeV and for < 3000 GeV^ so that the heavy 
flavors could be treated massively. The emphasis also differs, which for this fit is more on the 



valence quarks at high and high x while in | 135 | it is more on the gluon density at low 
and on the extraction of as- 
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where the sums run over the data i, various data sets j and sources k of correlated 
systematic uncertainties. The quantity uj stands for the number of standard de- 
viations for the relative normalization uncertainty of the data set j, and Si^k{sj,k) 
the relative shift of the data i induced by a change by Si^k standard deviations of 
the correlated systematic source k. 

The results of the fit are presented in Table |5.2| in which the is given for 
each data sets, together with their optimal normalization factors. The total 



Table 5.2: Results of the new NLO QCD Fit showing the used data sets 
and their contribution to the a-iid the optimal normalization factors. 



Data set 


e+p NC 


e+p CC 


e-p NC 


e-p CC 


Low 


Total 


data points 


134 


29 


115 


26 


65 


369 




81.9 


19.9 


95.6 


36.4 


55.3 


289.1 


normalization 


1.001 


1.006 


1.010 


1.005 


0.9777 





per degree of freedom (dof) is 289.1/(369-16)=0.82. The normalization factors 
are also well within the quoted luminosity uncertainty. The parameters of the fit 
are given in Table f).3[ 

Table 5.3: Fitted parameters of the quark and 
gluon distributions at Ql = 15 GeV^. The param- 
eters Au, Ad and Ag are constrained by the sum 
rules. 



xq 


A, 


B, 


C 


aq 




1.71 


0.463 


4.93 


12.25 




4.90 


0.872 


5.89 


-0.015 


xg 


3.39 


-0.162 


11.67 




x{u + c) 


0.140 


-0.268 


8.31 


3.54 


x(d + s + b) 


0.140 


-0.268 


13.54 


22.7 



The resulting valence quark distributions xu^ and xd^ are shown in Fig. ^.26 
labeled "NLO QCD Fit: HI only" in comparison with the results obtained from 
the local extraction method (Sec. ^.3.T ), and other parameterizations MRST |58 



CTEQ5 g, and the previous HI NLO QCD Fit g. 
To conclude, the following points are ready to be listed: 

• For the first time, the valence quark distributions at high x are constrained 
by the HI experiment alone with an experimental precision (shaded error 
bands in Fig. ^.26|) varying between 6% at a; = 0.25 and x = 0.4 and ~ 10% 



at x = 0.65 for xu^, and ~ 20% for xdy. The quoted uncertainties correspond 
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Figure 5.26: The valence quarks distributions xuy and xd^ determined both with 
the new NLO QCD Fit (shaded error bands, labeled as 'NLO QCD Fit: HI 
only", see text for the difference with the 'NLO QCD Fit" of g) using all HI 
data at high as well as the new low data JJ35|/ and with the local extraction 
method (data points with the inner and full error bars showing respectively the 
statistical and total errors ) in comparison with other parameterizations which use 
fixed-target data at low Q^. 
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to the experimental errors only. However, it has been studied with a few 
additional fits that other uncertainties are negligible: 



1. The cut has been varied from 20 GeV^ up to 200 GeV^ As ex- 
pected, the valence quarks at high x are constrained by the high 
data only. 

2. The massless treatment of the heavy flavors has been replaced with 
the massive one. The results are stable. 

3. The uncertainty of has been varied by ±0.0017, corresponding to 
the experimental precision of the latest HI determination ||135|| . The 
resulting change on xu^ and xd^ is again small. 

The valence quark distributions determined with the new NLO QCD Fit 
are in good agreement with the other parameterizations at all values of x 
shown except for xu^ at x = 0.65, where the new valence density is about 
~ 17% lower than the other parameterizations with little dependence on 
in the covered range. These other parameterizations all used BCDMS 



data to constrain xu^ at high x. However the discrepancy remains 
small, within about two standard deviations. It does point to a larger ratio 
of d/u, a preferable value as discussed by Yang and Bodek [p09|| . The only 



difference is that the larger ratio is achieved by a smaller xu^ here instead 
of a larger xdy. It is therefore very important that the issue is clarified with 
the future HERA high-precision data at high Q^. 

Finally, the results of the fit agree well with those obtained with the local 
extraction method. 
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5.4 Electroweak tests at HERA 



So far in this report, we have concentrated on the tests and improvements on 
the current knowledge of the strong sector of the SM by imphcitly assuming the 
vahdity of the electroweak sector. The electroweak parameters are fixed to the 
world average values provided mainly by e^e~ high precision experiments. 

In principle, the electroweak sector can also be tested at HERA. A first such 
test was performed by HI with its first 14 CC candidates obtained from 0.35 pb~^ 
of data taken in 1993 [ |188| |. For the first time, the CC cross section was observed 
to be damped at high energy due to the propagator mass of the exchanged W 
boson in the space-like regime (Fig. |5.27 ). 



tot 



1000 




TeV 



Figure 5.27: The energy dependence of the vN cross section. The crosses rep- 



resent the low energy neutrino data 1212^ while the full circle refers to the CC 
result of the 1993 HI data l \18^ , which for the purpose of comparison, has been 
converted to a vN cross section. The experiment at HERA corresponds to an 
equivalent fixed-target energy of about 50 TeV. The full line represents the pre- 
dicted cross section including the W propagator. The dashed line is the linear 
extrapolation from low energies. 
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This first test was made on the total CC cross section. Additional tests are 
now possible with differential CC cross sections. In fact, the double differential 
CC cross section can generally be expressed in terms of a propagator mass Mp^op 
of the exchanged particle: 



decree Gl ( Ml^^^ 



&CC (5.23) 



where acc is the reduced CC cross section defined in Eq.( |5.10| ). The normalization 
of the CC cross section is fixed by the Fermi coupling constant Gp while its 
dependence by Mprop. 

Using the e~^p data taken from 1994 to 1997, which represent an increased 
integrated luminosity by two orders of magnitude with respect to the first CC 
data from 1993, both HI and ZEUS have made a fit of the propagator mass Mprop 
to their cross section data [^, 203]. The propagator mass obtained is respectively 



80.9 ± 3.3(stat) ± 1.7(syst) ± 3.7(pdf) GeV HI (5.24) 
81.4+|^(stat) ±2.0(syst)l|;j(pdf)GeV ZEUS (5.25) 

where the last error is due to the uncertainty of the parton density distributions 
(pdf) dominated by the d valence quark. From this aspect and the fact that the 
e~p CC cross section at high is an order of magnitude larger than the e~^p 
cross section, a better precision is expected from the e~p data. Indeed, with an 
integrated luminosity of 16.4 pb~^ for the 1998-1999 HI e~p data, which is less 
than half of that from the e~^p data, a similar fit yields ||204|| 



79.9 ± 2.2(stat) ± 0.9(syst) ± 2.1(pdf) GeV . (5.26) 

The dominant experimental systematic error is from the hadronic energy scale 
uncertainty, which was estimated to be ±2% for both the e~^p and e~p analy- 



ses (Sec. ^.1.41) . However, the resulting uncertainty on the CC cross sections is 



smaller in e p collisions than in e~^p collisions due to the different dependence 



(Fig. |5.2^ ). In all these measurements, the pdf error is large than the experimen- 
tal systematic uncertainty demonstrating the importance and the necessity to 
improve our knowledge on the parton distribution functions. 

These measurements are in good agreement with the combined results Mw = 
80.427 ± 0.046 GeV and Mw = 80.448 ± 0.062 GeV respectively from LEP II and 



pp colliders [pl4 ]. To achieve such a precision at HERA would be difficult but 



is not impossible at the future high luminosity run |pi5|| . Nevertheless the fact 
that the W mass measured at HERA in the space-like regime agrees so well with 
those when it appears as a real boson with subsequent decays or as a virtual 
boson in the space-like regime indicates that there is little space left for other 
exotic contributions to the measured CC cross sections at HERA. 
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Other tests e.g. on the weak neutral current couphngs of quarks Vu, a^, I'd, 
and ttd are possible but again have to wait for the future high luminosity runs 



with polarized beams ||216 
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5.5 Searches for new physics beyond the Stan- 
dard Model 



From an early e~^p data sample of 1994-1996 corresponding to an integrated lumi- 
nosity of 14.2 pb~^ and 20.1 pb^^ respectively for HI and ZEUS, both experiments 
have reported I]!!], [1^ an excess of events at high Q^{> 15 000 GeV^) and at large 
masses M = a/xs, the center-of-mass energies of the electron-parton collision, 
with respect to the SM DIS expectation for the NC process. HI has also ob- 
served a less significant excess in the CC channel. In 1997, the HERA machine 
was particularly successful and both experiments have collected more e~^p data 
than they had before. Since then, new data from both e~p and e~^p collisions have 
been collected respectively in 1998-1999 and 1999-2000 at a higher center-of-mass 
energy ^/s ~ 320 GeV. These data have been used extensively to search for new 
physics phenomena beyond the SM, of which three examples are briefly described 
in the following subsections. 



5.5.1 Search for leptoquarks at HERA 



The ep collider offers the unique possibility to search for s-channel production 
of new particles which couple to lepton-parton pairs. Leptoquarks are one such 
example. Leptoquarks are color triplet bosons which appear naturally in vari- 
ous unifying theories beyond the SM such as Grand Unified Theories [ pi7|| and 



Superstring inspired Eq models [pl8|] , and in some Compositeness [pl9|] and Tech- 
nicolor ||220|| models. 

Leptoquarks are produced at HERA in the s-channel (Fig. ^.28| (a)). They 
appear also as exchanged bosons in the w-channel (Fig. |5.28] (b)). Leptoquark 




Figure 5.28: Diagrams (a) s-channel resonant production and (b) u-channel ex- 
change of a leptoquark (LQ) with fermion number F = 2. Diagrams involving a 
F = leptoquark are obtained from (a) and (b) by either changing the lepton- 
beam charge or exchanging q and q. 

states are classified according to the fermion number F = L + 3B, where L 
is the lepton number and B is the baryon number of the state. Thus for e~p 
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collisions, the fermion numbers are F = 2 [F = 0) for electron fusing with 
a quark (antiquark) from the proton whereas for e^p collisions they are F = 
[F = —2) for positron fusing with a quark (antiquark) from the proton. The F = 
2 leptoquarks would be produced with a larger cross section in e~p collisions while 
the F = states are best searched for in e~^p collision data. The unknown Yukawa 



coupling of the leptoquark to an eg pair is denoted by A in Fig. |5.28| . Leptoquarks 
coupling only to first generation fermions give e + q ot u + q' final states leading to 
individual events indistinguishable respectively from standard NC and CC DIS 
events, which become backgrounds for the searches considered here. Statistically, 
however, one expects for scalar leptoquarks produced in s-channel a resonant 
peak in the mass distribution and an isotropical da/dy distribution where y = 
(1 + cos^*)/2 with 6* being the decay polar angle of the lepton relative to the 
incident proton in the leptoquark center-of-mass frame. Events resulting from 
the production and decay of vector leptoquarks would be distributed according 
to da/dy oc (1 — y)^. These y spectra from scalar or vector leptoquark production 
are markedly different from the distribution {da/dy oc expected at fixed x for 
the dominant t-channel photon exchange in NC DIS events. This is illustrated in 
Fig.^:29l(a) for NC DIS events selected from e+p 1994-1997 data and in Fig.|539|(b) 
for a production of a scalar leptoquark of mass at 200 GeV, where the kinematic 
variables ye and Mg ■ 



isXe were defined using the electron method (See p. 6. if ). 
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Figure 5.29: Kinematics in the ye — Me plane of (a) the selected NC DIS can- 
didates from the e+p 1994-1997 HI data (two isocurves at = 2500 GeV^ and 
15 000GeV^ are plotted as full lines); (b) a scalar F = leptoquark of mass at 
200 GeV decaying into e + q for a coupling A = 0.05 generated with the LEGO 
event generator 



For this reason, a mass dependent y cut is applied to optimize the signal 
significance and background contribution. The resulting mass distributions are 
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shown in Fig. |5.3CI| for all e^p data taken by HI since 1994. The excess around 
200 GeV in the e+p 1994-1997 data is not confirmed by the new e'p 1998-1999 
and e+p 1999-2000 data. 

Since no significant evidence for a leptoquark signal has been observed, con- 
straints on the Yukawa coupling of leptoquarks were derived |p22| , |223| , |224|| . Two 
examples are shown in Fig. |5.31| for a scalar leptoquark with F = (S'i/2,l) from 
the e~^p 1994-1997 data and a scalar leptoquark with F = 2 (Sq^l) from the 
e~p 1998-1999 data in comparison with the indirect limits from HI p26| , p27 



(Sec. |5.5.^ ) and LEP ||228|| , as well as limits from the Tevatron experiments [|229|| . 
At LEP [ |228| , |230| , |231| , sensitivity to a high-mass leptoquark arises from effects 
of virtual leptoquark exchange on the hadronic cross section. In both examples 
shown, the HI limits extend beyond the present reach of other colliders. More- 
over the constrains obtained from the e~p data on F = 2 leptoquarks are more 
stringent than those from the previous higher statistics e^p data on F = states, 
due to the enhanced energy in the center of mass of the ep collision. Similar limits 
on F = leptoquarks have also been obtained by ZEUS using the e~^p 1994-1997 
data 



5.5.2 Search for supersymmetry at HERA 

Supersymmetry (SUSY) is one of the most likely ingredients for a theory beyond 
the SM. In particular the Minimal Supersymmetric extension of the Standard 
Model (MSSM) describes as well as the SM all experimental data, and in addi- 
tion it offers among its appealing consequences solutions for the cancellation of 
quadratic divergences occurring in the scalar Higgs sector of the SM and models 
beyond the SM. 

SUSY relates fermions to bosons and predicts for each SM particle a partner 
with spin differing by half a unit. For example selectrons ei, are scalar partners 
of electrons e^, cr, and similarly squarks {uLydi), ur, cIr are the partners of up 
and down quarks. Two Higgs doublets with vacuum expectation values V2, Vi are 
necessary to generate masses for up- type quarks (^2) and for down- type quarks 
and charged leptons (vi). The partners of the gauge bosons W^,Z^,'y and the 
two Higgs doublets are called gauginos and higgsinos. They can mix and form 
two charged mass eigenstates xf2 (charginos) and four neutral mass eigenstates 
Xi 2 3 4 (neutralinos) . 

Since supersymmetric particles are not observed at the masses of their SM 
partners, SUSY must be broken. In the MSSM, this breaking is achieved by 
adding extra mass parameters M2 and Mi for the SU{2) and U{1) gauginos. 
Thus the masses of charginos and neutralinos depend on Mi,M2,tan/3 = V2/V1 
and the higgsino mass parameter /i. 

i?-parity {Rp), defined as Rp = (_i)3-B+i+2S^ jg ^ multiplicative quantum 
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Figure 5.30: Mass spectra for NC DIS-like final states selected from tie e~^p 
1994-1997 data (b), e-p 1998-1999 data (b), and e+p 1999-2000 data (c). Tie 
data (symbols) are compared with DIS expectation (histograms) before and after 
a mass dependent y cut designed to maximize the significance of an eventual 
scalar leptoquark signal. The greyed boxes indicate the ±lcr uncertainty of the 
statistical and systematic errors of the NC DIS expectation. 
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Figure 5.31: Exclusion limits at 95% confidence level on the Yukawa coupling A 
as a function of the leptoquark mass for (top) a scalar leptoquark with F = 
and (bottom) a scalar leptoquark with F = 2 described by the Buchmiiller-Riickl- 
Wyler model J22^. Greyed and hatched domains are excluded. The HI indirect 
limits are derived in the analysis of contact interactions described in Sec\5.5.c . 
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number which distinguishes particles {Rp = +1) from SUSY particles {Rp = — 1). 
Here B, L, and S denote respectively baryon number, lepton number, and spin of 
a particle. In SUSY models with i?-parity conservation, supersymmetric particles 
can only be produced in pairs and the lightest supersymmetric particle (LSP), 
which is generally assumed to be Xi? is stable. At HERA the dominant Rp- 
conservation process is the production of a selectron and a squark via a t-channel 
exchange of a neutralino ep — eqX (Fig. ^.32|) . The e and q can then decay 
into any lighter gaugino and their SM partners. The decay involving Xi gives an 
experimentally clean signature of missing transverse energy plus an electron and 
a hadronic system. Such a search has been performed by HI [^33|] using 6.38 pb~^ 



of e'^p data from 1994-1995, and by ZEUS |p34|| using 46.6pb ^ of e~^p data from 
1994-1997. 




IL.R 



Figure 5.32: Selectron-squark production via neutralino exchange and the subse- 
quent decays into the lightest supersymmetric particle Xi- 



The most general SUSY theory which preserves gauge invariance of the SM 
allows, however, for Rp violating [fip) Yukawa couplings A, A', A" between one 
scalar squark or slepton and two SM fermions: 

W^, = KjkUL^Ek + K.kLiQjD, + X'rjj.D^Dk. (5.27) 

where k = 1,2,3 are generation indices, Li{Qi) are the lepton (quark) SU{2)l 
doublet superfields and Ei{Dj, Uj) are the electron (down and up quark) SU{2)l 
singlet superfields. Of particular interest for HERA are the fip terms X'LiQjDk 
as HERA provides both leptonic and baryonic quantum numbers in the initial 
state. The resonant squarks at HERA are thus singly produced (in contrast to 
the i?p-SUSY) in the s-channel (Fig. |5.33D with masses up to the kinematic limit 



163 



of ^/s ^ 300 — 320 GeV. From the theoretical understanding of unification, there 
is no clear preference between i?p-conservation and fip, it is thus mandatory that 
the latter possibility is also sought experimentally [^3^ 
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Figure 5.33: Lowest order s-channel diagrams for first generation squark pro- 
duction at HERA followed by (a), (c) decays and (b), (d) gauge decays. In 
(b) and (d), the emerging neutralino or chargino might subsequently undergo fip 
decays of which examples are shown in the dashed boxes for (b) the xt (d) 
the x^■ 



Compared to the squarks in i?p-SUSY, the squarks in /?p-SUSY can have addi- 
tional decay modes decaying via Yukawa coupling into SM fermions as illustrated 
with two examples in Fig. |5.33| (b),(d). Moreover, the LSP (again assumed to be 
Xi), which is no longer stable, decays via A'^^. into a quark, an antiquark and a 
lepton. 

In cases where both production and decay occur through a A'^^^, the squarks in 
-^p-SUSY behave as scalar leptoquarks (Sec. 5.5.I| ) and the constraints obtained 
on leptoquarks are also applicable for squarks. There is however an interest- 
ing difference, namely the exclusion phase space covered by HERA data can be 
relatively larger for squarks than for leptoquarks. This is because the mass con- 
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straints of 242 GeV (205 GeV) from Tevatron ||229|| correspond to a branching 
ratio Peg of the new particle into e + g of 1(0.5), which can be naturally small 
in ^p-SUSY framework given the competition with gauge decay modes of the 
squarks. 

With e~^p collisions, HERA is best sensitive to couplings A'^^^ among the nine 
possible couplings Xyi^, where mainly m-^ squarks are produced via processes in- 
volving a valence d quark. On the contrary, the e~p data will allow to better 
probe couplings A'^^^^ and squarks. 

Depending on whether the produced squarks undergo a decay or a gauge 
decay, there are many different final state event topologies, e.g. (a) a lepton plus 
a jet, (b) a neutrino plus a jet, (c) a right sign lepton plus multijets, and (d) 
a wrong sign lepton plus multijets. Topologies (a) and (b) are indistinguishable 
from NC and CC DIS events respectively as for leptoquarks. Topology (b) has 
only a low sensitivity with the e+ beam since the produced squarks couples 
to a sea quark u from the proton (Fig. |5.33| (c)), the density of which is small at 
high X. The main SM background for topology (c) is also from NC DIS where 
QCD radiation leads to multijets. Topology (d) has such a striking final state 
that it is essentially background free. 

Under the assumption that only one of the Yukawa couplings A'^^^^^ dominate^, 
mass dependent upper limits on these couplings are derived by combining all 
topologies. Combining all contributing channels with different topologies im- 
proves the sensitivity considerably, up to a factor of ~ 5 at lowest masses, com- 
pared to what would be achieved using only the contribution from topology (a). 
The results p36 , 237 1 are shown in Fig. 5. 34 . The HI limits are compared in 



Fig. ^.34| (a) to the most stringent indirect limits. The production of a m squark 



via a X'lii coupling is very severely constrained by the non-observation of neu- 
trinoless double beta decay [P38|| as shown by the dotted curve. The most sever 
indirect limit on the coupling A'^215 which could allow for the production of squarks 



c, comes from Atomic Parity Violation [235, 239 1 and is indicated by the dashed 



curve in Fig.|5.34|(a). For squark masses below ~ 240 GeV HERA limits signif- 



icantly improve this indirect constraint on A'^21 t>y a factor of up to ~ 3. In 
Ref. [p36| |, further limits are derived both in a constrained version of the MSSM 
by assuming a universal mass parameter rriQ for all sfermions at very high scale 
and in the framework of the Minimal Supergravity Model, a more constrained 
SUSY model. The limits in the latter model extend beyond the domain covered 
by other collider experiments ||240|| especially at large tan (3. 



^-'^This is not unreasonable as in the SM the top quark Yukawa coupling is almost a factor of 
40 larger than the bottom Yukawa coupling. 
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Figure 5.34: (a) Exclusion upper limits at 95% confidence level by HI ||236|| on 
the coupling A'^^^ as a function of the squark mass for tan/3 = 2, in the "phe- 
nomenological" MSSM. For each squark mass, a scan on the MSSM parameters 
M2 and n has been performed and the largest (lowest) value for the coupling limit 
is shown by the upper (lower) curve, (b) The similar limits from ZEUS |P37|| for 
tan (3 = 2 as well as for tan/3 = 10 (dotted lines). 



5.5.3 Search for contact interactions 

A broad range of hypothesized non-SM processes at mass scales beyond the HERA 
center-of-mass energy, ^/s = 300 — 320 GeV, can be approximated in their low- 
energy limit by effective four-fermion eeqq contact interactions, analogous to the 



effective four-fermion interaction describing the weak force at low energies [^41 
Here three sources are briefly discussed, conventional contact interactions arising 
from a substructure of the fermions involved, the exchange of a new heavy particle 
(e.g. leptoquark), low scale quantum gravity effects which may be observable at 
HERA via the exchange of gravitons coupling to SM particles and propagating 
into extra spatial dimensions. 

New currents or heavy bosons may produce indirect effects through the ex- 
change of a virtual particle interfering the 7 and fields of the SM. In the 
present of eeqq contact interactions which couple to a specific quark flavor (g). 
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the SM Lagrangian £sm received the following additional terms ||242|, E3 



a,b=L,R 

(5.28) 

where the sum runs over the left-handed (L) and right-handed (R) fermion he- 
licities. The three terms correspond to respectively the scalar, vector, and tensor 
interactions. Only the vector term is considered at HERA since strong limits 



beyond the HERA sensitivity have already been placed on the other terms p43 
The effective coupling coefficients rjl^^ is defined as 
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g 
Vab 



A^ 



(5.29) 



where e = ±1,0 with ±1 representing constructive and destructive interference 
with the SM currents, g is the overall coupling strength, and A^^ is the effective 
mass scale. In the study of possible fermion compositeness or substructure, g is 
conventionally chosen to be A/ivr, whereas for study of the virtual leptoquark, 
these parameters are related to the leptoquark mass Mlq and the coupling A: 



9_ 
A 



A 



Ml 



LQ 



(5.30) 



Based on the NC cross sections described in Sec. p. 2. 6] , the lower bounds on 
the scale parameters A^ obtained by HI ||226|| using the e'^p 1994-1997 data range 



between 1.3 TeV and 5.5 TeV at 95% confidence level for various chiral structures. 



The limits from ZEUS ||247|| are similar and range from 1.7 TeV to 5 TeV. The most 
restrictive range from both experiments is for the VV model, where all contact 
terms enter with the same sign. Combining with the new e~p 1998-1999 and e'^p 
1999-2000 data, the HI limits have been substantially extended and amount to 
1.6 — 9.2 TeV ||227|| . The HERA results of direct searches for (eg) compositeness 
are thus competitive with those of other experiments at LEP [p4^ 
Tevatron i25( 



2491 and 



The same data have been used by HI to derive lower limits on the ratio 
Mlq/A for the exchange of virtual leptoquarks for range between the center-of- 
mass energy ^/s and up to 1.7 TeV ||227|| . These measurements extend thus the 
direct leptoquark searches at HERA to high masses Mlq > ^/s. Two examples 
have been shown for Si/2,l and Sq^l in Fig. |5.31| . The most stringent limits, 
however, are those for vector leptoquarks with coupling to up quarks. 

It has been recently suggested that gravitational effects may become strong at 
subatomic distances and thus measurable in collider experiments [ p51|] . In such 



a scenario, which may be realized in string theory, gravity is characterized by a 
scale Ms ~ (9 (TeV) in 4 + n dimensions. The extra spatial dimensions n are 
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restricted to a volume associated with the size R and the scale in 4 + n and the 
ordinary 4 dimensions are related by 



~ R'M's^'' , (5.31) 

where Mp ~ 10^^ GeV is the Planck mass. An exciting consequence would be 
a modification of Newton's law at distance r < R, where the gravitational force 
would rise rapidly as F oc l/r^^" and become strong at the scale Ms- Experi- 
mentally, gravity is essentially not tested in the sub-millimeter range ||252|| and 



scenarios with n > 2 extra dimensions at large distances R < 100 (xm are con- 
ceivable. 

In models with large extra dimensions the spin 2 graviton propagates into the 
extra spatial dimensions and appears in the 4-dimensional world as a spectrum 
of massive Kaluza-Klein excitations with masses m^-'^ = j / R, including the zero- 

(i) 

mass state. The graviton fields Gjiu couple to the SM particles via the energy- 
momentum tensor T'^'^ 

^G = -^G!f2T^^ (5.32) 



Summation over the whole tower of Kaluza-Klein states j with masses up to 
the scale Ms compensates the huge 1/Mp suppression and results in an effective 



contact interaction coupling p53 



where A is the coupling strength of order unity. 

Lower limits from HI on the scale parameter Ms are derived again from fits to 
the NC cross sections (Sec. 5. 2.^ ). For the e~p data stronger bounds are obtained 



for positive coupling than for negative coupling. The opposite behavior is ob- 
served in e^p scattering. Both lepton-beam charges thus complement each other 
and a combined analysis ||227|| of all e^p data yields limits on Ms of 0.63 TeV 



for positive coupling A = 1 and 0.93 TeV for negative coupling A = — 1 as illus- 
trated in Fig. p.35| . Similar investigations of virtual graviton exchange in e'^e~ 



annihilation provide comparable limits [248, 254 . 
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Figure 5.35: NC cross sections da/dQ'^ normalized to the SM expectation for 
the combined e^p and the e~p scattering data. The curves show the effects of 
graviton exchange in large extra dimensions given by a common fit to the scale 
Ms (95% confidence level lower limits) with coupling A = 1 and A = — 1. The 
errors represent statistics and uncorrelated experimental systematics. 
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Chapter 6 

Summary and outlook 



The ep collider HERA has opened up a completely new kinematic domain un- 
explored in the previous fixed-target experiments. The results of the first mea- 
surement of the proton structure function at small x down to 0.5 x 10~^ came 
as a surprise for many: the measured x dependence of the structure function for 
a fixed momentum transfer revealed a strong rise towards low x refiecting an 
increasing number of slow partons being probed within the proton. This first 
measurement was based on about 25 nb^^ of data taken in 1992. One year later, 
the data have been increased by more than a factor of 10. The data have allowed 
the structure function measurement to be extended to values of beyond those 
measured at fixed-target experiments. By making use of the hard initial radiative 
events, it was possible to extend the kinematic region to lower and high x. 
With more data and the improved apparatus, the precision measurements are re- 
alized at HERA reaching a statistical precision below 1% and a systematic error 
of typically 3%. 

The structure function measurement of 1993 has allowed a first determination 
of gluon density at low x based on the NLO DGLAP evolution equations. A 
strong-rise behavior toward low x is also found and is believed to be responsible 
for the observed behavior of the structure function. The subsequently improved 
measurements of the structure function have resulted in an ever- increasing knowl- 
edge of the gluon distribution at low x. Today the gluon density is better known 
at low X than at other ranges, a situation completely different to several years ago, 
before the HERA start. The latest precision measurement at low and x by HI 
has not only provided a best determination of the gluon but also together with 
the BCDMS data one of the most precise measurements of the strong coupling 
constant a^. 

The increased integrated luminosity has also allowed to explore the very high 
g2 region up to ~ 30 000GeV^ With the earlier data from 1994 to 1996, an 
excess of events at high and high x was reported by HI and ZEUS. With 
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new data available taken for different lepton-beam charges and higher center-of- 
mass energy, no excess was confirmed. The NC and CC cross sections have been 
measured both double differentially and single differentially in x and Q^. The 
comparison of the e^p NC cross sections has revealed the presence of the 
interference contribution at high and has allowed a first measurement of the 
structure function xF^.i.x^Q'^) at high . The comparison of the e^p CC cross 
sections has shown a huge difference in the cross sections reflecting the different 
quark flavors probed by the bosons in the scatterings. The comparable 
NC and CC cross sections at high has confirmed electroweak unification in 
the space-like regime. 

The measurement of the NC cross sections covers a range of about four 
orders of magnitude, over which the cross sections fall with by more than eight 
orders of magnitude, a behavior well described by the Standard Model expectation 
based mainly on low data providing thereby a nontrivial test of the validation 
the QCD evolution. The Standard Model expectations also describe all other 
measurements with two possible exceptions: 

• The NC cross sections at > 20 000GeV^ and x around 0.4 are higher 
than the expectation in both e~^p data samples taken at different center-of- 
mass energies. The excess remains however statistically not very significant. 

• The NC cross sections at x = 0.65, the largest x value measured, are for e^p 
and e~p at essentially all lower than the expectation which is mainly 
based on the BCDMS data. 

A NLO QCD fit using all high data shows the data at x = 0.65 are about 17% 
lower than the expectation and there remains a discrepancy of about two standard 
deviations. The fit has demonstrated that the quark distributions at high x and 
high are constrained by the HERA data alone with an experimental precision 
of ~ 6% at X = 0.25 and x = 0.4 and ~ 10% at x = 0.65 for the u valence 
quark and about 20% for the d valence quark. The results are in good agreement 
with that obtained by the local extraction method. These determinations are free 
from any nuclear correction in contrast to the current knowledge on the d valence 
quark which was mainly obtained from fixed-nuclear-target data at low Q^. 

The high data have also been used to test the electroweak sector of the 
Standard Model. The derived propagator mass of the W boson in the space-like 
regime agrees well with other direct or indirect measurements of a real or virtual 
boson in the time-like regime. The determination at HERA is however already 
limited by the uncertainty of the parton density distributions demonstrating the 
importance and necessity to improve our knowledge on these parton distributions 
for future precision measurements. 

On the same data, various searches have been performed for seeking both 
resonance production of leptoquarks and squarks in i?-parity violating supersym- 
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metric models, and contact interactions of a substructure of the involved fermions, 
exchange of a new virtual heavy particle, as well as low scale quantum gravity 
effects. No significant deviations from the Standard Model expectations were 
found and various limits have been derived. These limits either extend beyond 
or are comparable with those provided by indirect measurements or by direct 
searches at other colliders. 

At high Q^, the precision of essentially all current measurements are still 
limited by the statistical error. In order to improve the precision, to reach even 
higher and to fully explore the discovery potential of HERA, it is important 
to have as much data as possible and in a relatively short period of time with 
stable running conditions to minimize additional systematic uncertainties. This 
can only be achieved with a luminosity upgrade. The upgrade has started since 
the beginning of September 2000 and the new physics run of the HERA phase two 
is foreseen in September 2001. With this upgrade program for both the machine 
and the detectors together with the polarized beams, HERA, being one of few 
high energy machines running in the next years, will have an exciting and bright 
future. 
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